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of  Tttoodoroon  and  Loevy  type  lift  deficiency  function*.  The  forward -flight 
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There  v*  fiv  aercmechanlcal  iyitm  that  operate  in  a more  hostile 
aerodynamic  environment  than  the  helicopter  rotor.  In  hover,  the  wakes 
shed  by  the  rotor  blades  can  tend  to  pile  up  beneath  the  rotor  and  lead 
to  unfavorable  blade  aeroelastlc  response.  In  forward  flight  the  rotors 
are  required  to  provide  lift  and  propulsive  force  while  the  blades  are 
being  subjected  to  rapidly  varying  aerodynamic  forces.  The  advancing 
blades  can  operate  at  high  tip  Mach  numbers  and  might  experience  adverse 
compressibility  effects,  while  the  retreating  blades  experience  both 
reversed  flow  and  stall  effects . Since  the  rotor  blades  are  relatively 
flexible,  these  aerodynamic  phenomena  can,  under  certain  conditions,  lead 
to  aeroelastlc  Instability.  Typical  of  these  are  advancing  blade  static 
pitch  divergence  and  classical  flutter,  while  the  retreating  blades  can 
encounter  stall  flutter  and  drag  divergence.  Other  phenomena  that  lead 
to  oscillations  of  the  tip  path  plane  are  also  possible.  In  hover,  the 
shed  wakes  can  induce  blade  flutter  and  the  tip  vortices  shed  by  the 
blades  can  cause  subharmonic  blade  oscillation  and  can  also  result  in 
the  blades  going  out  of  track. 

In  addition  to  the  above,  other  instabilities  can  be  encountered  that 
are  a direct  result  of  blade  kinematic  coupling.  Among  these  are  classical 
pitch-flap  and  pitch-lag  instabilities,  caused  by  adverse  coupling  betveen 
the  pitching  and  out-of-plane  and  inplane  motions  of  the  blades.  Also, 
the  fact  that  the  blades  are  relatively  flexible  can  give  rise  to  fre- 
quency relationships  that  result  in  a class  of  instability  generally 
referred  to  as  flap-lag  instability. 

Many  of  the  phenomena  encountered  by  helicopter  rotors  are  little 
affected  by  motions  of  the  rotor  hub.  These  can  generally  be  classed  as 
rotor  instabilities.  However,  a number  of  unstable  phenomena  that  have 
been  experienced  are  a direct  result  of  coupling  betveen  rotor  motions 
and  hub  motions.  The  most  classical  of  these  is  ground  resonance,  or 
more  aptly,  mechanical  instability.  This  can  occur  when  the  frequencies 
and  damping  of  the  modes  of  the  helicopter  on  the  ground  bear  a certain 
relationship  to  the  rotor  blade  inplane  frequency  and  damping,  when  this 
inplane  frequency  is  less  than  the  rotor  speed.  This  instability  is 
purely  mechanical,  that  is,  it  requires  no  aerodynamic  influence.  A 
similar  phenomenon  that  does  require  aerodynamics  occurs  when  the  heli- 
copter is  airborne.  This  is  generally  referred  to  as  air  re.unsnace  and 
normally  involves  the  airframe  rigid- body  modes,  bending  modes,  or  local 
transmission  modes.  When  a rotor  is  in  high-speed  axial  flight,  for 
example,  a propeller  rotor,  an  instability  that  causes  whirling  of  the 
rotor  hub  can  occur.  The  precise  form  of  this  instability  and  the  nature 
of  the  whirling,  advancing  or  regressing,  depend  on  the  degree  of  flexi- 
bility of  the  blades,  or  more  exactly,  the  blade  out-of-plane  and  inplane 
frequencies.  Another  type  of  whirling  instability  that  involves  rotor 
coning  and  blade  flapping  has  recently  been  encountered  during  an  NASA- 
sponsored  test  program  at  Sikorsky  Aircraft. 
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Rotor  blade  pitch  control  systems  have  a wide  variety  of  daalgna  and 
dynamic  characterics . They  invariably  include  servo  actuator*  and  can 
significantly  affect  the  rotor  system  aeroelaatic  response  • The  blade  pitch 
frequency  can  be  substantially  reduced  by  coupling  vith  the  control  system 
which  can  lead  to  an  otherwise  stable  system's  becoming  unstable  if  the 
reduction  in  frequency  results  in  adverse  blade  modal  interactions.  Air- 
ferns  motion  feedback  through  the  servo  supports  or  input  valve  can  also 
affect  the  rotor  response. 

Based  on  the  foregoing , it  seems  clear  that  if  predictions  relating 
to  rotor  stability  are  to  be  believed,  then  the  mathematical  models  of  the 
system  being  analysed  should  include  adequate  descriptions  of  not  only  the 
rotor  blades  and  aerodynamics , but  also  of  control  systems  and  hub  impe- 
dances when  there  is  any  possibility  that  these  may  influence  the  results . 

The  analysis  described  in  this  report  provides  this  capability  by  incor- 
porating realistic  descriptions  of  the  dynamics  of  fully  coupled  rotor/ 
airfTame/control  systems.  It  also  provides  a refined  description  of  the 
aerodynamics  that  includes  stall,  compressibility,  and  unsteady  effects. 

The  analysis  can  be  used  to  study  the  stability  of  main  or  tail  rotors 
or  propellers  in  hover,  forward  flight,  or  pure  axial  flow.  Since  a major 
objective  of  the  contract  under  which  this  work  was  completed  was  to 
provide  the  Army  vith  the  analysis  and  an  operational  computer  program, 
the  method  of  developing  the  analysis,  its  important  features,  and  the 
equations  of  motion  are  given  in  detail  in  Appendixes  A and  B.  A program 
user's  manual  has  been  provided  as  a separate  document. 

Also  included  in  the  report  are  the  results  of  perametrlc  and  sensitivity 
studies  of  semi  articulated  and  hingeless , composite  bearingless  type  tail 
rotors  wing  the  present  analysis . It  should  be  noted  that  not  all  of  the 
options  available  in  the  analysis  were  exercised  during  these  studies. 
Specifically,  blade  motions  were  described  by  only  three  modes,  no  control 
system  dynamics  were  Included,  hub  motions  were  described  by  only  two  modes, 
certain  forward  flight  dynamic  effects  were  not  incorporated,  and  unsteady 
aerodynamics  were  not  used.  These  facts  should  be  taken  into  consideration 
when  conclusions  are  drawn  from  the  results  of  the  studies . 
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Two  baseline  tall  rotor*  that  ara  representative  of  ourrant  Sikorsky 
Aircraft  designs  vara  salactad.  Tha  physical  and  structural  characteris- 
tics of  tha  rotor  blades  and  rotor  hubs  are  presented  respectively  in 
Tables  1 and  2.  Both  rotors  vara  nonart iculated  in  tha  Inplane  direction. 
The  tail  rotor  that  is  articulated  in  the  out-of-plane  direction  vill  be 
referred  to  as  tail  rotor  no.  1,  vhlle  the  tail  rotor  having  a rigid  root 
boundary  condition  for  out-of-plane  notion  vill  be  referred  to  as  tail 
rotor  no.  2.  It  should  be  noted  that  the  hinge  offset  for  out-of-plane 
notion  listed  for  tail  rotor  no.  2 in  Table  1 represents  an  equivalent 
hinge  offset  • or  the  radial  location  at  vhich  the  rotor  blade  flexibility 
can  be  assuned  to  becon*  finite  in  the  flatvlse  direction. 

Additional  rotor  blade  characteristics  are  presented  for  both  rotors 
in  Figures  1 through  6.  These  show  respectively  the  radial  distributions 
of  chord,  tvist,  veight,  flatvlse  and  edgevlse  area  moments  of  inertia, 
and  torsional  inertia.  The  blade  flatvlse  mass  moment  of  inertia  vas 
assumed  to  be  sero,  vhlle  the  edgevlse  mass  moment  of  inertia  vas  assumed 
to  equal  the  torsional  mass  moment  of  inertia.  The  blades  vere  assumed 
to  be  torsionally  Inelastic,  the  torsional  frequency  depending  only  on 
the  control  system  stiffnesses  and  the  blade  torsional  inertia. 

All  conditions  investigated  assumed  sea  level,  standard  day  operating 
conditions  (air  mass  density  vas  .002378  slug/ft3  and  speed  of  sound  vas 
1116  ft /sec). 

Tvo  coupled  flatvlse/edgsvise  blade  modes  vere  used  for  each  rotor. 
The  first  mode,  vhich  vas  primarily  flatvlse  for  both  rotors,  represents 
rigid-body  flapping  for  tall  rotor  no.  1 and  flatvlse  bending  for  tail 
rotor  no.  2.  The  second  mode  vas  mostly  edgevlse  and  described  a bending 
mode  for  both  rotors. 
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TABU  1.  TAIL  ROTOR  BLADE  CHARACTERISTICS 


Fkraittir  Description 

Tail  Rotor 
*>.  1 

Tall  Rotor 
Ho.,2 

Units 

Radius 

10.0 

5.5 

ft 

Hint*  offset  for  out-of- 
plan*  aotlon 

.4160 

.2917* 

ft 

Hiaft>*r  of  blades 

4 

4 

Rotor  sp**d 

700.0 

1214 

rpo 

Tip  loaa  factor 

.97 

.97 

Prolag  angle 

3.33 

0.00 

deg 

Precone  angle 

0.00 

0.00 

deg 

Airfoil  section 

NACA  0012 

NACA  0012 

Young's  Modulus 
Pitch  horn  radial  loca- 

107 

106 

psi 

tion 

Chordirlse  distance  from 
elastic  axis  to  pushrod 
positive  toward  leading 

9.5 

7.375 

in 

edge 

4.5 

5.0 

in 

Pitch  bean  arm  length 

10.64 

8.91 

in 

Critical  pitch  during 
Critical  blade  bending 

6 s 

2. 

% 

■ode  duping 
Blad*  root  flapping 

.5 

1.5 

% 

boundary  condition 
Blad*  root  inplane 

artic. 

rigid 

boundary  condition 
Radial  location  of  inner 

riged 

rigid 

snubber 

Radial  location  of  outer 

13.5 

in 

snubber 

30.0 

in 

Weight  at  blade  pushrod 
Stiffness  of  pitch  bean 

2.1 

.75 

lb 

arm 

Dasplng  associated  with 

15 ,200 

32,600 

lb/in 

pitch  bean  an 
Stiffness  of  actuator 
shaft  for  pufe  moment 
applied  at  pitch  beam 

0.00 

0.00 

lb  sec/in 

end 

Humber  of  blade  modes 

754,000 

935,000 

in-lb/rad 

(flatvise/edgevlse) 
"equivalent  hinge  offset 

2 

IS 

2 

\ 


IS 


«**.«•*" !*■  ■ 
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TABU  2.  TAIL  ROTOR  HUB  CHARACTERISTICS 


Tail  Rotor  Tail  Rotor 


Parameter  Description 

Ho.  1 

Ho.  2 

Units 

Humber  of  modes 

2 

2 

Generalized  mass  of  node  1 

6439 

44.6 

in-lb-sec 

Generalised  mass  of  node  2 

11739 

75.1 

in-lb-sec1 

Frequency  of  mode  1 

1750 

2530 

cpm 

Frequency  of  mode  2 

502 

2160 

cpm 

Damping  of  mode  1 

3.4 

4.0 

% 

Damping  of  mode  2 

4.1 

4.0 

% 

Lateral  displacement  of 
mode  1 

0.00 

0.00 

in 

Lateral  displacement  of 
mods  2 

0.00 

0.00 

in 

Longitudinal  displacement 
of  mode  1 

-37.6 

-12.77 

in 

Longitudinal  displacement 
of  mode  2 

0.00 

0.00 

in 

Vertical  displacement  of 
mode  1 

0.00 

0.00 

in 

Vertical  displacement  of 
mode  2 

37.4 

12.77 

in 

Year  displacement  of  mode  1 

1.00 

1.00 

rad 

Y nr  displacement  of  mode  2 

0.00 

0.00 

rad 

Pitch  displacement  of  mode  1 

0.00 

0.00 

rad 

Pitch  displacement  of  mode  2 

1.00 

1.00 

rad 

BLADE  CHORD, IN. 


BO 
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TAIL  ROTOR  NO. I 


\ 


TAIL  ROTOR  NO.  2 


0 .2  .4  .6  .8  1.0 

RADIAL  DiSTANCE/ROTOR  RADIUS 


Figure  1,  Tail  Fotor  Hade  Chord  Distribution. 
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BLADE  TWIST 


0 2 .4  .6  .8 

RADIAL  DISTANCE /ROTOR  RADIUS 


Figure  2 . Tall  Rotor  Blade  Twist  Distribution. 


BLADE  WEIGHT,  LB/IN 


Figure  3.  Tail  Rotor 


2 4 6 8 


RADIAL  DISTANCE/ ROTOR  RADIUS 


Figure  5.  Tail  Rotor  Blade  Edgewise  Area  Moment  of  Inertia  I'istribution. 


-SEC* 


0 2 .4  .6  .0  1.0 

RAOIAL  DISTANCE/ROTOR  RADIUS 


Figure  6.  Tail  Rotor  Blade  Torsional  Tne'";is.  "istricution. 
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AEROELASTIC  STABILITY  CHARACTERISTICS  IN  HOVER  AMD  FORWARD  FLIGHT 


The  aeroelastic  stability  characteristics  of  the  tvo  baseline  tail 
rotors  were  developed  for  advance  ratios  of  0,  0.25,  and  0.50.  The  blade 
collective  pitch  at  the  75  percent  radial  location  vas  varied  from  -30  to 
30  degrees  for  the  hover  condition,  while  the  forward  flight  range  vas  from 
-20  to  20  degrees.  The  modal  damping  and  frequency  are  presented  in  Figures 
7 through  13  for  tail  rotor  No.  1 and  in  Figures  14  through  20  for  tail 
rotor  No.  2.  For  each  rotor,  the  stability  characteristics  of  the  following 
modes  are  shown  as  a function  of  blade  collective  pitch: 


1.  blade 

2.  blade 

3.  blade 
U.  blade 

5.  blade 

6.  blade 

7 . blade 

8.  blade 

9.  blade 

10.  rotor 

11.  rotor 


flatwise  symmetric  mode 
flatwise  forward  whirl  mode 
flatwise  backward  whirl  mode 
edgewise  symmetric  mode 
edgewise  forward  whirl  mode 
edgewise  backward  whirl  mode 
torsional  symmetric  mode 
torsional  forward  whirl  mode 
torsional  backward  whirl  mode 
hub  pitch  mode 
hub  yaw  mode 


It  should  be  noted  that  the  modal  damping  given  for  the  blade  whirl 
modes  in  this  report  use  a blade  frequency  referred  to  a rotating-axis 
system.  Thus,  the  fixed-axis  forward  whirl  blade  frequency  is  decreased  by 
the  rotor  speod,  while  the  backward  whirl  frequency  is  increased  by  the 
rotor  speed.  In  most  cases,  when  the  blade  flatwise  and  edgewise  whirl 
frequencies  are  transferred  to  a rotating-axis  system,  they  agree  quite 
well  with  cne  values  for  the  symmetric  modes. 


The  ho/er  stability  response  of  the  blade  flatwise  symmetric  and  whirl 
modes  is  similar  for  both  rotors,  as  seen  in  Figures  7 (a)  and  lU  (a). 

These  show  that  the  rotor  system  stability  characteristics  are  essentially 
symmetrical  about  the  zero  blade  collective  pitch  line.  This  result  is 
substantiated  further  by  an  inspection  of  Figures  7 (b)  and  lU  (b)  for  the 
blade  edgewise  modal  damping . It  is  also  noted  that  the  edgewise  response 
of  baseline  tail  rotor  No.  2 is  much  more  sensitive  to  changes  in  collective 
pitch  than  tail  rotor  No.  1.  Both  rotors  show  reductions  in  stability  at 
zero  collective  pitch  and  at  high  values  of  blade  pitch.  Figures  7 (c)  and 
Ik  (c)  Illustrate  that  the  stability  of  the  blade  torsional  modes  is  not 
significantly  affected  by  a large  change  in  blade  collective  pitch. 


The  stability  characteristics  of  the  rotor  hub  pitch  and  yaw  modes  in 
hover  appear  in  Figures  7 (d)  and  lU  (d)  for  tail  rotors  No.  1 and  No.  2 
respectively.  As  seen  from  Table  2,  the  first  hub  mode  represents  a rigid- 
body  rotation  about  the  airframe  yaw  axis  which  results  in  a longitudinal 
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displacement  at  the  rotor  head.  Similarly , the  second  huh  mode  is  rigid- 
body  rotation  about  the  pitch  axis  which  yields  a vertical  displacement  at 
the  rotor  head.  The  damping  of  the  hub  modes  Increases  with  collective 
pitch  up  to  10  degrees.  As  blade  pitch  is  Increased  further,  a reduction 
in  stability  is  evident. 

The  modal  frequency  response  with  rotor  blade  collective  pitch  in  hover 
is  presented  in  Figure  8 for  tail  rotor  No.  1 and  in  Figure  1$  for  tail 
rotor  No.  2.  The  two  rotors  exhibit  similar  behavior,  although  rotor  No.  1 
shows  larger  variations  in  modal  frequency  with  blade  pitch.  It  1b  seen 
that  increasing  or  decreasing  blade  pitch  from  zero  results  in  reduction  in 
the  flatwise  and  edgewise  modal  frequencies.  The  torsional  whirl  mode  of 
tail  rotor  No.  1 at  first  increases  with  blade  pitch  up  to  -20  degrees  and 
then  decreases.  The  blade  torsional  symmetric  frequency  for  both  rotors  is 
not  plotted  in  Figures  8 and  15,  being  much  higher  thar.  all  others:  5.^3 

cycles/rev  for  rotor  No.  1 and  13.85  cycles/rev  for  rotor  No.  2.  It  is 
unaffected  by  pitch  changes  since  it  is  mostly  a function  of  physical  param- 
eters not  associated  with  blade  pitch  or  thrust.  The  frequencies  of  the 
rotor  hub  pitch  and  yaw  modes  do  not  vary  with  blade  pitch.  It  is  noted 
that  the  edgewise  frequency  for  both  rotors  is  above  one  cycle/rev;  such 
rotors  are  classified  as  stiff-inplane  rotors. 

The  aeroelastic  stability  characteristics  of  the  two  baseline  tail 
rotors  in  forward  flight  are  illustrated  in  Figures  9 through  12  for  rotor 
No.  1 and  in  Figures  16  through  19  for  rotor  No.  2.  Advance  ratios  of  .25 
and  .50  were  investigated.  These  advance  ratios  correspond  to  forward 
velccities  of  109  and  217  knots  for  rotor  No.  1 and  10U  and  207  knots  for 
rotor  No.  2.  Only  the  modal  damping  of  the  blade  flatwise  and  edgewise 
modes  is  presented  since  the  blade  torsional  and  rotor  hub  modal  response 
is  not  affected  significantly  by  advance  ratio.  To  facilitate  comparisons 
of  stability  trends  with  advance  ratio,  the  modal  damping  of  the  blade 
flatwise  symmetric  mode  and  of  the  edgewise  forward  whirl  mode  lb  plotted 
in  Figures  13  and  20  for  rotor  No.  1 and  2 respectively.  These  plotB 
indicate  generally  a reduction  in  system  stability  of  both  modes  with 
advance  ratio.  This  effect  is  greatest  near  blade  pitches  of  -10  degrees. 
The  region  of  highest  stability  of  the  blade  flatwise  symmetric  mode  is 
near  a pitch  angle  of  zero  degrees.  The  damping  decreases  as  advance  ratio 
increases  for  both  tail  rotors.  A similar  behavior  is  exhibited  by  the 
flatwise  backward  and  forward  whirl  modes. 

A comparison  of  the  modal  frequency  response  illustrated  in  Figures  8, 
10,  and  12  for  rotor  No.  1 and  in  Figures  15,  17  and  19  for  rotor  No.  2 
reveals  that  the  effect  of  advance  ratio  on  frequency  is  small  for  all 
modes  except  the  blade  torsional  whirl  mode,  which  shows  a slight  decrease 
in  frequency  with  advance  ratio. 
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Figure  7 . Modal  Damping  Response  With  Blade  Collective  Pitch  for  Tail 
Rotor  No.  1 in  Hover. 
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Figure  8.  Modal  Frequency  Response  With.  Blade  Collective 
Pitch  for  Tail  Rotor  No.  1 in  Hover. 


Figure  10.  Modal  Frequency  Response  Vith  Blade  Collective 
Pitch  for  Tail  Rotor  No.  1 at  p = .25. 
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(b)  Blada  EdgwiM  Forward  Whirl  Noda 


Figure  13.  Effect  of  Advance  Ratio  on  Modal  Damping  for  Tail  Rotor  Ho.  1. 
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I LA  DC  COLLECTIVE  PITCH , OS 


(a)  Uadi  Flatwisa  Nods 


BLADE  COLLECTIVE  PITCH  . DEG 

(b)  Bladt  Edgawita  Noda 


Figure  ll».  Modal  Damping  Response  With  Blade  Collective 
Pitch  for  Tail  Rotor  No.  2 in  Hover. 
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BLADE  COLLECTIVE  PITCH,  OEG 

(*)  Bl«da  FUtwIia  Noda 


BLADE  COLLECTIVE  PITCH,  DEC 

(b)  Blada  Edgawlta  Noda 


Figure  16 . Modal  Jamping  lesponse  With  Blade  Collective 
Pitch  for  Tail  Rotor  No.  2 at  u » .25. 


Modal  Frequency  Response  With  Blade  Collective 
Pitch  for  Tail  Rotor  Ho.  2 at  y = .25. 
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(b)  llada  EdgawlM  Ptodt 


Figure  18.  Modal  Damping  Response  With  Blade  Collective 
Pitch  for  Tail  Rotor  No.  2 at  p ■ .50. 
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BLADE  COLLECTIVE  BITCH,  OEG 


(a)  Bltd*  Fit  twit*  SywatHc  Mod* 


BLADE  COLLECTIVE  PITCH, OEG 


(b)  Blada  Edgwls*  Forwtrd  Whirl  Nod* 


Figure  20.  Effect  of  Advance  Ratio  on  Modal  Damping  for  Tail  Rotor  No.  2 
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The  effects  of  independent  variations  of  several  design  parameters  on 
the  stability  characteristics  of  the  tvo  baseline  tail  rotors  were  investi- 
gated in  hover  and  at  an  advance  ratio  of  .$0.  The  blade  collective  pitch 
vas  8 degrees  throughout  the  study.  The  following  parameters  were  examined: 

1.  control  system  asymmetric  stiffness,  K 

2.  control  system  symmetric  stiffness, 

3.  blade  torsional  natural  frequency,  u^2 

U.  pitch- flap  coupling, 

5.  blade  edgewise  natural  frequency, 

6.  rotor  hub  frequencies  (pitch  and  y&v),  Up  u>Y 

7.  blade  precone  angle, 

8.  pitch-lag  coupling, 

The  results  of  the  parametric  study  are  presented  in  a nondlmenslonal 
form  to  enhance  their  general  applicability.  All  system  frequencies  are 
nondimenslonalized  by  the  rotor  speed. 

The  variations  in  control  system  asymmetric  stiffness  are  illustrated 
in  terms  of  an  equivalent  blade  torsional  frequency  defined  by 

“ti  * (D 


where 


L2  * chordwise  distance  from  elastic  axis  to  the  pushrod,  positive 
toward  leading  edge,  in. 

? 

1 0 * blade  torsional  mass  moment  of  inertia,  in-lb-sec 
K3  * control  system  asymmetric  stiffness,  lb/in. 

The  control  system  asymmetric  stiffness  is  defined  by 


where 


K3  * 2K2Ki/(2K 2 + 


K2  « stiffness  of  actuator  shaft  for  pure  moment  applied  at  pitch 
beam  end,  in-lb/rad 
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K.  ■ control  systoi  symmetric  stiffness  (or  pitch.  bear.  arm 
stlffnsss,  lb/in. 

V ■ number  of  rotor  blsdts 

L ■ pitch  boss  in  length,  In. 

Ths  rsrlsbls  in  Equation  2 la  the  actuator  abaft  moment  stiffness  K . The 
baseline  values  for  are  2.32  cyeles/rev  for  tall  rotor  No.  1 and  5.U2 

cycles/r ev  for  tall  rotor  No.  2. 

The  parametric  trends  obtained  for  the  control  system  symmetric  stiff- 
ness are  also  expressed  in  terms  of  an  equivalent  blade  torsional  frequency, 
which  is  defined  by 

"T2  * L2KJ/Ie  (3) 

Here  the  variable  is  the  pitch  beam  arm  stiffness  . The  baseline  values 
of  are  5.1»8  cycles /rev  for  tall  rotor  No.  1 and  13.85  cycles/rev  for 

tall  rotor  No.  2. 

The  blade  torsional  natural  frequency  is  derived  from  Equation  3 when 
ynmetric  mode  stability  characteristics  are  presented  and  from  Equation  1 
when  unsymmetric  mode  stability  characteristics  are  presented.  Thus,  a 
comparison  can  be  made  of  the  parametric  trends  developed  as  a function  of 
the  blade  torsional  frequency  due  to  Independent  variations  in  stiffness 
or  Inertia. 

The  rotor  blade  pitch-flap  coupling  Is  defined  below  for  tail  rotor 
No.  1,  which  is  articulated  In  the  flatwise  direction. 

TAN  «,  - (L  -e)/L2  (U) 

where  3 3 

L3  * pitch  horn  radial  location.  In. 

e » blade  flapping  hinge  offset.  In. 

For  the  hingeless  tall  rotor  No.  2 an  effective  pitch- flap  coupling  can  be 
evaluated  from 


TM  S3 


B/L2[*F1  ♦ (♦F1-*P2)(r1-L3)/(r2-r1)] 


(5) 
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where 


R ■ blade  radius,  in. 

rj  ■ radial  location  of  inner  snubber,  in. 

T2  ■ radial  location  of  outer  snubber,  in. 

■ flatwise  node  shape  at  r} 

■ flatwise  node  shape  at  r2 

Variations  in  the  blade  edgewise  natural  frequency  were  made  by 
changing  the  radial  distribution  of  the  blade  edgewise  area  moment  of 
inertia  which  appears  in  Figure  5 • 

The  values  of  the  rotor  hub  pitch  and  yaw  frequencies  that  are  referred 
to  as  reference  hub  frequencies  are  defined  in  terms  of  the  generalized 
stiffness  Kq  i and  mass  Mg  ^ of  the  1th  fixed  system  mode  as  follows : 


“i  ■ ’'“a/Vi  <6) 

Variations  in  the  fixed  system  reference  frequencies  were  made  by 
changing  the  generalized  stiffness.  The  reference  values  were  varied 
simultaneously  by  similar  amounts. 

In  summaxy,  the  variables  associated  with  each  of  the  eight  parameters 
investigated  are: 

1.  actuator  shaft  moment  stiffness,  K2 

2.  pitch  beam  stiffness,  Kj 

3.  blade  torsional  mass  moment  of  inertia,  I# 

4.  chordwise  distance  from  elastic  axis  to  puahrod,L2 

5.  blade  edgewise  area  moment  of  Inertia,  1^ 

6.  rotor  hub  frequencies.  Up  and  u>y 

7.  blade  precone  angle,  0q 

8.  pitch-lag  coupling,  a2 


The  baseline  values  of  all  parameters  investigated  are  listed  in 
Table  3 for  both  tail  rotors.  All  baseline  values  apply  at  advance  ratios 
of  zero  and  0.50. 
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TABLE  3.  TAIL  ROTOR  BAS  ELUTE  PARAMETERS  FOR  PARAMETRIC  STUDY 


Parameter  Value 


Number 

Parameter 

Symbol 

Tail  Rotor 
No.  1 

Tail  Rotor 
Ho.  2 

Units 

1 

*3 

2725 

4990 

lb/in. 

2 

K| 

15200 

32600 

lb/in. 

3a 

“Tl/fl 

2.32 

5.42 

cycles /rev 

3b 

“T2/Q 

5.48 

13.85 

cycles /rev 

4 

*3 

45. 

40.5 

deg 

5 

"W“ 

2.06 

1.74 

cycles /rev 

6a 

Wp/Q 

0.72 

1.78 

cycles /rev 

6b 

Wy/tJ 

2.50 

2.08 

cycles /rev 

7 

*0 

0.0 

0.0 

deg 

8 

<*1 

0.0 

0.0 

deg 

Effect  of  Control  8yatcm  Asymmetric  Stiffness 

The  effect  of  control  system  asymmetric  stiffness  on  the  aeroelastic 
stability  characteristics  of  the  baseline  tail  rotors  is  shown  in  Figures 
21  through  26  at  advance  ratios  of  zero  and  0.50. 

Only  the  blade  vhirl  or  unsymmctrlc  modes  are  presented  in  these 
figures,  since  the  symmetric  blade  mode  equations  are  not  affected  by 
changes  in  the  actuator  shaft  moment  stiffness.  Control  system  asymnetric 
stiffness  increases  with  actuator  shaft  moment  stiffness.  This  is  evident 
from  Equation  2,  when  expressed  in  the  following  form: 

K3  - Ki/(1+NKjLi/2K2)  (2a) 

The  blade  torsional  frequency  is  proportional  to  the  control 
system  asymmetric  stiffness,  as  can  be  seen  from  Equation  1. 

Variations  in  actuator  shaft  moment  stiffness  resulted  in  a range  of 
blade  torsional  frequencies  of  1.25  to  4.37  cycles/rev  for  tail  rotor  No. 

1 and  from  2.88  to  10.65  cycles/rev  for  tail  rotor  No.  2.  Four  values  of 
stiffness  were  used  to  generate  the  curves . The  baseline  value  of  the 
parameter  being  investigated  is  indicated  by  an  arrow  along  the  abscissa 
of  all  plots. 
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From  an  Inspection  of  Figure • 21  and  22,  it  is  observed  that  as  the 
moment  stiffness  of  the  actuator  shaft  increases,  the  hover  stability  of 
the  blade  flatvise  whirl  node  increases,  while  the  nodal  damping  associated 
with  the  blade  torsional  and  edgewise  whirl  nodes  and  the  rotor  hub  pitch 
and  yaw  modes  decreases  for  both  baseline  tail  rotors.  However,  it  appears 
that  tho  effect  of  variations  in  actuator  shaft  moment  stiffness  on  rotor 
stability  is  more  pronounced  for  * ail  rotor  No.  1,  especially  for  the  blade 
flatwise  whirl  mode.  Thi«  mode  Incomes  unstable  as  the  blade  torsional 
frequency  approaches  1.25  cycles /rev. 

The  trends  at  an  advance  ratio  of  0.50,  shown  in  Figures  23  and  2k , 
are  basically  the  sane  as  in  hover  except  that  for  tail  rotor  No.  1 at  the 
lower  blade  torsional  frequencies  a strong  interaction  seems  to  exist 
between  the  Dlade  edgewise  and  torsional  backward  whirl  modeo.  Hie  edge- 
wise mode  exhibits  a large  increase  in  modal  damping,  while  the  torsional 
response  becomes  less  stable.  Also,  the  hub  yaw  mode  for  rotor  No.  1 
shows  a slight  increase  in  stability  with  stiffness  at  an  advance  ratio 
of  0.50,  while  it  remained  largely  unaffected  in  the  hover  condition. 

The  effect  of  control  system  asymmetric  stiffness  on  hover  modal 
frequency  is  illustrated  in  Figures  25  and  26  for  tall  rotors  No.  1 and 
2 respectively.  It  is  observed  that  the  modal  frequencies  of  the  blade 
flatwise  and  edgewise  whirl  modes,  as  well  as  the  hub  modes,  are  not 
affected  by  large  variations  in  control  system  asymmetric  stiffness.  As 
might  be  expected,  the  blade  torsional  whirl  frequency  for  both  rotors 
has  nearly  the  same  value  as  the  equivalent  torsional  frequency  o>T  , 
defined  by  Equation  1.  1 

The  parametric  trends  of  modal  frequency  predicted  for  all  parameters 
investigated  in  this  study  for  both  baseline  rotors  were  not  affected  by 
a change  in  advance  ratio  from  zero  to  0.50.  Thus,  only  the  frequency 
results  in  hover  are  shown. 
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Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Actuator  Shaft  Moment  Stiffness)  on  Modal  Lamping  for 
Tail  Rotor  No . 2 in  Hover. 
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Effect  of  Control  By  tea  By— trie  Stlffneas 

Variations  in  the  control  system  symmetric  stiffness , otherwise 
referred  to  as  the  pitch  beam  an  stiffness,  yielded  the  stability  charac- 
teristics presented  in  Figures  27  through  32.  The  blade  torsional  frequency 
Wj2 * defined  by  Equation  3*  ranged  fro*  2.74  to  15.50  cycles/rev  for  tail 

rotor  Vo.  1 and  fro*  6.90  to  39*20  cycles /rer  for  tall  rotor  Ho.  2.  Four 
values  of  stiffness  vere  used  to  generate  the  curves.  Fro*  Equation  2,  it 
can  be  seen  that  a change  in  control  system  symmetric  stiffness  Kj  results 
in  a nav  value  of  control  system  asymmetric  stiffness  K3  unless  the 
actuator  shaft  *o*ect  stiffness  K2  is  varied  accordingly.  This  vas  done 
in  order  to  separate  the  effect  of  the  control  system  asymnetric  stiffness 
discussed  in  the  previous  section  fro*  that  of  the  control  system  symmetric 
stiffness.  Urns,  the  modal  doping  and  frequency  for  the  blade  whirl  nodes 
and  the  rotor  hub  modes  remain  unchanged. 

Hie  modal  damping  and  frequency  of  the  blade  symmetric  modes  are  shown 
as  a function  of  pitch  beam  atlffhess  or  equivalent  blade  torsional  fre- 
quency, «T2,  in  Figures  27  through  32.  An  Inspection  of  Figures  27 

through  30  reveals  that  the  doping  of  the  blade  torsional  mode  decreases 
while  the  damping  of  the  flatwise  mode  increases  with  torsional  frequency. 
Hie  effect  on  the  edgewise  mode  is  very  slight.  Hiese  results  apply  to 
both  baseline  rotors  at  advance  ratios  of  zero  and  0.50.  Figures  31  and 
32  show  that  the  flatwise  and  edgewise  blade  symmetric  frequencies  remain 
constant.  The  blade  symmetric  torsional  frequency  is  almost  the  same  as 
that  defined  by  Equation  3.  Hie  effect  of  variations  of  pitch  beam  stiff- 
ness on  the  blade  sysnetric  modes  is  generally  similar  to  the  effect  of 
actuator  shaft  moment  stiffness  variations  on  the  blade  whirl  or  unsym- 
metric  inodes. 
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Figure  28-  Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Pitch  Beam  Arm  Stiffness ) on  Modal  Damping  for  Tail 
Rotor  No.  2 in  Hover. 
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Figure  31 . Effect  of  Blade  Torsional  Frequency  (*s  a Function  of 
Pitch.  Beam  Am  Stiffness)  on  Modal  Frequency  for  Tail 
Rotor  No.  1 in  Hover. 
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Figure  32.  Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Pitch  Beam  Ana  Stiffness)  on  Modal  Frequency  for  Tail 
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Effect  of  Blade  Torsional  Natural  Frequency 

The  effect  of  changes  In  bind*  torsional  moment  of  Inertia  on  the 
nodal  dating  and  frequency  of  tb*  two  baaelln*  tall  rotor*  Is  Illustrated 
In  Figures  33  through  38.  All  blade  nodes  are  shown  In  these  figures 
sine*  the  blade  torsional  moment  of  inertia,  Ig,  appears  In  both  Equations 

1 and  3,  which  affect  respectively  the  blade  uns /metric  and  symmetric 
degrees  of  freedon. 

The  blade  torsional  frequency,  «T2,  defined  by  Equation  3,  is  used 

to  illustrate  the  results  for  the  blade  symmetric  modes.  This  frequency 
rsnged  from  2.11  to  7.60  cycles /rer  for  tall  rotor  No.  1 and  from  b.9b  to 
17.60  cycles/rev  for  tail  rotor  No.  2.  Six  values  of  frequency  were  used 
to  generate  the  curves. 

The  modal  damping  and  frequencies  for  the  blade  whirl  modes  and  the 
rotor  hub  modes  are  presented  as  a function  of  the  blade  torsional  frequency, 

as  calculated  from  Equation  1.  This  frequency  ranged  from  O.89  to 

3.22  cycles/rev  for  tail  rotor  No.  1 and  from  2.98  to  6.90  cycles/rev  for 
tall  rotor  No.  2.  Six  values  of  frequency  were  used  to  generate  the  curves. 

A comparison  of  the  stability  characteristics  presented  in  this 
section  for  the  blade  symmetric  modes  can  be  made  with  the  results  given  in 
Figures  27  through  32,  which  show  the  effect  of  the  pitch  beam  arm  stiffness. 
Similarly,  the  behavior  of  the  blade  whirl  and  hub  modes  can  be  compared 
with  those  in  Figures  21  through  26,  which  show  the  effect  of  actuator 
shaft  moment  stiffness. 

The  modal  damping  of  the  blade  symmetric  modes  in  hover  appears  in 
Figures  33(a)  and  3b(a)  for  tall  rotors  Nos.  1 and  2 respectively.  It  is 
noted  that  in  the  cannon  frequency  range,  2.7b  to  7.60  cycles/rev  for  tail 
rotor  No.  1 and  6.90  to  17.60  cycles /rev  for  tail  rotor  No.  2,  the  para- 
metric trends  for  both  rotors  are  similar  for  the  blade  torsional  and 
edgewise  modes.  The  flatwise  mode  of  tail  rotor  No.  1 shows  a slight 
decrease  in  stability  with  decreasing  torsional  inertia.  Figure  33(a), 
while  the  same  mode  becomes  more  stable  (Figure  27)  as  pitch  beam  stiff- 
ness increases.  Rapid  changes  in  stability  are  predicted  at  low  torsional 
frequencies,  between  2 and  3 cycles/rev,  for  tall  rotor  No.  1 in  Figure 
33(a).  The  blade  torsional  mode  shown  an  abrupt  reduction  in  stability, 
while  the  edgewise  and  flatwise  modes  exhibit  large  Increases  in  damping. 

This  result  may  be  explained  by  the  modal  interaction  indicated  in  Figure 
37(e)  for  the  edgewise  and  torsional  modes  as  they  both  approach  a value 
of  2 cycles/rev.  For  tail  rotor  No.  2,  the  torsional  frequency  is  well 
removed  from  the  edgewise  and  flatwise  frequencies,  as  shown  in  Figure 
38(a). 
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The  stability  characteristics  of  the  blade  symmetric  modes  at  an 
advance  ratio  of  0.50  art  illustrated  in  figures  35(a)  and  36(a)  for  tall 
rotors  Vos.  1 and  2 respectively.  The  parametric  trends  shown  In  these 
figures  are  similar  to  those  already  discussed  for  the  hover  condition. 

The  modal  frequencies  of  the  blade  symmetric  flatwise  and  edgewise 
modeii  are  not  affected  by  changes  in  blade  torsional  inertia,  as  shown  in 
figures  37(a)  and  36(a)  for  both  rotors.  The  blade  symMtric  torsional 
frequency  is  practically  the  same  as  the  equivalent  torsional  frequency, 
u>j2 . These  results  are  similar  to  those  presented  in  figures  31  and  32 

for  pitch  beam  arm  stiffness  variations . 


The  variations  in  modal  damping  with  blade  torsional  inertia  for  the 
blade  unsymmetrlc  modes  and  the  rotor  hub  modes  are  presented  in  Figure 
33(b)  for  tall  rotor  No.  1 and  in  Figure  34(b)  for  tail  rotor  No.  2 for 
the  hover  condition.  These  results  can  be  compared  with  those  from 
Figures  21  and  22  in  the  common  frequency  range  of  1.25  to  3.22  cycles/rev 
for  tall  rotor  No.  1 and  2.96  to  6.90  cycles /rev  for  tail  rotor  No.  2. 

Good  agreement  exists  in  the  trends  shown  by  ell  blade  and  hub  inodes  for 
toil  rotor  No.  1 and  by  the  blade  edgewise  and  torsional  modes  for  tail 
rotor  No.  2.  However,  tail  rotor  Ho.  2 indicates  a slight  decrease  in 
the  stability  of  the  blade  flatwise  whirl  rode  as  the  blade  torsional 
moment  of  inertia  decreases  (Figure  34(b)),  while  in  Figure  22,  the  damping 
of  this  mode  increases  with  actuator  shaft  moment  stiffness.  For  this 
rotor,  it  is  also  noted  that  the  behavior  of  the  rotor  hub  modes,  especially 
the  yaw  mode,  is  quite  different  for  the  two  parameter  variations  as  the 
blade  torsional  frequency  is  reduced  below  5 cycles/rev. 

As  advance  ratio  is  increased  to  0.50,  the  stability  trends  predicted 
for  the  blade  whirl  and  rotor  hub  modes  as  a function  of  blade  torsional 
inertia  are  generally  the  seme  as  in  the  hover  condition  for  both  tail 
rotors.  These  results  are  illustrated  in  Figures  35(b)  and  36(b)  for  tail 
rotor  No.  1 and  2 respectively.  However,  some  differences  in  behavior 
exist  in  the  blade  edgewise  mode  for  tail  rotor  No.  1 and  the  blade  flat- 
wise forward  whirl  mode  for  tall  rotor  No.  2.  These  differences  probably 
result  from  the  coupling  of  the  symmetric  and  unaymmetric  blade  modes  which 
exists  in  forward  flight. 


The  variations  in  modal  frequency  with  blade  torsional  frequency  for 
the  blade  whirl  and  rotor  hub  modes  are  shown  in  Figure  37(b)  for  tail 
rotor  No.  1 and  in  Figure  38(b)  for  tail  rotor  No.  2 in  hover.  The  results 
for  tall  rotor  No.  1 indicate  a strong  interaction  betveen  the  torsional 
and  edgewise  whirl  modes  as  the  blade  torsional  frequency,  w^/fi  , approaches 

2 cycles/rev.  A further  increase  in  blade  torsional  moment  of  inertia, 

which  lovers  w /ft,  results  in  another  interaction  between  the  torsional 
Ti 

and  flatwise  whirl  modes  near  a value  of  co^/ft  of  1.70  cycles /rev.  It  is 

noted  that  for  tail  rotor  No.  2,  no  such  Interactions  in  blade  whirl 
frequencies  occur  since  the  blade  torsional  whirl  frequency  is  well 
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restored  from  *11  other  system  frequencies  in  the  range  of  investigated. 

Both  rotors  exhibit  negligible  changes  in  the  hub  pitch  and  yaw  modal 
frequencies.  A comparison  of  the  frequency  response  as  a function  of  blade 
torsional  moment  of  inertia,  Figures  37(b)  and  38(b) , or  actuator  shaft 
moment  stiffness.  Figures  25  and  26,  shows  similar  trends  for  both  tail 
rotors.  However,  the  interactions  between  the  blade  whirl  modes  for  tail 
rotor  Ho.  1 are  mure  pronounced  in  Figure  37(b)  than  in  Figure  25. 
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Figure  33*  Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Damping  for  Tail 
Rotor  No.  1 in  Hover. 
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Figure  34 . 


BLADE  TORSIONAL  FREQUENCY  / ROTOR  SPEED,  wTt  /O 


(a)  Blade  Sjaaetrlc  Nodes 


(b)  Blade  Milrl  and  Hub  Nodes 


Effect  of  Blade  Torsional  Frequency  (.as  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Damping  for  Tail 
Rotor  No.  2 in  Hover. 


63 


BLADE  TORSIONAL  FREQUENCY/ ROTOR  SPEED  , w Tt/Q 


(a)  Had*  Symmetric  Nodes 


BLADE  TORSIONAL  FREQUENCY /ROTOR  SPEEO  , wT|  /fl 


(b)  Blada  Whirl  and  Hub  Nodes 


Figure  35.  Effect  of  Blade  Torsional  Frequency  fas  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Damping  for  Tail 
Rotor  No.  1 at  y ■ .50. 


Figure  36. 
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(a)  Blade  Sjmetrlc  Nodes 


BLADE  TORSIONAL  FREQUENCY  t ROTOR  SPEED  , «T,/n 


(b)  Blade  Whirl  and  Hub  Nodas 


Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Damping  for  Tail 
Rotor  No.  2 at  p ■ .50. 


65 


■ » «-. ■, -r-—  — -r**“r*r* '• ' t-rr -tJwjsrAn.*  r — JHn'.ia^' 


0 2 4 8 1C 

BLADE  TORSIONAL  FREQUENCY  / ROTOR  SPEED  , WTI/H 


(«)  Blade  Syawtrlc  Node 


8 

t 


s 


(b)  Blade  Whirl  and  Hub  Nodes 


Figure  37 . Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Frequency  for  Tail 
Rotor  No.  1 in  Hover. 
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BLAOE  TORSIONAL  FREQUENCY /ROTOR  SPEED  , UJt /fl 


(•)  Blade  Sywetrlc  Nodes 


(b)  Blade  Whirl  and  Hub  Nodes 


. Effect  of  Blade  Torsional  Frequency  (as  a Function  of 
Blade  Torsional  Inertia)  on  Modal  Frequency  for  Tail 
Rotor  No.  2 in  Hover. 


67 


Effect  of  Pitch-flap  Coupling 

The  effect  of  blade  pitch-flap  coupling  on  the  stability  characteristics 
of  the  tvo  baseline  tail  rotors  is  illustrated  in  Figures  39  through  M for 
advance  ratios  of  zero  and  0,50.  The  range  in  pitch-flap  coupling  investi- 
gated vas  from  -30  to  60  degrees  for  both  tail  rotors . Five  values  vere 
used.  Equation  k was  used  to  calculate  the  coupling  for  tail  rotor  No.  1, 
which  is  articulated  in  the  flatwise  direction,  while  an  effective  value 
of  pitch-flap  coupling  was  computed  from  Equation  5 for  tail  rotor  No.  2, 
which  is  nonarticulated.  The  stability  characteristics  in  hover  are 
discussed  first. 

From  Figures  39  and  40,  it  is  noted  that  pitch- flap  coupling  stabilizes 
the  blade  torsional  symmetric  and  unsymmetric  modes  and  the  blade  edgewise 
whirl  modes;  this  effect  becomes  more  pronounced  as  coupling  is  increased 
above  30  degrees  for  both  rotors.  A strong  destabilisation  of  the  blade 
flatwise  modes  is  evident  throughout  the  entire  range  of  pitch-flap 
coupling  investigated.  An  instability  of  the  blade  flatwise  whirl  mode 
of  tail  rotor  No.  1 is  predicted  at  a pitch-flap  coupling  of  60  degrees 
(Figure  39).  The  blade  edgewise  symmetric  mode  for  tail  rotor  No.  1 
becomes  marginally  stable  as  pitch- flap  coupling  approaches  60  degrees, 
while  no  change  in  stability  is  shown  by  the  same  mode  for  tail  rotor  No.  2. 
From  the  results  presented  in  Figure  39  and  bO,  it  can  be  concluded  that 
the  stability  of  the  blade  torsional  and  flatwise  inodes  is  considerably 
more  sensitive  to  variations  in  blade  pitch-flap  coupling  for  tall  rotor 
No.  1 than  for  tail  rotor  No.  2. 

The  modal  danping  of  the  rotor  hub  modes  for  tail  rotor  No.  1 is 
generally  unaffected  by  pitch- flap  coupling  variations.  On  the  other  hand, 
tall  rotor  No.  2 exhibits  a slight  reduction  in  the  damping  of  the  rotor 
pitch  mode,  while  the  yaw  mode  shows  a large  increase  in  stability  at  the 
higher  values  of  pitch- flap  coupling . 

The  parametric  trend  of  modal  damping  with  pitch-flap  coupling  at  an 
advance  ratio  of  0.50  is  illustrated  in  Figures  1*1  and  1*2  for  tail  rotors 
No.  1 and  2 respectively.  A comparison  of  these  results  with  the  results 
presented  in  Figures  39  and  U0  for  the  hover  condition  shows  similar 
effects  of  pitch-flap  coupling  on  the  stability  of  the  rotor  blade  and  hub 
modes.  The  only  difference  in  modal  damping  at  an  advance  ratio  of  0.50 
is  shown  by  the  blade  edgewise  mode  for  tail  No.  1.  The  stability  of  this 
mode  remains  unchanged  throughout  the  range  of  coupling  examined. 

The  modal  frequency  response  with  blade  pitch-flap  coupling  is 
illustrated  in  Figures  1*3  and  1*1*.  From  these  figures,  it  can  be  seen 
that  the  blade  flatwise  frequency  for  the  symmetric  and  unsymmetric  modes 
Increases  with  pitch- flap  coupling  for  both  tail  rotors.  On  the  other 
hand,  the  presence  of  pitch-flap  coupling,  positive  or  negative,  results 
in  lower  blade  torsional  frequencies,  as  indicated  in  Figures  1*3  and  1*1*. 

The  frequencies  of  the  blade  edgewise  and  rotor  hub  modes  are  little 
affected  by  pitch-flap  coupling  variations. 
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Figure  39.  Effect  of  Pitch-Flap  Coupling  on  Modal  Damping 
for  Tail  Rotor  No.  1 in  Hover. 
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PITCH- PUP  COUPLING,  I,  .DCS 


(•)  Blade  SjOMtrlc  Nodes 


PITCH  -n.  *P  COUPU  HO , I , , DCO 


(b)  Blade  Whirl  and  Hub  Nodus 


Figure  1|0 • Effect  of  Pitch-Flap  Coupling  on  Modal  Damping 
for  Tail  Rotor  No.  2 in  Hover. 
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Effect  of  Pitch-Flap  Coupling  on  Modal  Damping 
for  Tail  Rotor  No.  1 at  p * .50. 


TORWARO  WHIM. 


Cb)  Blade  Whirl  and  Hub  Modes 
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Figure  UU . Effect  of  Pitch— Flap  Coupling  on  Mod&l  Freouency 
for  Tail  Rotor  No.  2 in  Hover. 
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Effect  of  Blade  Edgarlee  Natural  Frequency 


The  blade  edgewise  natural  frequency  was  varied  by  changing  the  radial 
distribution  of  the  edgevlse  area  moment  of  inertia.  These  variations 
gave  a frequency  range  of  0.66  to  2.63  cycles/rev  for  tail  rotor  Ko.  1 and 
0.68  to  2.66  cycles/rev  for  tail  rotor  No.  2.  Four  values  were  used. 

The  nodal  damping  and  frequencies  obtained  for  the  tvo  rotors  in 
hover  and  at  advance  ratio  of  0.50  are  presented  in  Figures  45  through  50. 

From  Figures  45(a)  and  46(a),  it  is  observed  that  the  hover  stability 
of  the  blade  synmetric  flatvise  and  torsional  modes  is  little  affected  by 
changes  in  the  edgevlse  natural  frequency.  The  blade  symmetric  edgevlse 
mode  exhibits  a slight  reduction  in  modal  damping  for  edgewise  frequencies 
below  1.8  eycles/rev  for  both  rotors.  However,  as  the  edgewise  frequency 
is  Increased  further,  this  mode  becomes  considerably  more  stable. 

The  frequency  response  of  the  blade  symmetric  modes  for  variations  in 
edgewise  natural  frequency  appears  in  Figures  49  and  50.  The  frequencies 
of  the  torsional  and  flatvise  modes  of  both  rotors  remain  constant  through- 
out the  range  of  edgewise  frequency  investigated.  The  frequency  of  the 
symmetric  edgewise  mode  is  almost  equal  to  the  edgewise  blade  natural 
frequency  for  both  rotor  systems . A coalescence  of  the  blade  symmetric 
flatvise  and  edgevlse  modal  frequencies  occurs  at  an  edgevlse  natural 
frequency  of  1.3  cycles /rev. 

The  hover  stability  characteristics  of  the  blade  whirl  and  rotor  hub 
modes  are  illustrated  in  Figures  45(b)  and  46(b)  for  tail  rotors  No.  1 
and  2 respectively.  From  these  figures,  it  is  observed  that  the  trend  of 
modal  damning  with  blade  edgevlse  natural  frequency  is  quite  different 
for  the  tvo  baseline  rotors.  Tail  rotor  No.  1 exhibits  large  changes  in 
system  stability  as  the  edgewise  natural  frequency  is  increased  above  1.6 
cycles/rev,  while  the  stability  behavior  of  tail  rotor  No.  2 shows  signif- 
icant variations  for  edgewise  natural  frequencies  below  1.6  cycles/rev. 

The  backward  whirl  mode  of  the  tail  rotor  No.  1 is  marginally  stable  in 
the  range  of  w^/fl  from  0.66  to  1.80  cycles /rev.  Both  the  forward  and 

backward  whirl  edgewise  modes  become  quite  staole  at  higher  edgewise 
natural  frequencies,  while  the  torsional  response  shows  a large  decrease 
in  stability.  This  behavior  may  be  caused  by  the  interaction  between  the 
torsional  and  edgewise  inodes  shown  in  Figure  49(b).  The  blade  flatwise 
whirl  mode  and  the  rotor  hub  pitch  and  yaw  modes  for  tail  rotor  No.  1 are 
not  significantly  affected  by  variations  in  blade  edgewise  natural 
frequency.  The  stability  characteristics  presented  in  Figure  46(b)  for 
tail  rotor  No.  2 indicate  that  a strong  coupling  is  present  betveer.  the 
blade  flatwise  and  edgewise  forward  whirl  modes  and  the  rotor  hub  pitch 
mode  for  edgewise  natural  frequencies  below  1.6  cycles/rev.  The  stability 
of  the  blade  flatwise  forward  whirl  and  the  hub  pitch  inodes  decreases, 
while  that  of  the  blade  edgewise  forward  whirl  mode  increases  in  this 
frequency  range.  From  Figure  50(b),  it  is  noted  that  the  frequencies  of 
these  modes  are  all  near  each  other.  The  response  of  the  blade  torsional 
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■ode  1*  not  affected  by  chances  In  the  blade  edgewise  natural  frequency. 

This  result  is  not  surprising,  since  the  blade  torsional  whirl  frequency 
is  well  removed  from  all  other  system  frequencies,  as  shown  in  Figure  50(b). 

The  stability  trends  with  blade  edgewise  natural  frequency  presented 
in  Figures  1*7  and  48  at  an  advance  ratio  of  0.50  show  behavior  similar  to 
that  discussed  for  the  hover  condition.  The  only  exception  is  the  blade 
edgewise  forward  whirl  mode  for  tail  rotor  So.  1,  which  is  unstable  for 
ttES/0  leM  than  one  cycle/rev.  It  nay  be  added  that  the  variations  in 

stability  with  blade  edgewise  natural  frequency  are  not  as  pronounced  for 
the  forward  flight  condition  as  for  hover. 
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(•)  Had*  Sywttrlc  Nod** 


(b)  Had*  Whirl  and  Hub  Nodas 


Figure  k6.  Effect  of  Blade  Edgewise  Natural  Frequency  on 
Modal  Damping  for  Tail  Rotor  No.  2 in  Hover. 
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(«)  Blade  Syattrk  Mode* 


(b)  Blade  Whirl  and  Hub  Modes 


Figure  **7.  Effect  of  Blade  Edgewise  Natural  Frequency  on 
Modal  Damping  for  Tail  Rotor  No.  1 at  y = .50. 
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(a)  Slide  Symetrlc  Nodes 


(b)  lied*  Whirl  end  Hub  Nodes 


9.  Effect  of  Blade  Edgewise  Natural  Frequency  on 
Modal  Frequency  for  Tail  Rotor  No.  1 in  Hover. 
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(•)  Blade  Synatrlc  Aode* 

8 


(6)  Hid*  Whirl  and  Hub  Node* 


Figure  50.  Effect  of  Blade  Edgewise  Natural  Frequency  on 
Modal  Frequency  for  mail  Rotor  No.  2 in  Hover. 
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Effect  of  Rotor  Hub  Frequencies 


The  effect  of  variations  in  reference  rotor  hub  frequency  on  the 
stability  characteristics  of  the  tvo  baseline  tail  rotors  is  illustrated 
in  figure  $1  through  56  in  hover  and  at  an  advance  ratio  of  0.50.  The 
generalised  stiffness  of  the  fixed  system  yav  and  pitch  modes  vaa  varied 
to  give  a range  in  reference  yav  frequency  of  0.625  to  20.0  cycles/rev 
for  rotor  Ho.  1 and  0.52  to  16.6U  cycles/rev  for  rotor  Ho.  2.  Pour  values 
of  stiffness  were  used.  In  Figures  51  through  56*  the  reference  rotor  hub 
pitch  frequency  is  also  shown;  it  ranges  from  0.18  to  5.76  cycles  /rev  for 
rotor  Ho.  1 and  from  0.W»5  to  lk.2h  cycles /rev  for  rotor  No.  2.  The  modal 
damping  and  frequency  for  the  blade  vhirl  and  rotor  hub  modes  only  are 
shovn  in  these  figures  since*  for  the  modes  used  here,  the  fixed  system 
motions  do  not  affect  the  blade  symnetric  modes . 

The  stability  characteristics  in  hover  are  presented  in  Figures  51 
and  52  for  tall  rotors  Ho.  1 and  No.  2 respectively.  From  these  figures, 
it  can  be  seen  that  there  is  no  appreciable  change  in  the  modal  damping 
of  the  blade  vhirl  modes  at  reference  rotor  hub  frequencies  above  U cycles/ 
rev.  Belcv  this  frequency,  hovever,  the  interaction  between  the  fixed 
system  motions  and  the  blade  flatwise  and  edgewise  modes  becomes  more 
pronounced*  especially  for  tail  rotor  No.  2.  This  effect  may  be  expected 
from  an  inspection  of  the  modal  frequency  characteristics  presented  in 
Figures  55  and  56.  At  low  reference  hub  frequencies,  the  blade  flatwise 
backward  whirl  mode  for  tail  rotor  No.  2 has  a sharp  reduction  in  modal 
damping , while  the  blade  edgewise  backward  vhirl  becomes  more  stable.  For 
both  rotors,  the  stability  of  the  rotor  hub  pitch  and  yav  modes  increases 
significantly  at  low  reference  hub  frequencies. 

The  stability  characteristics  at  an  advance  ratio  of  0.50  are  shovn  in 
Figure  53  for  tail  rotor  No.  1 and  in  Figure  51*  for  tall  rotor  No.  2.  A 
comparison  of  these  forward  flight  results  with  the  hover  results  shows 
no  major  differences  in  the  parametric  effect  of  reference  hub  frequency 
on  modal  damping. 

The  frequency  characteristics  of  both  baseline  rotors  as  a function 
of  reference  hub  pitch  and  yav  frequencies  appear  in  Figures  55  and  56. 

It  can  be  seen  that  reference  hub  frequency  has  generally  a negligible 
effect  on  the  blade  vhirl  frequencies.  Hovever,  as  previously  discussed, 
an  interaction  between  the  hub  and  blade  modes  is  present  for  tail  rotor 
No.  2 at  reference  hub  frequencies  below  h cycles/rev.  Linear  variations 
in  modal  frequency  for  the  rotor  hub  pitch  and  yav  modes  with  reference 
hub  frequency  are  shown  in  these  figures  for  both  baseline  tail  rotors. 

The  calculated  frequencies  are  slightly  lover  than  the  reference  values 
due  to  the  mass  and  inertia  effects  of  the  rotors. 
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Tail  Rotor  No.  1 in 
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Figure  52 . Effect  of  Rotor  Hub  Frequency  on  Modal  Damping 
for  Tail  Rotor  No.  2 in  Horer. 
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Figure  53.  Effect  of  Rotor  Hub  Frequency  on  Modal  Damping 
for  Tail  Rotor  Ho.  1 at  y ■ .50. 


Figure  54 . Effect  of  Rotor  Hut  Frequency  on  Modal  Damping 
for  Tail  Rotor  No.  2 at  y = .50. 


for  Tail  Rotor  No.  1 in 
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Figure  56.  Effect  of  Rotor  Hub  Frequency  on  Modal  Frequency 
for  Tail  Rotor  Ho.  2 in  Hover. 


Effect  of  Blade  Precone  Angle 


The  effect  of  blade  precone  angle  on  the  atability  characteristics 
of  tall  rotor  Ho.  2 Is  presented  in  Figures  5T  through  39.  Only  tall  rotor 
No.  2 vas  considered  for  this  portion  of  the  study  since  it  is  nonarticu- 
lated  In  the  flatwise  direction  and  can  thus  be  preconed.  A range  in  blade 
precone  from  -6  to  6 degrees  vas  examined.  Five  values  vere  used.  The 
results  from  Figure  57(a)  in  hover  show  that  blade  precone  has  little  effect 
on  the  modal  damping  of  the  blade  symmetric  flatwise  and  edgewise  modes , 
while  the  torsional  mode  exhibits  an  increase  in  stability  as  the  blade  is 
preconed.  However,  from  Figure  38(a)  at  an  advance  ratio  of  0.30,  the 
blade  symmetric  flatwise  mode  indicates  an  increase  in  modal  damping  vith 
blade  precone,  while  the  torsional  mode  becomes  more  stable  than  the  hover 
case. 

The  variations  in  modal  damping  with  precone  for  the  blade  whirl  and 
rotor  hub  modes  appear  in  Figures  37(b)  and  38(b)  for  advance  ratios  of 
zero  and  0.30  respectively.  From  these  figures,  it  is  observed  that  the 
blade  torsional  whirl  mode  is  least  stable  for  zero  precone,  a behavior 
similar  to  that  of  the  symmetric  torsional  mode.  An  increase  in  precone 
reduces  the  stability  of  the  blade  edgewise  vhlrl  mode,  while  the  flat- 
wise vhirl  mode  becomes  more  stable;  these  effects  are  generally  greater 
for  the  forward  flight  condition.  The  rotor  hub  modes  show  generally  small 
variations  in  mode  damping  response  with  blade  precone  except  for  the 
pitch  mode  in  hover.  The  dasping  of  the  pitch  mode  decreases  from  6 to  U 
percent  as  precone  is  increased  from  -6  to  6 degrees. 

A negligible  effect  of  blade  precone  on  modal  frequency  is  exhibited 
by  the  results  shown  in  Figure  59(a)  for  the  blade  symmetric  modes  and  in 
Figure  39(b)  for  the  blade  unsymmetric  and  rotor  hub  modes. 


95 


MOOAL  DAMPING.  PERCENT 


BLADE  PRECONE  , 0 .DEG 
0 


(a)  Blade  Symtrlc  Nodes 


BLADE  PRECONE  , , DEG 


(b)  Blade  Whirl  and  Hub  Nodes 


Figure  58 . Effect  of  Blade  Precone  on  Modal  Damping 
for  Tail  Rotor  No.  2 at  p = .50. 
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(•)  Blade  Synetrlc  Nodes 


(b)  Blade  Whirl  and  Hub  Nodes 


Figure  59 • Effect  of  Blade  Precone  on  Modal  Frequency 
for  Tail  Rotor  No.  2 in  Hover. 
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Effect  of  Pitch-Lag  Coupling 


Variation?  in  pitch- lag  coupling  over  the  range  -0.2  to  0.2  for  both 
baseline  rotors  in  hover  and  at  an  advance  ratio  of  0.5  produced  no  dis- 
cernible change  in  stability.  The  fact  that  both  rotors  are  relatively 
s+.iff  inplane  is  the  cause  of  this  effect.  However,  it  is  known  that 
pitch-lag  coupling  can  appreciably  affect  the  stability  of  soft  inplane 
rotors.  Therefore,  in  such  cases  the  parameter  should  be  accurately 
defined  and  incorporated  in  any  analysis. 
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COHCLUDIHO  REMARKS 


The  major  objectives  of  this  study  were,  first,  to  develop  a computer- 
ized eigenvalue  analysis  that  had  a sufficiently  complete  dynamic  and 
aerodynamic  description  of  coupled  rotor/airframe/control  systems  to  allov 
accurate  stability  estimations  and,  second,  to  demonstrate  some  of  its 
capabilities  through  parametric  and  sensitivity  studies  of  representative 
semiarticulated  and  hingeless  tail  rotor  systems.  Both  of  these  objectives 
were  achieved. 

The  studies  described  in  the  foregoing  were  necessarily  restricted  to 
simplified  systems,  tut  they  do  serve  to  illustrate  the  potential  of  the 
analysis.  The  inclusion  of  forward  flight  effects  in  the  analysis  is  a 
big  step  forward  and  considerably  extends  the  envelope  within  which  the 
analysis  can  be  applied.  Although  not  exercised  in  these  studies,  the 
unsteady  aerodynamics  serve  to  further  refine  the  modelling  of  aerodynamic 
effects,  thus  permitting  flutter  and  vake  interaction  studies.  The  models 
used  for  the  control  systems  can  accurately  represent  the  dynamics  of  most 
existing  systems,  thereby  reducing  the  number  of  effects  not  properly 
accounted  for  in  the  overall  system.  The  overall  usefulness  of  the  analysis 
was  further  extended  by  inclusion  of  the  gimbaled  rotor  capability. 

Since  a number  of  unstable  phenomena  involve  coupling  between  rotor 
system  and  fixed  system  motions,  it  is  important  that  the  analysis  have 
the  facility  to  predict  these.  This  is  adequately  accounted  for  by  inclu- 
sion of  an  airframe  (or  fixed  system)  model  that  can  represent  rigid  tody 
and  elastic  modes,  including  local  transmission,  etc.,  modes. 

The  methods  employed  in  the  development  of  the  analysis  are  given  in 
the  appendixes,  together  with  the  final  equations  of  motion. 

A final  comment  on  the  computer  program  of  the  analysis:  although 

this  in  its  entirety  gives  a refined  description  of  coupled  rotor /air frame/ 
control  systems,  it  was  realized  that  there  are  many  instances  when  studies 
esploying  simple  systems  are  desired.  To  this  end,  the  program  is  so 
designed  that  the  user,  by  use  of  simple  control  cards,  can  specify  the 
precise  degree  of  freedom,  type  of  control  system,  mode  of  flight,  and 
type  of  aerodynamics  that  he  wishes  to  use.  From  this  point,  the  program 
assumes  executive  control  and  automatically  performs  the  desired  functions. 
This  feature  allows  systems  with  as  little  as  one  degree  of  freedom  to  be 
analyzed  and  makes  it  possible  for  inexperienced  users  to  run  the  program 
with  little  difficulty. 
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RECOMMENDATIONS 


The  analysis  described  In  this  report  represents  a fairly  complete 
description  of  a coupled  rotor/airframe/control  system.  Its  capacity  to 
treat  rotor  systems  In  hover  and  forward  flight  has  been  demonstrated. 

However,  the  effects  of  unsteady  aerodynamics  and  cyclic  inputs  to  the 
dynamic  system  equations  have  not  been  explored.  In  addition,  the  represen- 
tation of  the  blade  pitch  control  system  In  the  parametric  and  sensitivity 
studies  did  not  account  for  multi -degree -of -freedom  effects  or  the  effects 
of  the  control  servos  and  servo  feedback  parameters.  The  representation 
of  hub  motions  was  also  idealized. 

Clearly,  In  practice  such  exclusions  and  simplifications  may  be  neither 
warranted  nor  desired.  If  ve  wish  to  define  the  overall  system  stability 
characteristics  to  the  best  of  our  ability,  and  within  the  constraints  of 
the  accuracy  of  the  analysis  ve  are  employing,  it  behooves  us  to  model  that 
system  as  accurately  as  possible.  Therefore,  ve  must  include  in  our  studies 
all  of  the  elements  that  ve  feel  could  play  an  important  part  in  our  predic- 
tions . 

It  is,  therefore , recommended  that  the  effects  of  inclusion  of  unsteady 
aerodynamics  and  cyclic  inputs  to  the  dynamic  system  equations  be  established. 
Recent  studies  using  a simplified  Floquet  analysis  have,  for  example,  shown 
that  the  inclusion  of  cyclic  pitch  had  a pronounced  effect  on  the  pitch  lag 
stability  of  a rotor.  In  this  case,  the  classical  hover  expressions  were 
nonconservative . 

It  is  also  recommended  that  the  analysis  be  subjected  to  extensive  use. 

It  is  only  by  vising  the  analysis  and  possibly  even  questioning  certain 
aspects  of  the  modelling  that  the  user  can  become  familiar  enough  with  its 
general  characteristics  to  effectively  incorporate  improvements.  Also,  it 
is  only  through  use  that  a complete  understanding  of  the  real  meaning  of 
the  output  can  be  gained. 

In  its  present  form,  fhe  computer  program  of  the  analysis  relies  on 
some  unlinked  transient  aero elastic  analysis  for  production  of  the  deriva- 
tive terms  required  for  forward  flight  studies.  The  analysis  currently 
used  for  this  purpose  is  the  Y-200  Normal  Modes  Analysis  provided  to  the 
Amy  under  Contract  DAAJ02-71-C-002U . Unfortunately,  this  analysis  has 
some  inconsistencies  with  the  E-927  Aeroelastic  Rotor  Stability  Analysis . 

Not  the  least  of  these  are  the  facts  that  Y-200  can  accept  only  one  blade 
angle  distribution  and  that  it  does  not  permit  representation  of  control 
system  dynamics.  It  is,  therefore,  recommended  that  these  and  any  other 
inconsistencies  be  eliminated  and  that  Y-200  be  linked  to  E-927  to  facilitate 
use  of  the  program. 
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Finally,  any  analysis  is  only  as  good  as  the  confidence  ve  place  in 
its  predictions.  Obviously,  such  confidence  can  only  be  obtained,  or 
Increased,  by  correlating  the  analysis  vith  test  data.  Unfortunately,  only 
Halted  correlation  of  the  present  analysis  has  been  completed.  Much 
reaains  to  be  done.  To  close  this  gap,  it  is  suggested  that  consideration 
be  given  to  the  initiation  of  a model  tast  program  to  provide  the  needed 
data.  Such  a program  could  be  designed  tc  accurately  model  not  only  the 
rotor  blades  but  also  the  control  system  and  the  test  rig  characteristics. 
By  making  parametric  variations  of  the  system  dynamic  characteristics 
and  conducting  hovering  and  forvard  flight  tests,  very  valuable  data  could 
be  acquired.  Since  the  analysis  is  capsble  of  accurately  modelling  the 
systems,  very  meaningful  correlation  could  be  obtained. 
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APPENDIX  A 


HOVER  ANALYSIS 


This  appendix  gives  details  of  the  development  of  an  analysis  designed 
to  study  the  stability  of  helicopter  sain  or  tall  rotors  in  hover,  or  under 
conditions  of  pure  axial  flow.  It  can  also  be  used  to  study  the  stability 
of  propellers. 

Important  assumptions  are  noted  belov  and  the  coordinate  system 
employed  is  shown  in  Figure  Al.  The  coordinates  used  for  the  control 
systems  are  defined  later  in  the  text. 

ASSUMPTIONS 

1.  Dynamic  and  aerodynamic  effects  assume  small  perturbations  about  steady 
initial  values  of  the  system  generalized  coordinates. 

2.  Aerodynamic  forces  are  developed  using  strip  theory. 

3.  Radial  flow  effects  are  neglected. 

U.  Products  involving  up  to  the  squares  of  steady  displacements  and  a 
perturbation  displacement,  velocity,  or  acceleration  are  retained. 

5.  Rotor  speed  is  constant. 

6.  No  small-angle  assumptions  are  applied  to  the  blade  collective  pitch 
or  twist. 

7.  Because  of  the  rotating  system  generalized  coordina.^  transformations 
employed,  only  rotors  with  three  or  more  blades  can  be  treated. 

8.  Stall,  compressibility , and  unsteady  effects  are  included. 

9.  Inflow  is  constant  over  the  rotor  disc. 

10 . Blades  have  an  elastic  axis . 

11.  Blade  flap  and  lag  hinges  are  coincident. 

12.  The  blade  feathering  bearing  can  either  remain  in  the  plane  of  the  hub 
or  follow  the  blade  root  out-of-plane  slope.  A similar  statement  can 
be  made  regarding  inplane  motion. 

13.  The  rotor  shaft  experiences  no  torsional  deformations. 


» r 4w*->4* 
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lU.  Only  those  aerodynamic  forces  associated  vith  the  rotor  system  are 

considered.  Rotor/airframe  Interference  effects  and  airframe  lift/drag 
are  not  Included. 

DEVELOPMENT  OF  HOVER  ANALYSIS 


Ejynamic  System  Equations 

The  development  of  the  equations  of  motion  of  the  total  dynamic  system 
follows  the  classical  path  of  first  defining  the  kinetic  and  potential 
energies  and  dissipation  potentiala  of  the  rystem,  then  establishing 
Lagrangefe  equations  in  the  form 


d,3T  > 3T  3V  + 3D 

' aij  34j  ’ 


CAl) 


In  what  follows,  to  preserve  clarity,  the  total  dynamic  system  ib 
broken  down  into  four  major  subsystems.  These  are: 

1 . Hie  rotor  blades . 

2.  The  fixed,  or  airframe,  system. 

3.  The  control  systems. 

h.  The  servo  systems. 

Each  of  these  Is  treated  separately,  but  it  will  be  apparant  that 
their  dynamic  responses  are  all  close  coupled  through  the  motions  of  the 
respective  coupling  generalized  coordinates. 

Equation  A1  is  indeed  a simple  statement  of  fact,  but  the  development 
of  the  terms  on  the  left-hand  side  does,  for  a multi-degree-of- freedom 
system,  present  a very  formidable  task  that  requires  extensive  expansion 
and  manipulation  of  large  arrays  of  matrices.  In  the  interests  of  documen- 
ting the  analysis  in  a not  too  cumbersome  fashion,  obvious  intermediate 
steps  are  omitted,  but  sufficient  detail  is  retained  to  3how  how  the  final 
forms  of  the  equations  of  motion  were  established. 

In  the  development  of  the  equations , a number  of  transformation 
matrices  and  vectors  are  employed.  These  are  defined  below  for  easy  refer- 
ence. 
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Relations  A2 

through  A9  are  the 

coordinate  transformation 

matrixes . 

A10  is  the  matrix  of  the  blade  flatwise,  torsional,  and  edgewise  mass 
moments  of  inertia  about  the  local  center  of  gravity.  All  defines  the 
uncoupled  blade  bending  coordinates  including  the  radial  shortening  caused 
by  bending.  A12  is  a velocity  vector  that  includes  the  effectB  of  bending 
and  pitching  velocities,  steady  bending  displacements  and  pitch  angles, 
and  center  of  gravity  offset.  A13  defines  the  local  radial  velocity. 
Relations  AlU  through  A27  are  self-explanatory  definitions  of  various 
coordinates  and  vectors. 

The  coordinates  were  transformed  as  follows : 

(1)  Translate  through  CG. 

(2)  Rotate  through  uncoupled  blade  flatwise  and  edgewise  slopes  and  inte- 
grate to  define  the  vector  All. 

(3)  Rotate  through  blade  pitch  and  define  the  coupled  coordinates  given 
by  AlU  and  A15 . 

(U)  Rotate  through  blade  rigid-body  lag. 
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(5)  Rotate  through  blade  rigid-body  flap. 

(6)  Translate  through  offset. 

(7)  Rotate  through  azimuth  tingle. 

(8)  Rotate  through  hub  pitch. 

(9)  Rotate  through  hub  roll. 

(10)  Translate  through  hub  vertical,  lateral,  and  longitudinal  displacements. 
Derivation  of  Subsystem  Dynamic  Characteristics 

Rotor  Blades 


Consider  first  the  kinetic  energy.  This  can  be  broken  down  into  two 
parts:  that  associated  with  translational  motion  and  that  associated  with 

pure  rotational  motion. 

Trans lantional  Kinetic  Energy 

Using  the  coordinate  system  shown  in  Figure  A1  and  performing  the 
required  coordinate  rotations  and  translations,  it  can  be  shown  that  the 
absolute  translational  motion  of  a blade  element  i at  blade  station  rn  j 
on  any  blade  n is  given  by 


A^YAqexA*,n(A0,nAY,nqn,i  + V + Q 


(A28) 


where 


Ae  .n.i^n.i 

{ " *E,n,i’  rn,i * qF,n,i  } 


(A28a) 


Taking  the  first  time  derivative  of  ^ ^ the  total  kinetic  energy 
arising  from  translational  motion  of  all  elements  of  all  blades  can  be 
written 
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Rotational  Kinetic  Energy 

Since  angular  velocity  follows  the  rule  $a  ■ B + where  $a  b 
denotes  the  angular  velocity  of  the  coordinate  system  A relative  to  that 
of  the  coordinate  aystem  B,  the  angular  velocity  components  ol  the  blade 
can  be  expressed  as  the  sim  of  the  relative  angular  velocities  between  each 
of  the  coordinate  systems  given  in  Figure  Al.  Thus,  the  absolute  angular 
velocity  of  any  blade  element  i at  blade  station  r . on  any  blade  n is 
given  by  n‘* 


NE 


NE 
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The  various  components  of  this  velocity  vector  are  readily  developed 
if  we  recognize  that  the  ik,  Jk,  k^  are  the  unit  vectors  associated  with 
their  corresponding  rotational  transformation  matrices.  Therefore,  to 
obtain  we  must  express  the  ik,  Jjj.,  kk  in  terms  of  i2»  J2»  k2  . Using 

the  transformation  matrices  defined  earlier,  we  can  write 
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Substitution  of  the  above  in  Equation  A30  expresses  the  absolute 
rotational  velocity  in  terms  of  the  generalized  coordinates  and  their 
derivatives . 

. These  relationships  allow  us  to  write  the  absolute  velocity  vector 

. in  terms  of  the  appropriate  unit  vectors  at  a blade  element  1 . 
n,i 

Thus,  the  total  kinetic  energy  arising  from  pure  rotational  motion 
of  all  elements  of  all  blades  can  be  written 


Total  Kinetic  Energy 

The  total  kinetic  energy  arising  from  translational  and  pure  rotational 
motion  of  all  elements  of  all  blades  is  given  by  the  sum  of  Equations  A 29 
and  A32,  namely. 


BTOT 


BT 


BR 


(A33) 


Potential  Energy 

The  blades  have  two  direct  sources  of  potential  energy:  that  associ- 

ated with  blade  bending  and  that  associated  with  twisting  deformations. 

The  bending  potential  energy  of  a blade  can  be  written  in  the  classical 
form  which  expresses  this  energy  'n  terms  of  the  fourth  spanvise  derivative 
of  the  flatwise  and  edgewise  deformations.  However,  since  in  this  analysis 
a modal  approach  is  employed  and  it  is  assumed  that  all  modal  informs'  ion 
relating  to  the  blades  is  available  for  inclusion  in  the  elgenanalysi  .,  the 
total  potential  energy  associated  with  bending  of  all  blades  in  simply 
written  as 
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The  potent ial  energy-  derived  from  torsional  deformations  of  all  blades 
is  given  by 
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In  this  expression  Kn,i  is  the  torsional  stiffness  of  the  blade  n at 
the  blade  station  rn>i  . It  is  obtained  from  the  blade  polar  second  moment 
of  area  distribution  and  the  blade  torsional  modulus  of  elasticity  using 
the  relationship 


(A36) 


Further  sources  of  potential  energy  are  contained  in  the  root  springs 
that  may  be  used  when  the  blades  are  considered  inflexible.  These  root 
springs  restrain  both  flapping  and  inplane  motion  and  appear  explicitly 
only  in  the  rigid-body  flap  and  lag  equations . The  total  potential  energy 
associated  with  the  springs  is  given  by 


N 


It  should  be  noted  that,  when  the  blade  elastic  modes  are  used,  the 
potential  energy-  from  these  springs  is  implied  in  Equation  A3U  since  the 
springs  determine  the  blade  root  boundary  conditions . 

Total  Potential  Energy 

The  total  potential  energy  arising  from  bending,  twisting,  and  rigid- 
body  motions  of  all  blades  is  given  by  the  sum  of  Equations  A3b,  A35,  and 
A37.  namely. 
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Dissipation  Potentials 

We  can  immediately  recognize  tvo  sources  of  energy  dissipation  associ- 
ated with  the  blades:  that  resulting  from  bending  deformations  and  that 

resulting  from  torsional  deformations.  Unfortunately,  no  realistic  analy- 
tical definition  of  structional  damping  of  any  kind  has  ever  been  established. 
Nevertheless,  it  is  known  that  structural  damping  does  exist,  no  matter  how 
small.  Therefore,  this  msy  be  included  as  a percentage  of  critical  damping. 

We  must,  however,  remenfcer  that  in  specifying  a percentage  of  critical 
damping,  we  must  also  specify  the  frequency  upon  which  it  is  based.  In  the 
case  of  blade  bending,  this  presents  no  problem  since  the  modal  frequencies 
are  obvious  choices.  However,  in  the  case  of  blade  torsion,  the  choice  is 
not  so  obvious  since  coupling  with  the  control  system  can  cauBe  radical 
changes  in  the  torsional  natural  frequency . To  circumvent  this  problem, 
the  torsional  damping  is  based  on  the  rotor  rotational  speed,  which  iB  always 
known.  It  is  then  the  privilege  of  the  user  of  the  analysis  to  choose  the 
damping  level  that  he  feels  most  closely  satisfies  the  condition  being 
analyzed. 

Based  on  the  above,  we  can  then  write  the  total  dissipation  potential 
arising  from  blade  bending  as 
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and  that  arising  from  blade  torsion  as 
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Energy  is  also  dissipated  by  the  lag  dampers  that  act  on  the  rigid- 
body  inplane  degree  of  freedom.  This  dissipation  potential  is  given  by 
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Total  Dissipation  Potential 

The  total  energy  dissipation  arising  from  bending,  twisting,  and  rigid- 
body  inplane  motion  of  all  blades  is  given  by  the  sum  of  Equations  A39,  AbO, 
and  Abl,  namely. 
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Fired,  or  Airframe,  System 

Neglecting  aerodynamic  interference  effects  betveen  che  rotor  system 
and  any  part  of  the  airframe  and  also,  for  the  present,  neglecting  mechanical 
control  coupling,  the  only  coupling  that  can  exist  betvr  eu  rotor  motions 
and  airframe  motions  is  that  transmitted  through  the  rotor  hub.  There- 
fore, knowing  the  dynamic  characteristics  of  the  airframe  as  seen  at  the 
rotor  hub  allows  us  to  define  a dynamically  coupled  rotor/airframe  system. 

Such  an  approach  is  used  in  this  analysis. 

It  is  assumed  that  the  airframe,  or  fixed  system,  dynamic  characteristics 
are  available  either  from  shake  test  data  or  from  a separate  analysis.  The 
data  required  for  each  of  the  NA  airframe  modes  <f  are: 

1.  The  generalized  mass,  MA> 

2.  Hie  modal  frequency,  <*>A. 

3.  The  percentage  of  critical  damping,  cA> 

k.  The  modal  components  of  motion  at  the  rotor  hub,  $x,  $z , 

♦ex*  *6Y’ 

With  this  Information,  we  are  able  to  define  the  total  kinetic  and 
potential  energies  and  the  dissipation  potential  associated  with  the  air- 
frame system  as  follows: 

Total  Kinetic  Energy 
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Total  Potential  Energy 
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The  modal  components  at  the  rotor  hub  provide  the  motion  Interface 
between  the  rotor  and  airframe  systems . 

Control  Systems 

This  analysis  cor  aiders  both  main  and  tail  rotor  systems.  Normally , 
the  control  systems  associated  with  these  have  little  in  common.  For 
example,  the  conventional  tail  rotor  has  only  collective  pitch  inputs.  An 
assortment  of  cables , quadrants , control  rods , and  bell  cranks  may  also 
feature  in  the  design.  The  control  servos  can  be  close  to,  or  some 
considerable  distance  away  from,  the  rotor.  On  the  other  hand,  the  conven- 
tional main  rotor  has  both  collective  and  cyclic  pitch  inputs,  necessitating 
the  use  of  a svash  plate  arrangement.  The  primary  control  servos  normally 
act  directly  on  the  stationary  svash  plate. 

To  accommodate  the  peculiar  characteristics  of  both  types  of  control 
system,  each  is  modelled  and  analyzed  separately. 

Hie  equations  for  the  main  and  tail  rotor  servos  are  identical. 
Therefore,  they  are  excluded  from  this  section  and  developed  later. 

Tail  Rotor  Control  System 

Figure  A 2 shows  the  model  used  to  describe  this  type  of  control 
system.  The  model  simulates  the  collective  pitch  spider  beam  and  actuator 
shaft  backed  up  by  three  spring,  mass,  damper  systems  that  can  be  used  to 
describe  the  dynamics  of  control  rod/quadrant /cable  arrangements.  The 
subsidiary  spring,  mass,  damper  systems  can  he  used  to  describe  the 
dynamics  of  appendages  such  as  control  surfaces  and  rudders. 

In  Figure  A2,  X is  the  motion  at  the  pushrod  of  blade  n.  It  is 
defined  by  1 ’ 
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Figure  A 2.  Tail  Rotor  Control  System  Model. 
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It  is  assumed  that  the  actuator  and  collective  pitch  spider  beam  move 
with  an  amplitude  X2  + qz  in  the  Z direction;  the  effects  of  all  other 
components  of  hub  notion  are  ignored.  It  is  also  assumed  that  airframe 
motions  can  be  neglected  in  the  remainder  of  the  control  system. 

The  energies  associated  vith  the  collective  motion  of  the  tall  rotor 
control  system  are  then  readily  written  as  follows: 
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Total  Dissipation  Potential 


TTOT 


■ HA  . , 

HA 


°2(i2  * ' *»)2  * °’(i3  ‘ i*’ 


A \2 


♦ CA  - *S>'  * Cs(i5  - V2 


1 i2 


. • • .7  < • • . 

+ c (X  - X ) ♦ C (X  - X ) 

31v  31  3 4l'  41  4 ' 


♦ C51(t51  - i5>2) 


(AU9) 


For  pure  cyclic  notion  of  the  rotating  system  generalized  coordinates, 

it  is  assumed  that  the  actuator  shaft  is  subjected  to  a pure  moment  at  its 

point  of  connection  to  the  collective  piten  spider  beam.  Under  these 

circumstances , the  control  system  stiffness  to  be  associated  with  the 

motion  X,  is 
1 »n 


c,n 


* 2Kma> 


(A50) 


Therefore,  the  total  potential  energy  associated  with  cyclic  motion 
of  the  control  system  is  given  by 


VTCT3T  * ^i^c.n^.n^ 


(A51) 


Main  Rotor  Control  System 

The  model  used  to  describe  the  main  rotor  control  system  is  shown  in 
Figure  A3.  This  is  depicted  as  a rigid  svash  plate  supported  on  three 
spring/servo  systems . The  servos  are  situated  at  a radius  Rs  from  the 
center  of  the  swash  plate  with  the  forward  servo  being  set  at  an  angle 
<5pg  relative  to  the  Y axis  as  shown.  The  blade  pushrods  are  situated 
at  a radius  Rg  from  the  center  of  the  swash  plate.  The  point  p locates 
the  position  of  the  pushrod  of  blade  n,  and  the  angle  that  positions 

this  point  relative  to  the  Y axis  is  given  by 


B,n 


SB,n  " «re  ♦ » ♦ *(n  - l>/»  ♦ «pRfn 


(A52) 
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Figure  A3.  Main  Rotor  Control  System  Model. 
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vh  era 


- tan_1[L,  /( e + L,  tan6  )] 

* n |H  ^ |Q  Q ^ |Q  3 1 n 


(A53) 


This  definition  neglects  steady  and  dynamic  deformations . The  servos 
at  XL  and  X*  are  positioned  90  degrees  and  180  degrees  respectively  in 
the  yaw  right  direction  from  the  servo  at  Xp. 

Hie  displacement  of  the  pushrod  of  blade  n,  X.  n,  is  defined  in 
exactly  the  same  manner  as  in  the  tail  rotor  control  ‘system.  Equation  AU6 . 
The  displacement  of  the  point  P on  the  swash  plate  is  derived  assuming  that, 
in  addition  to  its  cvn  dynamic  displacements  Xj,  , X^  , and  X^  , the  svash 
plate  also  follows  the  motions  of  the  hub.  * 

With  these  displacements  defined,  the  energies  associated  with  the 
main  rotor  control  system  can  be  shown  to  be  as  fol1  jws: 

Total  Kinetic  Energy 
N 


MTOT 


y £ 

♦ W^anA  - V * E 

“ 8*i  na 

+ /,2Rq^2\  " *A  ~ + £ ^0Y,sqs^ 

'S  NA 

+ V^A  + V + E (*Z,sqs}]2} 


NA 


(A5»») 
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Total  Potential  Energy 
N NA 

!‘(^LK1>„P1.n-  |:Uz,„V  -yXA  + XF) 

- ^/RS(1/2(XA  - XF><co8,nooB!n  - sln*nsin«n)) 

- / - XA  - XjXelD*.nC086n  ♦ cos¥n8in5n))]‘ 

D 

* V*F  - V"  * KA<XA  ’W' 


* kl(il  - V > 


(A55) 
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where 


*n  - fit  + 2ir(n  - 1)/N 


(A56) 


end 


6 

n 


FS 


+ 6 


PR,n 


(A57) 


Total  Dissipation  Potential 


NA 


MTOT 


« 5lCl»nCil.n-  £ <*Z,s*s>  ' * V V 

- (1/2(Xa  - Xp)  (cos^cosd^  - sini|/nsin6n)) 

- ^/p  (1/2(2Xl  - XA  - Xp)  (sini/»ncos6n  ♦ cos^sinfi^ )]“ 

s 

* cF(Xp  - X0F>2  ♦ oA(xA  - Xqa)2 

* CU,(iL  - W2’ 


(A58) 


Servo  Systems 


Only  servos  with  first-order  transfer  functions  were  considered. 

Main  and  tail  rotor  servos  were  assumed  to  be  governed  by  the  same  equation 
of  motion;  however,  they  Lay  all  have  different  dynamic  properties.  The 
model  used  to  describe  the  servo  system  is  shewn  in  Figure  Ak.  In  this 
model  the  displacement  X,  in  the  case  of  the  tail  rotor  control  system, 
corresponds  to  X5,  and  K corresponds  to  K5;  see  Figure  A 2.  Similar 
analogies  can  be  made  for  the  forward,  lateral,  and  aft  main  rotor  servos; 
see  Figure  A3. 

From  Figure  Al»,  we  can  see  that  the  displacement  of  the  valve  spool 
relative  to  the  servo  housing  is  cy  = X - X . This  can  be  written  in 
the  general  form 
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VALVE  SPOOL 


t 


r'igure  kk.  Control  System  Servo  Model. 
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U,  -HL.  + RX  + RX. 
* 1 ” 2 0 3 c 


(A59) 


where  the  servo  linkage  ratios,  Rj  , R2  sad  R , depend  on  the 
of  the  servo  linkage  arrangements.1 

Plow  through  the  valve  Is  given  by 

% - Vv-Vp 

■ VXC  - iO>  ♦ CLSSP  + (V*"s><P 

and  since 

KC(XC  - V * VP 


giving 

XC  - VP/KC  ' XAF 

and 

*C  ■ VP/ltC*XAF 


it  can  be  verified  that 


C R 


X,—,  + C R Xn 
1 vF  q 2 0 


C R XA„ 
q 3 AT 
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kinematics 


(A60) 


(A61) 


(A62) 


(A63) 


- [«ts  ♦ 'v>/*P  - <vvkc:ivp 

- vX  ♦ yfcllgA^ApJp  - 0 
If  ve  express  ty  in  the  fora 

eV  ■ V \ *VF  + *0  + *\V 
2 2 

and  let  R ■ oQ  such  that 
2 S 

6V  ■ °S(  ‘/«gX*F  tX0  + "’Vc* 
then  using  this  fora  for  ey  and  the  Identities 


■ t(CLS  + CP)/C,*P  * Rj«(J  (*6T) 

U2  - [l/fC  ♦ Vj/llNg*!]  (A68) 

together  with  the  equalities 

Ap«P  - - K(X0  - X)  (A69) 

VR  - - K(X0  - X)  (A70) 


in  Equations  A60  and  A6U,  it  can  be  shown,  after  some  algebra!'*  manipulation, 
that  the  servos  are  governed  by  the  equation 

(V»,)(VV,,(1  * * <V»,  ♦ K>x0 

- <as/ui><V°sVtV  - « + (''S/",)(V°3)XVF 

* (Vl‘1)(R3/t,S>XAF  * <“S/>‘1)(V'>SC,)XAI'  ’ 0 (*T1) 


(*6M 


(*65) 


(A66) 
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Here,  (ag/y^)  is  the  "Servo  Power-On  Stiffness  With  Locked  Input," 
(Ap/agCq)  is  the  "Servo  Time  Constant  or  Inverse  of  Velocity  Gain,"  and 
(14)  is  the  "Reciprocal  of  Coufclned  Actuator  Hydraulic  and  Structural  Series 
Stiffness ." 

In  this  analysis  Xyj,  and  were  not  considered  to  be  independent 
degrees  of  freedom.  Rather,  they  are  assumed  to  be  the  result  of  airframe 
vibration.  could,  for  example,  be  construed  to  be  the  result  of  vibra- 
tion at  the  pilot  station  being  transmitted  through  the  control  system, 
whereas  gives  a measure  of  the  effect  of  the  servo  support  vibrations. 
Since  accurate  measures  of  these  quantities  are  not  always  available,  it 
was  considered  that  they  could  best  be  represented  by  expressing  them  as 
functions  of  the  airframe  motions  at  the  hub  using  feedback  factors  for 
the  support  and  valve  to  define  relative  participations  as  follows: 


(8FV,X*X,-*s  + SFV,Y*Y,sqB  + SFV, Z*Z,sqB 

+ 8FV,ex*0X,sqs  + SFV,6Y^0Y,sqB^ 

NA 

J5.  ^SFa»x*x»8<18  + SFA,Y*Y,s^s  + SFA,Z^Z,b<1b 

+ SFA,0X^0X,sqs  + ®FA,0Y^0Y,s^b^ 

NA  ^ 

£ (SFA,X*X,sqs  + SFA,Y^Y,BqB  + SFA,Z*Z,sqs 
s-1 

+ SFA,0X*0X,s*g  + SFA,0Y^0Y,BqB) 


(A72) 


(A73) 


(A7»*) 


These  expressions  are  substituted  directly  in  Equation  A71  to  give 
the  fined  form  of  the  servo  equation. 


Final  Font  of  Dynamic  By  ten  Equation! 

Final  foraa  of  the  dynamic  system  equations  of  motions  are  obtained 
by  formally  expanding  all  of  the  total  kinetic  and  potential  energy  and 
dissipation  potential  expressions  derived  in  the  preceding  sections,  and 
then  performing  first-order  perturbations  of  all  of  the  generalized 
coordinates  about  specified  steady  initial  values.  This  procedure,  as 
mentioned  earlier,  requires  an  extensive  amount  of  matrix  manipulation. 
However,  the  methods  employed  are  veil  known  and,  provided  accuracy  is 
maintained,  there  is  no  problem  in  principle  in  arriving  at  the  final 
equations . For  this  reason , none  of  the  expansions  are  repeated  here . 

The  equations  that  follow  are  very  lengthy.  Therefore,  to  improve 
readability,  they  are  all  written  in  exactly  the  same  format.  Also,  since 
two  types  of  control  system  are  treated,  their  contributions  to  the 
equations  are  vritten  separately.  Finally,  the  servo  equations,  as  applied 
to  each  type  of  control  system,  are  written. 

The  first  block  of  equations  contains  all  of  the  terms  that  arise  from 
the  coupled  rot  or/ airframe  system.  It  also  includes  all  of  the  terms 
contributed  by  the  motion  at  the  pushrod  as  defined  in  Equation  AU6 . Each 
equation  first  defines  acceleration  coefficients,  then  velocity  coefficients, 
and  then  displacement  coefficients. 

In  every  equation,  all  generalized  coordinates  and  physical  properties 
associated  with  the  rotor  system  have  a subscript  n,  the  blade  number. 

This  subscript  is  omitted  from  the  equations  without  ambiguity.  Also,  to 
give  added  clarity,  the  equations  are  presented  in  integral  form. 
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Airframe  Mode  Equations 

* ♦n.A  * rU>  + * rt.]) 

* ♦z.j^z.i  * ♦ex,i[(*2  - rYo),ln*  + <•  ♦ r + Vo  - Vo)co“*l 
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* ♦z.i^6  + r + a2Yo  “ ^V008*  * (*2  ~ ***0  )®in'^ 
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+ b2(b2  + 2r0o)])}+  Ixdr  { ♦0Y#J4,EO(“  ♦0y>ico80o8ina|,  - ^^0080^082*) 
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* ♦ex,i'“<l,]!Oco,eo  * *,roBlnflo),lB2*)  + ♦eY,j(*eY,i(co82*  - Vin2*} 
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“T7"  ~4Z^ m «CT7 v . ■ TT^ z — 


♦ *'?0tia*o)9in2*  * ♦ex/’  q,EOco§eo  + «l,FOBineo)coB2*]) 

♦ ^♦eY.j^ey/"  ^ 0 C0B2*  ~ 8in2^  + ♦ex,i^“  c0 B2*  + 2Y0Bin2'l')] 

♦ ♦0x,j^eY,i(“  c°82*  * 2Y0Bin2*)  ♦ ♦ex/2yo0082^  + Bln2*)])J 

♦ Izdr(211<l-roC*eifJ[»eIii»in60c0.2*  - ♦ex,1»i"«„«o»2*]  * Vjt 

- ♦eYtioin0o8in2*  - ♦ex,!810^0082^  + n(^0Y^[4eytl(sin2eo8in2i|(  + (2yoBin20( 

+ 0 sin28  )coa2i| i)  + . (sin20  cob2\P  - (2y  sin20  + 0 sin20  )sin2^) 

0 0 0 A 1 1 0 O O O U 

♦ ♦ex,iCO820o  ] + ♦ex.J  [*0Y,i(Bin20oCO8211'  -(2YoBin20o  + 0o8in29o)8in2’l') 

+ $flv  .(-sin20  sin2ip  - (2y  sin20  + 0 sin20  )cos2<p) 

0X,1  o o o o o 

♦ VtCOs2io  ]>  >]  ! * I 2tA.A,A,J  ’j  | 

+ j E^R_e  mdr^2n  (♦x.J  ((b2Y0  + a2$0)c"S<'  + b28in*> 

+ ♦Y./Vo  + *2*o)Bin*  - V0**  + ♦eY,J(b2^t2 
+ regain*  + Wo  + Vo)c08,(')  * ♦0X„J(_b2(b2Yo  + Vo*81"* 

+ b2(b2  + r0Q)co8ip))}  + Ixdr{f3(^0Y^[co820o8ini);  + Cyocob20q 

- 0o8in20o)cO8*]  + >0X  j[cos20ocos^  - (yocos20o  - 0osin20o)sin^  ])} 

+ IYdr{nUey  j(y0co8i|»  + simp)  + <pQX  j (cosip  - y^irnp))} 

+ I dr{fl  (<pnv  .[-coa20  sirup  - (y  oos20rt  - 0 Bin20  )cos\p] 

2 01  o o o o o 

+ *0xj[  " cos20ocoBtp  + (yocos20o  - 0o8in?0o  )sin<p]) } J J 

♦ lit*"  6 mdr{2fl($x  j(b2  + rB0)c°sip  + $y  j (b2  + r0o)simp 

+ ^Y,?2tb2  * 2rB0)c°8*  - ♦0X,Jb2(b2  + 2r0O)8in’,')} 

+ Ixdr{n((P0y  jCOB20Qco8(p  - $QX  jOos28osinip) } + Iydrin($0y  ^ccsip 
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- ♦eXij.l-*»  + Izar(“C  - * ♦ex„)c"'2Vln','))]  * | 

♦ j f;[^-emdr{20(*X|  j[  - (»2  - nro)co8i|»  ♦ Ca^YQ  + r)sini|i] 

+ *Y,j£  " (Y°  + r)C0B*  " (a2  " rY0)8in*^  * *0Y  J[(r2fi0  + b2(V° 

+ r))sin*  - (b2(a2  - ryQ)  + a^cos*]  + *0X,jfCr2Bo  + ^(a^  + r))c°M' 

+ Cb2(a2  - ryQ)  + a^a^sin*]) } ♦ Ixdr{2nq'EOUeY>JBin0osiii<|» 

+ *0X  j8in0oco8*)  + n(*eY  jE»in20ocosf»  - CYoain20o  - 60)Bini|»] 

+ " 8in20o8intp  - CYosin20o  - 6q )coa<»3 ) } 

+ Iydr{2n(«eYJ(  - q’EOain0o  - 4,poeo80o)8in#  ♦ 

‘ <l,FOcO80o)co8*)  * n(^0YfjBo8ln*  “ ♦0X,jeoCO8,|,)} 

♦ + ♦0XJCO80oCO8*)  + flC*0Y,J(  - 8in26oC^ 

+ (Yosin20o  + 3Q)sin^)  + *eX,j(8in2e08in'J'  + (YoBin20Q  + Bo)cos^)}]  y | 

+ {R"endrf2fl{^,J^E,i(YoC08*  + 8in*)  + vi.iCO“’*' 

♦ ^F,i0oC°8^  + *Y,j£“  ♦E,iC<-08’1'  " V1^  * vl>i8in,|' 

+ ♦F,i0o8in^  + ♦0Y,J^E,i(b2  + reo)8in’»'  + b2UF,ieo 

♦ ♦ ♦ex.jt  - \L*T9i*o  + ♦E,i”o)8in*  + ♦E,iCb2  + rBoW'"] 

+ ♦w.jVV*  ' ♦ex,jvi,ib28in*)}  + ^^y.jK^V08* 

♦ Cq’E(^in20o  - YoBin0Q  ♦ ^cosO^Bin*)*’^ 

♦ (CO80QCO3\()  + (q'E0  - Yocob0o  ♦ Bo8in0o)siru(t)*,F^i] 

+ ♦oi.jft  - «in0o8in*  + U,JOBin20o-YoBin0o  + f^ose^coB*)^ 

♦ C - co80osin*  + Cq’K)-Yoco80o  + B^^Jcosij;)^^])} 
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♦ IYdp{n( ♦ex , J ^E,i(  ■ 8ln0oCO8*  + C " 2q,F0CO826o  ’ q,E08in26o 
“ 0OCO80O  + Vin0o)8in*)  * ♦,F,/CO80ocO8’,,  + ( • q,FX)8in20o  + q'E0CO826o 

- eoBineo  - Yocoseo)flini(.)]  ♦ ♦0XtJ[4»EjitBineosini|/  + ( - 2q'FX)co820o 

- q'EOBin26o  - 0ocob0o  Yo»in0o)coB*)  + *'p  J(-  coB0osin*  + (q'EC)oos20o 

- q,poain20o  - 80Bin8o  - yocos0o)cob*)]) } + Izdr{n(*QY  j[(-  sin0Qcop* 

+ C2q' eocos20o+  Yosin0o  + Boco80o)Bin*)*'E  ^ + ( - cos0ocos* 

+ (q'FX)8in20o  + YoCO80o  + 0o8in0o)sin'l')*,F,i^  + *0X,  j K8in0o8in* 

+ ^2q,fOCO820o +Yo8in0o  + e0C08e0)c08*)^E  i + (co80oBin*  + (q'FO8in20o 

- q,EOcoB20o+Yoco8  9o  + 8osin0o)  cob*) 


1 N IA  r p . 

* EEj{RW^  ♦ [♦„/- 

>n*l  i=r 

- b20o)cofl*)  + *ey  - (e  + r + a^  - b^Jsin* 

* U2  ' n,0)c0a*)J)]  h j * j Vj^A.^j  ( 

+ j “ ^b2Y0  + a26o)8in*  + b2c°8^ 

+ <frY  ^[(b  Y-  ♦ a 3n)cos*  + b sin*]  + [Ca  (a  y.  + r + e - 2b  B) 

I,J2°2°  2 0*tJ  2 2 0 2° 

- b2y  Jsin*  + ( - a (a  - nr„)  + b (b  + rB  ))cos*J 

2°  22°  22° 

+ *fl  [(a  (a  - ry_ ) - b (b  + rB  ))sin*  + Ca  (a  y.  + r + e - 2b  8 ) 
y^»J  22  0 22  0 22°  2° 

- b^Yo)cos*])}  + Ixdr{n2q'EC£*0yjjCOS0osin*  + *0x^cos0ocos*) 

+ fl2[*ftV  ,(cos20  cos*  - (y„co820„  - B sin20„ )sin*)  + *ftV  .(  - cos20  sin* 

HI  0 0 0 0 0 OA,J  0 

- (yocob20o  - eosin20o)coa^)]>  + IYdr{n2[*0y  (-  q'EOcos0o  + q'^s^  )sin* 
+ ♦er,<j("  q’Eocos0o+  q,F(fineo)co8*:,}  + Izdr{n2q,F0(_  *0Y,j8in0o8in* 


(a^  - nr0)sin*  - (e  + r + a^Y0 
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- ♦0x^aineocos*)  ♦ 02[  ♦gyjC  - cob20ocob*  + Cyocob20o  - 0osin20o)Bimp) 
+ *0x,j^  wOb2  {y^coB2()^  0^8in20^)co8|p)3}  J ^0®^| 

+ | /^[i£R  9 ■drtn2C  " + r6o),in*  + ^Y,J^b2  + re0)CO8’J' 

♦ ♦©Y.j^r  + e)  (Vo  + r " b26o}  + r(-Vo  " Vo*  " V2r60  + b2^Bin^ 
+ [ - a2^a2Y0  + r “ Vo*  " b2^b2  + ^o*  + r2Y0Jc03lM 

+ ♦0X,j(ta2(a2YO  + r ' Vo*  * b2C2r60  + b2*  - r2Y0l,in* 

+ ECr  + e)(a2YQ  + r - b^)  + r(a2Y0  - b^)  - b2(b2  + 2r0Q) Jcos*)]} 

♦ ■cxdr^2<1,EO^“  ^0Y,JCO88oCO8^  + ^0X  jcos®0sin’>)  + ft2E^0Y  j(sin20osin^ 

♦ Yoco820Q cosip)  + $0J  j(sin20o  <,osip  + Yocos20osirup)]} 

+ IYdr^2^0Yfj^,EOco80o  - q ,posin0o)cocV  + *0XJ(-  *'mcoa*0 
+ <l,FOBin0o)8irnp  + ♦qy  j(  - Yocosip  - sin<>)  + ♦qxj (Y0sinip  - cos*)]} 

+ Izdr{Q2q'FO(^0Y^8in0ocoBip  - <p0x^sin0osin<p)  + fl*E  4>0y  j (cos20osinip 
+ yq  8in20ocosip)  + $.,.v  j(cos20ocos^  - Yosin20osiiup)]}  J 0 | 


iir 


mdr{  n2E  ( (a2  - rY0)sinip  + {&2',r  + r)cosip) 


+ ♦yj  E(»2y0  + r)sinip  - (a2  - rY0)cosip)  + <p0y EM®^  - b2YQ) 

- VVo  + r*  + a2^  r + e3®0  + V " r2eo)8in*  + ^b2(a2  " n0) 

+ b2(a2Y0  + r)  ‘ Wo  - re0}  + r2Bc)c°8^  + *0X,J[(-b2(a2  “ ^0  * 

- VVr  + r*  + Wo  _ rPo*  " + (r(a260"  Vo)  " b2(a2Y0  + r) 

+ a ([r  + e]0  + b ) - r20  )cosip]]} 

2 02  0 
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•r 


♦ Ixdr{fi2q'E((^eYtj8in0ocoB(|»  - ♦ex^«ineosin*)  ♦ fl2[  ♦^(fl^osZe^ou* 

- feln20  slnt)  + $ ( - 0 cofl2e  Bintp  - ^in20  coa^)  J) 

O 0.1,J  OO  0 

+ IY4r{n2t*8Y.J(-  »'E0,ln9o  - <»kf0”eo)co“*  * 98X»J(^'jO,ln9o 

+ ^'F0co89o^8^n9  - ♦sr./o0089  * 9ex,j9o8l°9^) 

♦ izar{0Vro(*ej(Jco»e0co.*  - ♦ei>Joo«80»in*)  ♦ n2[*6Y ^Catose^ia* 

+ B08in260cos*)  + 8flI  , (*5Bln280coB*  - 0oain28oBin8)]}J  y | 

i it  trn  ' 


1 n=*l  i 


{E"8  «»lrta2[*Xi,(  - BE  , (Y0Bln*  - coa»)  - Vj  jSin* 


- ♦Pil0osln(/)  + ♦yjj*9e,1*81“*  * x0cos8)  + v^jCOBB  + 8p  ^cob*) 

+ ♦«Y,J(UF.l(Vo  + r * e - ai#1>1  ‘ ♦E,lVo)8in9  + <-  ♦f.i'S  ’ 'V 

+ * r9o))c°89)  + ♦0x,j((9F,l(a2  ' ' ♦E.i^z  + r9„>)8l"» 

* <♦  f,i(V»  * r * e " "V©1  * ♦E,iVc)c08'l'>  ‘ VjViV1'^ 

* ♦ex./l.^a'08^1  * IXdrt°2[96Y,j<  - 8ln9o8ln*  - (Vln0o 

- 0ocos0o)coB^)(|>'g^i  +4>eXj1(  - Bin0ocosi(;  + (YQsin0c  - &ocos0o)sini|i)$'E  i]} 

+ Ixdr{O2[4eYj(#'B^8in0osin*  + (-  q'E()sin20o  - Bocob0o  - q'FOcos20() 

+ y 8in0  )cosiJ))  + $ » . ( - cos0  sin^  ♦ {-  q'r,rtsin20rt  - B sinO^ 

0 0 Y , X 0 rv  0 0 0 

* ’’E0cos28o  .yocob60)cob8))  ♦ ♦0XfJ(8'E>1(BineocOB8  * 

+ 0qco80o  ,-»■  q,FQco820o  - Yosin0Q )sirup)  + ^(-  cos0ocost|i  + (-  q'EOco820o 

+ 6Q8in0o  + q'pQ8in20o  + yocos0q  )siiu(>) )]}  + Izdr{n2[$0y  C(cos0o8in<J> 

+ C-  q,EOc°820o  + Yocos0o  + PoBin0o)co8*)f.pjl  + ♦,Efiq’rocos20oco8*) 

+ ♦ox  j((cO8  0oCO8j)+  (<1,eoCO820o  ” YOCOB0O  " 0osin0o)sin’<')4’,F  i 
”*,E.i<l,FOC0828o8in^}]<lr,i|  “ ° 


CA75) 
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'r>>YKT  . ~*C r»5W."  ' 'V  Tj  "~Tf*  ?*”' 


Blade  Pitch  Equation* 

S[<!'"**4r<*I,lCtVo  * *2#0!*In*  - i2C0,*>  * *Y,1(  * <Vo 

♦ a20o)co*i|>  - b2»lni|i}  + a2^z,l  + *0Y,i^a2  * b2^Y0  ♦ a2(r  + e))aln^ 

♦ ( - a2(a2  - ry0)  - *2(l>2  ♦ rBQ))coe^  ♦ ♦ex,it(a2(a2  ’ ^ 

♦ b2^b2  + r0Q  ))■*“♦  ♦ C(®2  + ^|)y0  + a2(r  + e))co8\|»]) 

♦ IXdr{<1'B0(*0Y,iCO,eo,ln*  * ♦6X,lC0,e0e0,l|,))  + IYdrIV.lC",,B0CO,eo 

♦ 4'ro“ln6o,Bl”'1’  + ♦ei,i(  - «’jnco,eo  + 4Vln8o!coa*  * ♦eY,iCYo“in* 

- cos*)  ♦ ♦eXfi(Y0C08t(»  + sinifc)}  + Izdr{q'F0(  - *qY  jSinO^in* 

- ^gx,ib*qBoco&b^  J ♦e«i  J 

♦ [^11.’  ♦ if))  ♦ iYarU>]  *J ®T j •»  *,*!♦, J„ST  j 

X>2 WQ  + a2r>  + Xxdr{q'B3coB6o)  + Iydr{  - q'BOcoB80 

♦ ,'ro,ln6o  ♦ V 

♦ IZ'5r(  - 1'pO,1"°o)]  V | ♦ j - M,lft“53*8PRB  j 

♦jtf-  'emdr{b2r}  + Ixdr£q,EDBin0o}  + lydrf-q'^sii^  - q'^cos©^ 

♦ I2Br ^^FOCO*Bo ^ ] V * | - MlLitlMl*6PBY  | 

+ llS[^Wfa**»'‘  * Ve-1>  * * IYdr(  - 

+ IZdr{<l,po*,E,l}  J *0  ^T.i  \ + ^^l^^FPR.l 

♦ tL2t«»1/(B-«ij*ETil))qepR  qf>i  | 
iNA  T R-e 

■eodrJ2fl[^0y>1(a2(a2  - rv^sin*  + a2(a2y0  + r + e - b20Q )cos<J>) 
+ ^0X,1  ^ ” a2^a2Y0  + r + * - t>260)sin^  + a2(a2  -ryo)cos<>}]} 


I 


,>+<-,  r^-.frVjTfc'c,  sn-.  ** 


"V  ■ ■ 


1^(211(1 'EO[*8Y>1co.eoco.*  - ♦ex  lco.e0.in*]  ♦ a|>M,i<-  co’29o,ln* 

- (t0co*260  - B0.in2e0)co.*)  **#Xji(-  CO,260COB*  * (V°‘2eo 

- »o.in2eo).in*)])  » iYar{o[«eiil(-  * «'ro8in9o)2c°8* 

* 89x,i<i,'eoco,9o  - «’ro,ln6o>2*ln*  * ♦8Y,i(To“'*  * 8tn*>  * ♦ex.i'-  Vin9 

+ co.*)]}  ♦ Izto(2ft1'K1[-*eYjl.lne<)co.*  ♦ *eXjl.ine0.lD*] 

< n[*9I  1(co.260.in*  + (Y0C0.260  - 60«lc290)co.*)  ♦ ♦eXjl(co.2e()co.* 

- (yocos26o  - 0osin26o  Jain*)]}]  J + j (2C0I^n  + CjL2*8pr^t  [ 

[^R‘emdr{2n[-  b2(b2  + reQ)]}  + Ixdr{-  ncoa20o}  + Iydr{-  ft} 

+ Izdr{ncoa20o}]  | + j - J + ^R"emdr{2nta2(b2  + r6o)]} 

+ Ixdr{-  flsin20o}  + Izdr{nsin20o}]  ♦qY  | + | - C^^tana^^RY 

+ j E [iR'e"flr,2“(-  ♦f.lVo  * ♦*,1*2®. ’ ‘ 20Ul.lV  * IXdr<“[-  ♦'E.i81"6. 

-V  00«60]J  ♦ Ixdr{n[*'E  i.lne(>  - ♦'F>i=o8«0]>  * Izdr(“[*'E,1“ln8o 
* . i*  NE 

* ♦’f.!0089.^]  ♦e4p.ij  i^(*  ClL2(9FFB,i 

[L2t«K.1/(R-e)]*ETil))9ePRir,i  | + j [{"-■arfrfWV o - b28o) 

(b  ♦ re0)  + .2])  * Ixdr(-  n2co.2eoJ  * lzdr{a2eo.29(>)  ] *29t  J 
( iH_*K+e2  )eT  i ♦ I K,L2*28PB6T  I 

^R"Cmdr{fi2[a2(a2Y0  + r + e)  - b^b^  “ 2a2e0^}  + IXdr{n2^q'E0CO8&o 

- Yocos20o]>  + IYdr{Q2(-  <1'eocob0()  ♦ q'FOsin0o)}  + Izdr{n2[-  q'p()sin0o 

+ yocob20q]}]  $06  j + - K1L2tan63*0PRB  | + j \_l  mdr{n  ^2  + V0o 

+ ebj  } + Ixdr{-  n2eooo820o)  + Izdr{n20oco820o}]  *0Y  J 


+ 

- b. 


" 

I 


137 


r»:C-»'-*  .—  *—*r— / — «>■■«>  «rc-- 


i'rjT»-<r 


tD 


♦ j - K^«»lWj  ♦|^*--r{^[*F>i(.1  ♦ «0>  - ♦1>1(b2 

♦ (r  ♦ «)®0)])+Ix4p{-  G^^’j^cos^}  * IYdr*°2eo^,E,icol,eo  * ^'p,i,in6o^ 

♦ I2dr{02[.  ♦,g>1«'wco*290  ♦ (l'E0co.2eo  - 60«ir.eo)*-p  ]>]  ♦,  q,  J 

i HE  ' 

♦ L?  (-  “.VW  * £Vw1/<r-*)]W>>Wt,1  } 

Tail  rotor  control  system  contribution: 

♦ C,V«A  ♦ KlL2*8PRX2 

Main  rotor  control  system  contribution: 

♦ PsC,L  Vp  (-(sin«  + cosi)sint|*  + (cos6  - sin6)cos^) 

1 2 Hg 

♦ * a>«a  * PsCjI^^B/g^Usinfi  - cos6)sin^ 

- ( cos 6 + sin6)cos^)  + ^L^gppXp  + [(^L^/^cosfisin*  + sinficos^)]*^^ 

+ The  coefficients  of  XA,  Xp,  and  X^,  which  are  identical  to  those  of  iA, 

Xp,  and  X^  respectively,  with  replaced  by  K . 


- 0 


(A76) 
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• * trm  ■*’ ' * . >«*  -k  • 

: ✓ .-«*  t 


i 2’ ' ' ** 


Blade  Rigid- Body  Flapping  Equations 


lS[r^^2  * rSo)Bin*  - + re0)c08*  * *Z,i(r  + Vo 

• b20O>  + ♦«,lI(e(Vo  ■ b2fiO}  + r(r  + * + 2Vo}  + b2)8in*  + (“  a2(Vo 

- \»0)  - r(*2  - ry0))c08^]  + ^[(a^a^  - b^)  + r(a2  - rY0))sin* 

+ (e(a2Y0  - t>20o)  + r(r  + e + 2*JQ)  + b2)cos<|/]}  + Ix<ir{q'E0( 

- ♦ey  icoseocos^  + $0X  icos6oainij;)  + $0Y  ^ (cos20osin<|/  + (yocos20o 

- *58  sin20  )cos<jO  + $0V  . (cos20  cosi|>  - (y  cob2©^  - *s8  sin20  )sin^) } 

O O UAfX  O 0 0 0 0 

+ ly^^OY.i^'EO0080©;  *'FO8in0o  “ Yo)c08* 

+ *0X,i(‘  q,EOCOS0o  + ^,pO8in0o  + Yo)sin*}  + ^drCq'^C^^sinO^oa* 

- $av  .sin0  sirnj/)  + $ov  . (sin20  sintji  + (y  sin20  + >58  sin20  )cos)J)) 

vA|l  0 vl  )1  0 OOOO 

+ (>0X  i(sin2eocosi|<  - (Yosin20o  + ^8Qsin20o  )sini|i)  jjc^  j 
+ j ^R-emdr{((a2  + b2b0  + a^)}  + I^fq'^cosO^  + IY<lr{(-  <TEOeos0o 

+ <l'FO8in0o  + Yo}i  + IZdr{"  q,FO8ln0o}]^|+  | • ^“VePR^rj 
+ | ^R_emdr{r2  + b2  + 2a2rY0}  + Ixdr{cos20Q}  + lzdr{sin20( 

: ![<■-• 

>]vj  * jM^tan^tana^'  J 

* r)  * ♦e.iVo  * v,jl»2>  + IXdr(«'F,icosV  * IZdr(*  *'E,lslneo)]«T,lj 

( NE  * „i 

+ jEi(M1Van63(*FPR,i  + [L2taxwl/(R-®)3*ET,i))4p,i  [ 

+^^R"e,ndr{2n[^ey>i((r(r  + e + 2a2Y0)  + ©(a^  - b^)  + b2)cos* 

+ (a2(Vo  ‘ Vo}  " r2Y0)8in^  + ♦ex#i<<-  ^(r  + e + 2*^)  - efa^ 


+ | M^L2taa2i3B  | + 


+ Izdr{-  }83in20o} 


emdr{-  b2(a2  - ryQ ) } + Ixdr{*ssin20o} 


'mdr(*F,i(Vo 
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- b20o)  - b2)ain*  ♦ (a2(«2r0  - b^)  - r2Y0)co8*)]} 

+ Ixdp{20<i,Bo(*flY,iCOBeoBin*  * ♦ex,iCOBeoCOB',') 

♦ fi[#0Yii(cO8<>  - 2(yocos20o  - J56osin2eo)«in*)  + ♦0Xtl(-  sin*  - 2(yocob26o 

- Jrfo8in28o)coB*)]}  ♦ IYdr{fl[*0y^1(-  q'EOcos0o  + q^Bine^ain* 

+ ♦ex,i(-  q,EoCO80o  * q Vin0o)2c08* 

♦ ♦eY,i^2YoBin’1'  " C0B^  + ^X.i^e008*  + ♦ Izdr{2flq'ro[ 

‘ ♦0Y,i8ln0o8in<’  “ *0X,i8ineoCO8^  + n^0Y,i^CO8*  " ^o*in2»o 
+ 0oBin20o)sin*)  + *0X  A(-  sin*  - (2YQ8in20o  + Bosin20o)coB*)]}]q^  j 

+ j ^R"eodr{2flb2(b2  + r6Q ) > +Ixdr{ficoa20o}  + Iydr{n}  + Iz<lr{-  ncoB20o)J*0&,Ij 

+ j - C1L2tan«3*0pR  $TJ  + j C2L2tan263  0 J + j ^R-eodr{2nCb2(a2Yo  + r) 

+ r20Q)}  + Ixdr{fl(2q'EOain0o  + 6Q  - Yosin20o)}  + Iydr{n(-  2q'E()sin0o 


- 2^,FOCO80O  - 60)}  + IZdr{n(2q,FOC°80-  + YnBln20n  + 0n)} 


■M 


+ C L2tan6  tana 
i 1 2 3 


W«a[*E>i(b2  ♦ reo)]) 


M ilr 

+ Ixdr{fl[**Ejl(q,EO8in20o  + 8ocos0o  - Yo8in0Q)  + *,F,i(q’EO  + eo8in0o 

- YoCO80o)J}  + IYdr{nC*,E,iC"  q,EO8ln20o  " 2q,FOCO820o  + Vin0o  " 0oCO80o) 

+ *'F,i(<l,EOCO820o  “ q,FO8in20o  - VO80o  " 0o8in0o)J} 

+ IzdrM*,Etl(2q'FOcoB20o  + Yosin0Q  + 6ocos0Q)  + ♦ ,F,i(q,FX>8in20o  + Yocos0o 

+ 0oBin0o)^j^T»i|+j^CiL2tanfi3(^FPR,i  + tL2tanai/(R-#)^ET,i))qTfi| 

+ j ’^R'emdr{fi2E*0Yji(-  b2(b2  + 2r0o))cos*  + ♦0XflO>2O>,  + 2r0o))Bin*]}]qi  j 

+ j1  j^R_emdr{n2[a2(a2Y0  + r + e)  - b^b^  - Sa^)]}  + Ixdr{n2[q'E()cos0o 

- Yoc°b20o+  0oBin20o]}  + Iydr{n2[-  q'E()co80o  - a'FOsin0o]>  ♦ Izdr{Q2[q'p()8in0o 


140 


+ cob  20  - 0 sin20 

oooo 


* - kjl^3Wt 


* V 


j [f R-*“dr ®2-r0  - »260)  ♦ r(r  ♦ e)  - l*reot2  ♦ 2« 2n0  - b*]> 

lvdr{fl2flin20  } + Ivdr{-  fl2}  + I„dr{fl2coB20  }|0  f + j K L2tan26  0 
X ol  z oj  ) i 12  3 

+ | j^I<“emdr{n2rb2Ca2  - ryQ)  ♦ *^Q{2r  + e)]}  + Ixdr{ 

- %fl2sin20o)  ♦ Izdr{^2ain20o}]yJ  + | K^L2tan6  ^tano^y 

♦ Izrwt-  »E,ib2Yo  * ♦ r ♦ • ♦ 2Vo>  - v.,ib2]} 

'i^l 

+ Ixdr{-  n2^'EjlBin0o>  + IYdr{n2[^'E>i8in0o  + ♦,F,1cos6oJ} 

+ Izdr{-  «2U,F,ic08eoJ}}lr^[+  j E^KlL2tan63(*FPRti 

* [V“ai/(R-e»ET,l,),T,l  j 


Tail  rotor  control  system  contribution: 


-CL  tanfi  X,  - K L tani  X 
1 2 3 2 1 2 3 2 


Main  Rotor  Control  system  contribution: 

+ [-  ^Cjl^tanfi  3^B/p^(-(sin6  + cos6)sini|>  + (cos6  - sin6)cosiJ>) 

- 'sCjI^tanfi  ^ ]iA  * [-  ^I^tanS^/^UsinS  - co86)sin4»  - (cos6 

< sin6)cos<|;)  - JjC^tan^ ^Xp  +[-  CjL^tanfi^/j^CcosSsiniJf  + sin6cos\|/)]XL 

♦ The  coefficients  of  XA,  Xp,  and  XL,  which  are  identical  to  those  of 
XA,  Xp,  and  ^ respectively,  vith  C replaced  by  Kj . 


(A77) 
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•tc-« — — 


yr1" 


BIMe  Rigid-Body  Lagging  Equation! 


ft  R“®mdr{*XtlC-(a2  - rY0)«in*  - (»2yo  ♦ r)cos*)  + ♦y.I^Vo  + r)Bin* 
♦ <a2  - rY0)co«*)  ♦ iz%i*2»0  ♦ ♦0Y,i((Ce8o  ” b2](a2  " rYo))ain’1' 

+ ("  b2(Vo  + r)  - (a2  + r2)fi0)  ^ + ♦eX,it(b2(Vo  + r)  + (a2 
+ r2)0o)sin*  + ([e0o  - b2](a2  - rY0))cos*)>  + Ijdrfa'ggC-  ^ i8in0ocosi(» 

+ ♦0XtiBin0oBin*)  ♦ ♦eY^1(Jsain20oBin<»  + OsYosin20o  - 6osin20o)coB<i) 

+ ♦0X^1C%Bin20ocosi(»  - (%Yo8in20Q  - 0oBin20Q)Bin*) } ♦ IYdr^eY,i^q,EOsin0o 

+ q’poCOBOjcos^  ♦ ♦ex>1C-  q'EQ*^  - q,FOco80o)8in</}  + 

- ♦qy  iCob0ocob^  ♦ ♦eXficO80o8in*)  + ♦0Y|iC-  158in20o8in^  - O5YoBin20o 

+ Bocob20o)cob<»)  + 40x  1(-  %Bin20ocoBi|)  + (%YoBin20o  + &ocoB20o)sin\Ji)}j  ^ J 


mdr{b2r}  + ^drfq^BinO J + IY<lr{-  q'EOsin0o  - q'^cos©^ 


+ IZdr{q,FOCOS0o) 


(]Vt|+  j- V^aWT  ]♦  Ik 


R-e 


mdr{-  b2(a2 


- rr0)}  + Ixdr{5apin20o}  + Izdr{-  *ssin20o> 


L2tan6  tana  3 
‘ 2 3 1 I 


]“Mm> 

+ | ^R_emdr{a2  + r2}  + Ixdr{sin20o}  + Iz&r{co820o}  j Y I 
„ , , NE 

+ M1L|tan2a1Y  + ~ Vi,i}  + IXdr{*,p,i8in0o} 

1 !.i=1  / ( NE 

+ IZdr{*,E,iCOB0o}J  * T,i[  + J ^V^UfPR.I 
+ [^tana^R-e)^^))^  j +j^^R_emdr{2n^eYtl(a20osin* 

+ (r  + e)a  B coBt)  + $-flY  .(-(r  + e)a  g sinij/  + a20  cosip)]} 

4 o OAji  2 0 2 0 

+ Ixdr{n[-  ^0yti6oCOB20oBin^  - ^^cc^e^os^]}  + IYdr{neo(*0Y 

+ *0X,iCO8^}  + Iz<^r^0o^0Y,icos^0osin^  + ♦ex,icos29ocos'<')}]  % j 
+ j [^R"emdr{-  2fta2(b2  + r0Q)}  + Ixdr{nsin20o}  + Iz<ir{-  fisin20o}]  *0§T 
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X&" 


V .«"*  K.r-"?r 


♦ | - J 

t j ^R'eadr{-  2fl(r20o  + *2(*2V0  * r)>>  + Ixdr{fl(-  2q'EC)Bin0o  + Yoain20o 

- fl0)}  + Iydr{fl0o}  +Izdr{n(-  2<i’FOcos0o  - Yo»in20o  - B0)>]  $ | 

+ 1 CjljtinJjtittjj  | + | (2^1^  + CjL^tan^^y  J 

2«<Ve,1  + 1Vp,i)  - 2nrvlti>  * Ix<ir{n[^,gii(2q,j,oain2eo 

+ Vlneo>  - ♦,F,i8oCOfleo^}  + Vr{fl^,E#i(-  2q,EO8in20o  * ^F08in2eo 

- hBin60]  * ♦,F,i(<l,EO8ln20o  + 2<1'fOCO820o  + 8oCO80o}J} 

+ IZdr{fi^,E,i(<l,F0flln2eo  + 8o6in0o)  + ♦,F.i(~  2<1'fOCO820o 

< he 

- eocoB0o)J>l  i,,  ♦ (♦PPRfl  + [L^ano/CR-e)]^^))^^ 

, NAr  1 

j £[£R~emdr<nZ[*0Y,i(("  VVo  + r)  " r2Bo)sin^ +Cb2(az  - rYQ) 

+ a2rBo)cosi|»)  + b2(a2  - ryQ)  - &2r&Q)3ini>  + (-  b2(a2Y0  + r) 

- r2Bo)coa<()]}J  ^ J + J ^R"emdr{fi2[-  (b2  + a2)0Q  + eb2]}  + Ixdr{ 

- fl260co.ae(>)  * Izdr{02eoco.260>]  ♦,«,{  ♦ j - ^“Vera'l  j 

+ | ^R_emdr{n2[b2(a2  - ryQ ) + a20Q(2r  + e)]}  + Ix<lr{-  Jsn28in20()} 

+ Izdr{J5fi28in260>J  0 | + j K^tanfi^ana^  6 | + j |^R"emdr{fi2[e(r  - a2YQ ) 

- b2r0Q]}  + Izdr{-  fi2Bo8in20o}]Y  j+  j K^tan2^  J + j KyY 
♦Ei1(»260  - e)  - »280»Kfi]>«T,i  j 

* [L2t““l/(R'e):1*ET>i),,T,i  j 


Tall  rotor  control  system  contribution: 


- C1L2tana1X2  - K^tana^ 
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Main  rotor  control  aystem  contribution: 

♦ [-  *sC,L„tana  (-(sind  + coad)ain*i  + (coad  - aind)cos^) 

12  1 Rg 

- L tano  ]X4  + I-  \P  L tana  ((aind  - coad)aln<»  - (cosd 

12  1 A 12  1 *>g 

ind)coaf)  - ^CjL^tana^Jip  + [-  C^tana^^/^Ccoadsin^  + sindcoa^)]XL 


+ a 


+ The  coefficients  of  XA,  X^,  and  X^  which  are  identical  to  those  of 
XA,  Xp,  and  X^,  reapectively , vith  replaced  by  . 

. o (A70) 


144 


• ~ ***$*.  JBGI 


'i»l 

* ♦Y.it*E.J(-  V0**  - - vl,jC0“*  - ♦f..)I!oco,*J  * *z,i*r,j 

* ♦m,1I<*E„)Vo  * ♦f.j'Vo  ♦ r * 8»8ln*  * <-  ♦f.A  - 

- ♦t.j'S  ♦ rt0))co.*]  * ♦6Xfl[(«FiJ(»2  - rt0)  ♦ ♦EiJ(b2  ♦ r60)).in* 

+ l*E,jVo  * ,F,)(V.  * r * e)lco,*:)  * ♦l/.lVl.J8l,,»  * VlVl.j008*’ 

* ^‘♦'F.j^er.i'^V1"*’  * - 8o,ln6o  - i'eo1'05*' 

* ♦ex,i<co'9oco‘*  • <V089o  - 6o‘ln6o  - i'eo1s1”*)J) 

* Iidr(*'E,j[+eY,i',,F0c0'*  - Vi’W10*}*  * Izdr(*,E,jCnr,i(-  8ln6o”ln* 

- (vln90  * eo00,9o  * i'ro,co,*)  * ♦ex,il‘  Bin90COB*  + li,’1"',  ♦ V080,, 

* V^JsinF)])]  \\  *\  [r'"4rfVE,J  * * IxdH<1W’F.j> 


+ IjdrN  1'F0*'Eij ) ♦ ^(i'fO^'e.J1]  *e9T  * M1L29ePH(9Fl>R,J 


V“V*FPR.J 


t • I 

* [L2«M1/(S-.)]tET>J)  eT  j*  j [<"'e«8r<9F,j(»2Y0  ♦ x)  ♦ ♦EiJb2T0 

* Vl,J)  * IX9r^9'F,Jcos9o^  * IZdr<-  ♦E,j8i”6o)]  6 j*  j M,V“V*FPR,j 

* CL2t»n<>1/(R-')>ET,j)  9j  * | [<R-**irl9EJr  * Vl,j'  * IX9r^9'F,JS^n8o) 

* IZ<irf*,E,JCOseo>]  ^ ! + 1 M,l2tlnaiUFPB,j  * [V““l/("-e)]4ET>j)  Y 

(HE  ’ ' . 

* j E [^-e»drU9EiJ  ♦ ♦9f>j)  ♦ IxarUV./V.i1  * Izar{*,E,j9'E,i,J  iT.i 

t*1  NE 

+ (*FPB„)  * [t2tM»./ll,-'l>B,,jl  &<VW 

* p.2tmai/(R-e)>Erjl))qT>1  j ^ [{"''“drt2(1t*8Y,lt9F,jt*2  ‘ 

* * 1 * • * b2Bo,c08*)  * W ♦f,J(Vo  * x ♦ e - b260)sl„* 

* ’ rYo)c08*>^  * IxirW[*'F„)(*OT,i(oo86ocos't'  " (V089o  " I’eb 


ij* 


“ #o,in8o ),ln*)  * ♦ex,i*~  co>e0Bln*  - (Y0«.0B8e  ' q'^,  - eosin6o)co*<»)) 

♦ ♦,E,j(80Y.i(-  ^o001*  + (Vinflo  - 8oCO80o  - <i,E08in2eo)8ln*) 

♦ ♦ex4(,ineo8ln<.  ♦ (Yo«in0o  - flocos0o  - q'E08in26o)co.<») )]} 

+ IY4r^,E,JC00Yfl(8in8ocos*  + Cq'EO8in20()  - 2q'ro«ln280 

♦ 0oco«0o  - Yosin0e)«in*)  + *eXfit-  ■inO^in*  ♦ (q'Eosin20()  - 2q'proin28o 

+ eocoa0o  - yo8ineo)cc  -))  + c°8e0c°8*  + (-  Q’eO0082^ 

♦ I'po81112^  + 0o8ln0o  + Yoco80o)sin*)  + ♦eX,i(CO80o8in<'  ♦ (-  q'E()co820o 

♦ fi'F08in28o  + 80ain®0  ♦ ''oco80o)coa*))]}  + Izdr{n[**p  jU0T  ^cobO^ob* 

“(q,FO8in20o  ■ q,Eoco820o  + VO80o  + ®oflineo)8in*)  + ♦ex,i(-  ®°«e08ln* 

- (q,po8in20o  - q'  jjqCob20o  + YCO80Q  + 6osin0o)co8<») ) + ♦'g  j((-  sinO^OB^i 

+ C2q'F08in2eo  + Y08ine0  + ^cose^sin*)^^  ♦ ♦0Xfi(«in8o8in* 

4 C2q,Fo8in28o  + Yo8in8o  ♦ 8ocob0o)co8<-))]}]  ^ | + j [^R"emdr{2Q(^F>Jb2Bo 

' ♦E,ja20o)  + 20vi,jb2}  + IXdr{O(0,F,jCO89o  + ♦,EJ8in6o)} 

+ IYdr{n(-  ♦'g  jBin0o  + ^,F  jCO80o)}  + Izdr{n( 

" *'E,J8in6o  " 8'f,JCO80o)}]  *0®t|  + j"  CiL280PR^FPR.sJ 

+ r^tanaj/t R-e)]^)  0T  j + j [^'W-  * r0o)} 

+ Ixdr{Q[*'pj(Y0coe8o  - 0osin0o  - q'^)  + ♦,E>JtY0sin80  - 8ocos0o 

- <l,EO8in20o)^}  + IYdr{n^,E,J(<l,EO8in20o  + 2q,FOco820o  “ Vin0o 
+ 8ocob0o)  + ♦•PJ(-  q'BOco820o  + q,F08in2eo  + yoc°80o*  608i.neo)]> 

+ IZdr{QI0,E,J(-  2<1,FOCO820o  - Vin0o  “ 8OCO80O)  + *'FOBin2\ 
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■f  !T  * - 


* i r*  - r-  ' 


+ q,EOCOfl20o  - VO80o  “ 0o8ineo^}]  * ( + J C^L,tan43(^FFR,j 
* CV“,/(B-«»ET,J>  »(+  j [lfR'e"4r<2«Vl0  * >Vf.J  * 

+ Ixdr{n[V*  ,0ftcos0rt  - ,(ertBin0rt  ♦ 2qVBin2eJ]> 


f tJ  0 


E,J  o 


+ IYdr{fl[^EJ(2q*EOBin20o  + ^V1”20*)  ♦ Sosir0o)  ♦ ♦*FJ("  <1,EO8in20o 
- 2<l'FX)co82eo  • 0oCOB0o)J}  + IZdr{fl^,E,j("  q,FX)8in2eo  " 8o,ineo) 


C L tana  (*_„  . 
1 2 lvyFPR,J 


* ♦f,j'2<1’f00O“29o  * ‘»Vln260  + 6oC08V],]';  j8 

* d.3t<“a1/(R-,!»ET,j)J  j |i?1[il,'!"ar{2“[vi,J*E,l  - U1.1*E,J 

* V*v,j*E.l  - ♦f,1*E„))J}  * IXdrtn(80 

+ q'E08in9o^*,K>>Jl,,,E,i  ‘ 

+ ir^Cntq'Eo-^o  - 2q-FOcos0o  ♦ eoX*'EJ*'Ffi  - ♦'Pf^,E,1.» 

+ IZdr{n(q,EOSin0o  - 2q,FOCOB0o  " 8o){“  ♦‘f.J^E.I  + ♦' ’e.J^' *Fti3 1 }K,i I 
| NE  ’ 

+ J Ci^FPRtj  + C^ana^d*-*)]^^)  ^FPR.i 

+ CL2tan0i/(R-e)3*ET,i)qT,i  | + j 2S,JMq,A,J  *T.J  1 
+ j [^R’emdr{n2^Ffj(a2  + ey0]  ’ ♦E,J(0o(r  + e)  + V-1} 

+ Vr{-  02*'E,j(q,EOaln29o  + BoCO80o)}  + Iy^{«2^,E,j(<l,EOKin20o 

+ q,FOCO820o  + 8oCO80o)  + VF,J(-  q,EOcO820o  + q,FOsin26o  + Bo8in0o)]} 

+ Izdr{fl2[-  ♦•E>Jq,FOcos20o  + ♦,FJ(-  q'rosin20o  + q'E()cos  0Q 

- eo8in0o)3>J  ♦0eT  j ♦ j - K,VePR(*FPR,j  + EL2tan«1/(P-e)]^TfJ)eT  j 
+ 1 *b,j*2y0  + ♦f.jCvo  + r + e + 2b2Bo)  ■ vi,jV} 

+ R2^EJ8in0o}  + IYdr{R2(<>,E,j8ineo  - ^P,JCO80O)} 
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V.  f . 


K,L2t““l(*FPB.J 


* IZdr«l2*F,jco,9o»]B  j * | lt,L2t“<3(8FPR1J  * tV“Y<B-E>>ET,j>  6 | 

* j [^-marfU’t-  *E.j<Vo  - - ♦f,J‘2U>]1'  i 9 , 

+ CVanai/(B-e):l»Er,j)Yni^[iR‘'“lrN 

+ Ixlr(»2»lB290rEiJ»-Ejl)  * ^dr(^[4'EjJ(-  ♦'E>18l»2eo  - ht^yinSB^) 

* ♦V,J<^'E>l,In29o  - ♦'F,iC0826o)J)  * Iz^l°2[-  >s*'E,]*'F,l,ln!itl 

* ‘S*'E,l‘ln2eo  *^F.lco828o^,]’T,lj  * j K,f*FPR,j 

* [V"Vt"-*)}*n,j>  UFPR,i  + CV““l/(B‘e):l*EP.i,,1T,ij 


Tail  rotor  control  system  contribution: 

- °l(*FPR,j  * [V“°l/<R-*)]*ET.j)i2 

■ Kl*9FPR,J  * fL2tan“/*B''^*ET,J^X2 
Main  rotor  control  system  contribution: 

♦ [-  j + [i^tano^/fR-e)]*^^^  + BB/R  (~(sin6  + cosfi)sin<|» 

s 

+ (cosfi  - sin6)cosi|i))]XA  + [-  ^Uppp  j + [i^taruWCR-en^^Ml 
+ ^/^((sinfi  - cos6)sin^  - (cos6  + sin5)cosi|>))]Xp  + [-  CjUppR  j 
+ [L2tanai/(R-e)]^£T  j)(cos6sini|>  + sinficos^)]^ 

+ The  coefficients  of  X^,  Xp,  and  X^  which  are  identical  to  those  of  XA 
Xp,  and  X^  respectively,  with  replaced  by  K^. 


(A79) 
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Tall  Rotor  Control  System  Equation! 

m2x2 

s # ME 

+Z  {ClL2^epR®T  “ C1L2tan63$  - C1L2tanaiY  - £ ^C1^FPR  i 
n»l  * 

♦ [L2tana1/(R-a)]4gT^1)4p^1)  ♦ C!  X2)  + C2X2  - C2X3 

N 

(KiL240pR0T  - KiL2tani3$  - KjL2tanaiY 

“ Z (V**TR,i  + [L2tanai/(R-e)]4ETtl)qTti)  + K^}*  K2X2  - K2X3  - 0 

(A0O) 


m3x3  + (c2  + c3  ♦ c31)x3  - c2x2  - c31x31  - c3*4  + (K2  ♦ k3  + K31)X3 
- k2x2  - K31X31  - K3X4-  0 


(A8l) 


M4X4  ♦ (c3  ♦ c4  + c41)x4  - C3X3  - c41x41-  c,fx5  + (k3  + k4  ♦ k41)  x4 
- K3X3  - k31x31  - k4x5  = 0 

(A82) 

• • 

^5X5  + (C4  + C5  ♦ C51  )X5  - C4X4  - C51X51 

♦ (k4  ♦ k5  + k51)x5-  k4X4-  k xa  - K xo=  0 

CA83) 

When  the  servo  is  included,  the  damper  C5  is  not  used. 
m31x31  ♦ C3iX3j-  C31X3  + K3jX33-  K31X3  ■ 0 (A01*) 

*'  a • 

+ ^4lX4i”  C41X4  + K41X4  * 0 CA85 ) 

M51X51  + CsiXsi-CjiXs  + k51x51-  k51x5  - 0 U86) 

The  computer  program  of  the  analysis  allows  the  user  to 
include  or  omit  any  of  the  above  tail  rotor  control  system 
degrees  of  freedom  as  desired. 
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Main  Rotor  Control  System  Equations 

1 1 / / ▼ • t \ > *«  \ *■«#  ri./l  i 


W- + IL)  ♦ "s)]XA  ♦ WV^dj,  - W * * [ /2B|(-  VK 

* ^2R8**eX,lIFA  " ♦eY.iV  + ^,A,*1J  * ^ ^iV1  " 

+ cob5]s1ui(»  + [ain6  - cob6]cobi|<)  - *sC  L tan6  (1  - ^B/n  ([ain6 

vrtt  i 123  K_ 

♦ C086]8im|»  + [ainfi  - cob61cob!|>))£-  *sC  L tana  (l  - ^B/_  ([sin6  + co86]ain^ 

12  1 «g 

NE 

+ [flin6  - cob6]cosi|i))y  - Y ^(♦pPR.i  + (1 

- ^/^([sinfi  + cosfijsin^  + [ain6  - coB6]cosi|i))ir  i)  + (JsCi(RB/R^)2[-(Bin6 
+ coB6)sin#  - \bin6  - cos6)cosi|/  - *5COs26Bin2i{i  - *5sin26cos2^  + *5]+*i;C  )X.} 

1 A 

+ CaXa  + £ {(^C  (*B/p  )2£-  2cos6sini|<  - 2ein6cosi|j  + sin26Bin24>  - cob26cob2iJ>] 

a a n^_  1 «g 

+ *i;C  )X„}  + 2J  teC  (^/o  )2[co86Bin<|i  + sinficos^  + *5(cob26  - sin26)sin2<J> 

1 ^ n«l  1 *S 

- (sin6  + coa6)co868in2^  - (ain6  - coa6)Bin6cos2iJ/]xL  } 

+ Displacement  coefficients  which  are  identical  to  the  velocity  coeffici- 
ents, with  and  CA  replaced  by  and  KA  respectively. 


-kaxoa  = 0 


(A87) 


W '^h  - IFA)  + V^A  + W /^FA  + IL)  + ^ + [ /2R2(-  h^h 
NA  ^ N 

+ ^ S/  /2rs('  *0x»iIfa  " ♦ey.iV  + ^z.iV^  + ^^iV1 

- ([co86  - 8in6]sini|i  + [sin6  + cos6]coai|i))<f>  §„  - %C  L tan6  (1 

Hq  Urn  1 12  3 

- ^/{^([cosfi  - sin6]8inij/  + [sin6  + cos6]cos^))g  - ^C^tano^l 

. NR 

- ^/^([cosa  - 8in6]ainip  + [sin6  + cos6]cosi|>)  )y  - ^ ^Ci^FPR  i 
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♦ [L2t«nai/(R-e)]^ET  ^)(’  - ^/^([coad  - Bind]sin(|i  ♦ [sind  ♦ 

♦ cosd]cofli|i) )q^  ^ + {,cCl(RB/p^)2[-  2cotdBin^  - 2sindcoB^  •*•  Bin2dsin2<> 

N 

- coe2dcos2i|>]  ♦ kC  )X.}  ♦ 53  )2[(aind  - coBd)Bin<>  - (coad 

1 A n-1  1 *8 

+ sind)coBi|i  + 1jcofl26Bin2^  + ^810260082^  + %]  + HiC^Xp}  + CpXp 
N R_ 

+ 53  (^/o  )2[cosdain^  + sindcoB^  - ^(cos2d  •*■  Bin2d)sin2^  - (cob6 

n-1  1 *8 

- sindjcosdsin2^  - (cosd  + sindjaindcos2^]^} 

+ Displacement  coefficients  which  are  identical  to  the  velocity  coeffici- 
ents, with  and  Cp  replaced  by  and  Kp  respectively. 


(A88) 


“ KfX0F  “ 0 

^2R|(-  il'J'x* 

* f ~ IL^XF 

+ £ (C  L ( cosdsinij;  + 8indcosi|0$nDP6m  - CL  tand  ^B/n  (cosdsinij* 

n-1  1 2 *s  0PR  T i 2 3 Rg 

NE 

+ sinficos^jB  _ <^L  tana  ^/^(cosdsin^  + sin6cos^)Y  + ^(-  CjUppg  ± 

♦ [L2tanai/(R-e)]^ET  i) ( cosdsinij/  + sindcoBip))^  i + (^(^/p  )2( cosdsini(» 

+ 8indco8<p  + Js(cos2i5  - sin26)sin2^  - (sind  + cosd)cosdsin2i|;  + (cosd 
- 8ind)sindcoB2^)  )X^  + (*5C^(  ^B/p^)2(  cosdsin^  + sindcos^  - ^s(  cos2d 
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+ ain2$)sin2^  + (sin6  - cos6)cos6sin2^  - (cosfi  + Bln6)sin6cos2i|»)  )X,, 
+ (^C^C^B/-  )2(sin26sin2</  ■*■  2cos2Ssin2i|)  ♦ 2sin26cos2^) )X^}  ♦ 

D 


+ Displacement  coefficients  which  are  identical  to  the  velocity  coeffici- 
ents, with  Cj  and  replaced  by  K and  respectively. 


-VOL  ’ 0 (A89) 

When  the  servos  are  included,  the  dangers  C^,  Cp,  and  are  not  used. 
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+ 


+ 


+ 


->~>r 


Tall  Rotor  Servo  Equation 


(as/y1)(Ap/osC<i)(1  4 Ku2)*0 

(ag/^HV^VV 

HA 

(a8/yi)(Ap/°Bc(l)(  (SFa,X*X,1  + SFA,Y*Y,i  + 8FA,Z*Z 

8FA,0X*0X,i  + SFA,eY*8Y,i^  + ^aS/wi  + K^Xo  " KX 

(°S/Mi^R3/a8^  (8FA,X*X,i  + 8FA,Y*Y,i  * 8FA,Z*Z,i 

8FA,0X*0X,i  + SFA,0Y*0Y,i^ 

NA 

(as/Mi ) (Ri/aS^ ( ^SFV,X*X,i  + SFV,Y*Y,i  + 8FV,Z*Z,i 
8FV,0X*0X,i  + SF\pY*0Y,i^qi^  " ° 


(A90) 


In  this  equation  K represents  or  K^,  etc.,  in  the  tail  rotor  control 

system  shown  in  Figure  A2,  depending  on  which  elements  are  included  in  the 
model  employed  by  the  user.  The  computer  program  of  the  analysis  auto- 
matically incorporates  the  appropriate  value  and  also  allows  inclusion  or 
omission  of  the  servo,  as  desired. 
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Main  Rotor  Servo  Equations 
Ait  Servo  = 


(°S,A/pl,j^AP,A/oB,A  + KAm*,A^0,A 

" *08,A/m1  ,A  ^ *P  ,A /aB  ,ACq , A ^ KAW  2 , AXA 

NA 

“ ^a8,A/yl,A^AP,A/aS,ACqA^  £l  *8PA»X,A*X,i 
* SFA,Z,A*Z,i  * 8FA,0X,A*0X,i  + 8FA,0Y,A*0Y,i 

+ ^a8,A/yl,A  + V^A 


8FA,Y,A*Y,i 


‘kaxa 


NA 


“ ^°8,A^llltA^3,A^aS,A^‘2  ^8FA,X,A*X,i  * 


i«l 


+ 8FA,0X,A*0X,i  + 8FA,0Y,A*0Y,i^ 


NA 


+ ^aS ,A  /Ml,A^RlfA/aS,A^  ^ ^SFV,X,A*X,i 


SFA,Y,A*Y,i  + 8FA,Z,A,*Z,i 
+ 8FV,Y,A*Y,i  * SFV,Z/Z,i 


+ SFV,0X,A*0X,i  + ^XOY.A^Y.i^ 


0 


(A91) 


The  equations  for  the  forward  and  lateral  aervoa  are  identical  to  the 
above,  with  the  subscript  A replaced  by  F and  L respectively. 

The  computer  program  of  the  analysis  allows  the  user  to  omit  or  include 
any  or  all  of  the  servos,  as  desired. 
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Aero dynamic  Force; 


At  this  Juncture,  ve  have  developed  the  differential  equations  that 
govern  the  free  vibration  of  the  dynamic  system.  These  vere  developed 
using  Lagranglan  procedures  and  Lagrange's  equations  in  the  form  of  Equation 
Al.  We  nov  vish  to  add  the  influence  of  aerodynamic  forces  on  the  system 
such  that  Langrange's  equations  become 


Jk  If  j . |£  ♦ |S  ♦ IB 

dt  3qj  3qj  3qj  3q.j 


(A92) 


vhere  the  are  the  generalized  aerodynamic  forces. 

Since  ve  are  interested  in  the  stability  of  the  system  and  not  partic- 
ularly in  forced  response  phenomena,  the  Qj  will  be  those  forces  that  curiae 
from  motions  of  the  dynamic  system.  That  is.  Equation  A92  will  define  a 
set  of  homogeneous  differential  equations  in  the  generalized  coordinates, 
q.  . This  being  a linear  analysis , the  Qj  can  be  obtained  by  the  method 
or  small  perturbations . 


Although  the  analysis  is  linear,  it  is  worth  pointing  out  that  the 
coordinate  perturbations  are  performed  about  certain  steady  initial  value*? 
of  the  generalized  coordinates  that  do  not  necessarily  correspond  to  regions 
of  linear  lift  coefficient,  drag  coefficient,  etc.  Rather,  it  is  assumed 
that  the  perturbations  are  small  enough  for  the  behavior  of  the  system  to 
be  considered  linear  at  that  instant.  ThiB  allows  us  to  include  stall  and 
compressibility  effects  even  though  the  generalized  aerodynamic  forces 
developed  represent  quasi-steady  effects.  In  a later  section,  we  shall 
further  refine  the  aerodynamics  by  adding  circulatory  and  noncirculatora 
unsteady  effects  of  the  types  postulated  by  Theodorsen^'  and  Loevy.'^' 


(1)  Theodorsen,  T.,  General  Theory  of  Aerodynamic  Instability  and  the 
Mechanism  of  Flutter.  NACA  Report  ^96  f National  Advisory  Comml 
for  Aeronautics,  Langley  Aeronautical  Laboratory,  Langley  Field, 
Virginia. 

(2)  Loevy,  R.  G.,  "A  Two-Dimensional  Approximation  to  the  Unsteady  Aero- 
dynamics of  Rotary  WingB,"  Journal  of  the  Aeronautical  Sciences.  Vol, 
2U,  February,  1957. 
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Blade  Elemental  Aerodynamic  Forces 


The  elemental  aerodynamic  lift,  drag,  and  pitching  moment  on  any 
element  1 of  any  blade  n are  given  by 


“n.i 

“ **n,iUn,iCL,n,i6rn,i 

(A93) 

(A9U) 

"n.i 

^n.i  Un,iCM,n,i6rn,i 

(A95) 

Here,  d*^  ^ Is  the  aerodynamic  pitching  moment  at  the  airfoil  aero- 
dynamic center t and  dl^  ^ and  dl^^  act  In  the  directions  shewn  in  Figure 
A5.  The  lift  and  drag  forces  can  be  resolved  into  the  X5 , Y5 , Zg  axis 
system  to  give  elemental  T and  H forces  defined  by 


dT  . » dL  cost  4 - dD  4sint  , 

n,i  n,i  Tn,i  n,i  Tn,l 

dH  4 ■ dL  .sin*  . + dD  . cos*  . 

n njl  Q)X  Q|l  q |i 

where  the  local  inflow  angle  Is  defined 


(A96) 

(A97) 


♦n.i  " tan"I(UPfn,i/UT,n,i: 


(A98) 


of 


Thus,  to  evaluate  dT  dHn  and  dMjj^,  we  must 
Un,i  * %,n,i*  UT,n,i * %n,i’  ^D,n,i*  and  ^M,n,i* 


know  the  values 


cl  n i»  CD  n i»  n,i  depend  on  the  local  angle  of  attack,  and 

this  is* seen  fr&a* Figure  A5  to  be  given  by 


an,i  “ aP,n,i  - ®n,i  “ *n,i 


(A99) 
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Also  from  Figure  A5 , it  is  apparent  that 


UW> 


(A100) 


Up  n i and  UT  n>i,  being  the  local  relative  velocities  in  the  X , Y , Z 
axis  system,  are  obtained  as  follovs:  5 

Neglecting,  for  the  moment,  the  forward  velocity  conponent,  it  is 
readily  .established  by  inspection  of  Figure*  A1  and  using  Equation  A28  to 
define  ^ that  the  absolute  velocity  of  any  blade  element  1 on  any  blade 
n is  giveiS  by 


7n.l  * 7 - *n.l 


(A101) 


where  V - {0,  0,  - VA> 

The  local  relative  velocities  at  any  blade  element  are  obtained  by 
appropriate  transformation  of  Vr  i . Thus 


are 


{ TL  , TL  , - U } = [A-  A^-  A , A A Ac'  Ac'  ]T{V  4) 
T,n,i  R,n,i  P,n,i  q^  qQX  *,n  0,n  Y,n  q*,^  n,i 

(A102) 

Therefore,  in  the  Xg  , Yg , Z5  axis  system,  the  local  relative  velocities 


{UT,n,i*  UR,n,l ' ” UP,n,i}  * * Vn,i 'UP,n,i} 

(A103) 


Inserting  the  steady  initial  displacements  of  the  generalized  coordi- 
nates and  assuming  zero  initial  velocities  of  the  generalized  coordinates  in 
this  expression,  it  can  be  shown  that,  initially, 
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^T,n,i  " ®£®n  + rn,i  + an,iYo,n  “ ^n,i®o,n  ” ^an,i 

” rn,iYo,n^Yo,n  - " VABo,JYo,n  " q'EO,n,i^ 

( AlOU) 

Vn.i  " ‘ fl[(en  + rnfl)q,FO,n,i(Yo,n  ‘ ^EO.n.i''  + (an,i 

" rn,iYo,n^eo,n  + q'pO,n,i^  + VA  CA105) 

Since  the  angle  of  attack  is  measured  at  the  3A  chord,  the  appropriate 
definitions  of  aQ  ^ and  bQ  ^ in  Equations  AlOU  and  A105  are 


Jn,i 


! .A  - EA  )cos6 
n,i  n,i  on,i 

(A106) 

! . A - EA  . )sin0  , 
n,i  n,i  on,i 

(A10T) 

These  quantities  define  the  position  of  the  3/U  chord  on  the  deformed 
blade  relative  to  the  blade  feathering,  or  elastic,  axis . 

Ur  n i is  the  local  relative  radial  velocity  and  is  ignored  since  the 
effects  ’of  radial  flew  are  not  taken  into  account  in  this  analysis. 

We  sure  now  in  a position  to  calculate  the  blade  elemental  aerodynamic 
forces  and,  therefore,  the  generalized  aerodynamic  forces. 

Generalized  Aerodynamic  Forces 


If  we  apply  the  principle  of  virtual  work,  it  can  readily  be  established 
that  the  system  generalized  aerodynamic  forces  cure  given  by 


Vi 


NB 

£ 

i-1 


UdF  .}T 
n,l 


a,i 


{dff 


n,i 


}T  JL 


(A108) 
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J.T- •*■*-*•  T 


for  a la  the  rotating  system,  and 


N NB 

•EE 

n-1  1«1 


UdT  1)T  tt 
n,i 


a<1  <xn  + 

3qn,J  “-1 


{<W  .}T 

n,i 


Vj 


t6n  i» 

n,i 


(A109) 


i 


for  ^ j in  the  fixed  system. 

{dFQ  is  obtained  by  rotating  the  blade  elemental  T and  H forces  such 
that  they  act  in  the  directions  of  the  absolute  displacements.  Thus , 


} 


(A110) 


{dNQ  is  obtained  by  transforming  the  local  blade  elemental  pitching 
moment  rector  to  the  absolute  axis  system.  Thus , 


f 


‘“n.i!  * CVW*.nA6.A.n]f°'dVi01 


CA111) 


{On  i)  is  the  rotational  displacemt  vector  in  the  absolute  axis  system 

on  vhicfi  (dM  .}  acts.  It  Is  given  by 
n9i 


{0n,i}  “ ^ ’'^xVnVnVJ10*  “ 8n,i»0} 


(A112) 


If  ve  nov  perturb  the  generalized  forces  given  by  Equations  A106  and 
A10T  such  that,  in  general. 


I «j  - (All*) 

I 

where  the  (Qj)0  are  the  values  of  the  generalized  forces  at  the  steady 
initial  values  of  the  generalized  coordinates,  and  the  6Q.  are  increments 
of  the  generalized  forces  resulting  from  small  perturbations  of  the  gener- 
alized coordinates,  it  can  be  shown  that  if  we  make  a first-order  Taylor 
expansion  of  the  6Q.  about  the  steady  values  (Qj  )0t  then  the  perturbation 
forces  are  given  by  the  following  general  expression: 
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/ iSi  3T.3Qj.3H 
'3T  3Ut  3H  3Ut 

. 3QJ.  3T  3Q.1  3H 

' 3T  3 Up  3H  3 Up 


3Qj  3m  ,,3Ut.  . 

. iSi  sm  i iiSHLi  «, 

3M  3Up  0 3^ 


+ 


CA11U) 


where  the  summation  is  performed  over  the  number  of  generalized  coordinates. 

The  subscript  o in  this  expression  implies  that  the  quantities  are  to 
be  evaluated  at  the  steady  initial  values  of  the  generalized  coordinates. 

It  will  be  remembered  that  a summation  over  all  blades  is  implicit  when 
the  generalized  coordinates  are  in  the  fixed  system. 

Equations  A108  and  A109  give  the  values  of  the  generalizedjaerodynamic 
forces,  and  they  both  contain  the  absolute  displacement  vector  Xn>i.  When 
^ i as  defined  by  Equation  A28  was  developed,  we  were  at  that  time 
concerned  with  the  displacement  of  the  center  of  gravity  of  the  blade 
element  i.  However,  we  are  now  concerned  with  aerodynamic  forces;  and  if 
we  consider  these  forces  to  act  at  the  element  aerodynamic  center,  then  it 
is  the  displacement  of  this  point  that  must  be  defined.  This  is  easily 
accomplished  by  substituting  ACQ  ^ for  CGn  ^ in  the  expression  for  Xr 

It  now  only  remains  to  define  the  derivatives  given  in  Equation  AllU 
to  establish  the  values  of  the  generalized  aerodynamic  forces  Q.  required 
in  A72 . J 

The  following  list  of  aerodynamic  forces  and  derivatives  is  presented 
without  subscripts  to  preserve  clarity.  This  also  applies  to  the  generalized 
force  expressions,  which  are  presented  m integral  form.  Nevertheless,  it 
should  be  remembered  that  in  these  force  expressions,  all  generalized 

dinates  and  physical  properties  associated  with  the  rotor  system  have 
a subscript  n,  the  blade  number. 
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dT  - %pcuLjlut  - cDUp] 

(A115) 

dH  - JspeU[cLUp  + cDuT] 

(AI16) 

dM  - *3pc2U2Cu 

CA117) 

*1 

- 

^ - !»',el/utoL(u2  + u2)  - c„uTup  ♦ cLi0uTup 

" CD,aUp^  * ‘i(lc  ^v^Cl,MUT  " CD,MUTUP^ 

(A118) 

1 

4- 

S 

= W/Jcjm  - Optu2  ♦ uj)  - cL_au2 

* CuAV>  + !*>c1^Cl,MUTUP  - CD,M^ 

(A119) 

s 

m 

^ - *SPCU^CL,a,JT  - CD,aV 

(A120) 

hl 

- 

^ - **  VWp  * * U}>  * CL,aUP 

* cD,.«r"p]  * !i*1/TrcL,MV’p ♦ 

(A121) 

h 

2 

a 

^ = W/^cJu2  * u|)  * cDuTup  - cL><Aup 

” GD,aUT-^  + ’SPC  ^°L,MUP  + CD,MUTUP^ 

(A122) 

h3 

a 

^ * ^UCCL.«UP  * CD,aUT] 

(A123) 

“l 

a 

3Ut  " ’iPC  ^■GH2UT  + CM,aUP^ 

* ^Vm.mV 

(A12U 

m 

2 

= 

3 Up  ‘ *SPC  ^VUP  * CM,aUT^ 

+ ^pc2l^vCM,M^P^ 

(A125) 

L. 

3 

s 

^ ’ ‘!cc2cM,au2 

(A126) 
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. iHi  + *ut 

3q  3q 


BA 

- tI  (*9I,1(UPco8<'  * UFl"Y0  - I'eO3'1"*1 

‘ ♦ex,i(upIYo  - iVc0‘*  * Vln*))5i 

- (n(»e0  ♦ *•„»♦, *p 

- (il(b  ♦ r60)  ♦ Up(yQ  - I'^ns 

- <“'<%  - «V  * up6oh 

* ? *nt*'E>1(-  a - »Y0  ♦ (r  ♦ a)l'E0)  - ♦e.I^'eO  * *F,lBo 
i=l 

+ rt]  - + 4'^)  + - q'EO^}qT,i 

BA 

+ U0Y,i[b(Yo  ” qV8in'1'  ’ (b  + rVCOB^  + ♦eX,iCb(Yo  - 

+ (b  + rB0)sin<|>]  + ♦x>i[(YQ  - q’^sinij/  - cos*]  + ♦y>i[-(Y0 

- q'B0)cosi|(  - sin*]}^ 

+ {b)^00T 


♦ (b^0  - <i'eo))i5 

* <r  * “«'eo,y 


NE 

+ ^ friq'EO  + *E,i}qT,i 


BA 

t (♦ 

i*l 


3 Up 
34 


V 

( 


* q,F o5810^  - ♦ex,iCup(eo  + q,F0)co8*)>TJ 


+ <«**„  + 4'ro)40)Pr. 

- {fl(a  - nr0)  + up(eo  + q'ro))e 


V008* 


(A127) 
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♦ (fl(r0o  - eq'^b 


HE 

k 

NA 


{«[♦ 


• (_ 
F,i' 


a + 


rvo) 


*E,i^o 


q,po^  " 


V’Ffi(e0 


♦ q* 


FO 


+ £ <("(a  " nr0)coB*  + (r  + e + bcl,F0  ♦ aY0)Bin*)*QY  i 
i«l  * 

♦ ((a  - rr0)ain*  ♦ (r  ♦ e + bq^  ♦ hy0)cos*)*6x  ^ 

+ ♦x.i(#o  + qW8in*  ~ ♦y.i^o  + ^ F0^co8^ 

+ ♦z.i^ 


+ {a>^0dT 

+ {ayo  ♦ r + bq'^jg 
- 

NE 

+ £1  riq'pO  + *F,i>qT,i 

50  ■ If  * If  * {'1)Vt 


(A128) 

(A129) 


Equations  A115  through  A12 9 and  used  in  the  confutation  of  the 
generalized  aerodynamic  forces.  These  forces  comprise  the  right-hand 
side  of  Equation  A92  and  have  the  following  forma: 

( r N NA 
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♦ (ft.*'ro  “ b Yn)flin<;)dT]>^fl0 


0 T 


1 1 

+ {*0X,Jf((blSo  " r ' ai^'EO)8in^  + e(Yo  • 
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Unsteady  Aerodynamics 

Circulatory  Components 

Hie  importance  of  the  effect  of  unsteady  aerodynamics  on  rotor  stability 
has  been  demonstrated  by  a number  of  Investigators.  Anderson  and  Watts (3) 
give  a good  account  of  some  of  the  vork  that  has  been  done.  Their  discussions 
generally  relate  to  flutter  phenomena  experienced  by  unstalled  blades  in 
hover  conditions . Under  such  conditions , the  unsteady  aerodynamic  theories 
postulated  by  Theodors en  and  Loevy  can  be  used.  These  theories  differ  only 
in  that,  if  Theodorsen  functions  are  esq>loyed,  the  effects  of  wakes  from 
previous  blades  either  are  ignored  or  are  negligible,  whereas  the  use  of 
Loevy  functions  permits  the  effects  of  preceding  blades  to  be  taken  into 
account.  Theodorsen  analysis,  therefore,  necessarily  constitutes  a single 
blade  analysis,  while  Loevy  analysis  yields  a multiblade  solution. 

Both  approaches  define  circulatory  unsteady  effects  and  are  readily 
accommodated  by  employing  the  unsteady  aerodynamic  lift  deficiency  function 
C'»  f»  + iG'  defined  by 


C' 


H 


(2) 


1 


H ^(k)  + iH  ^ 
1 o 


(k)  + 2J1(k)W 

(k)  + 2[Ji(k)  + iJ0(k)]w 


(A135) 


This  expression  defines  the  Loevy  lift  deficiency  function.  When 
W*0,  Theodorsen  lift  deficiency  functions  are  defined. 

In  liquation  A135,  the  modifier  W,  which  1b  used  only  when  multiblade 
or  preceding  blade  wake  effects  are  considered,  is  given  by 


N-l  SzS. 

a + I [.khV2”]  " el»n 
a*I 

khN  i2im  , 
e e - 1 


(A136) 


u>c 

where  k=  2U  is  a reduced  frequency  par  Meter,  u>  is  the  flutter  frequency, 
c is  the  chord,  and  U is  the  total  Jnf.lov  velocity . m=w/n  is  the  ratio  of 
the  flutter  frequency  and  the  rotor  speed.  h*>*uUp/cfi  is  the  nondimensional 
vertical  wake  spacing.  $n  i8  the  interblade  phase  which  depends  on  the 
type  of  mode  being  analyzed. 


(3)  Anderson,  W.  D.,  and  Watts,  G.  A.,  Rotor  Blade  Wake  Flutter.  Lockheed 
Report  LR  26213,  December,  1973. 
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For  collective  modes,  ♦n 

For  reactionless  modes , 

For  advancing  modes,  $ 

n 

For  regressing  modes, 

Here,  n ■ 1,  2,  ...  (N-l) 


0 


irn 

(-»/2)n 

(ir/2)n 


Using  all  of  the  foregoing,C'  can  be  defined  at  every  blade  radial 
station . 


In  this  analysis,  it  is  assumed  that  the  rotor  is  in  some  steady-state 
condition  when  the  perturbations  of  the  generalized  coordinates  are  made. 
Therefore,  circulatory  unsteady  aerodynamic  effects  are  incorporated  as 
follows : 


We  have  previously  derived  expressions  for  all  of  the  generalized 
aerodynamic  forces  resulting  from  the  perturbations  of  the  generalized 
coordinates.  Equations  A130  through  A13b.  Clearly,  these  will  exist  whether 
or  not  unsteady  effects  are  present.  Therefore,  any  unsteady  effects  must 
be  added  to  these. 


Hence,  the  unsteady  aerodynamics  are  assumed  to  affect  only  the 
derivatives  of  the  lift  curve  slope  with  respect  to  angle  of  attack  and 
Mach  number,  and  those  aerodynamic  derivatives  associated  with  the  T and 
H forces.  Equations  A118  through  A123.  They  are  assumed  to  have  no  effect 
on  steady  lift,  drag,  or  pitching  moment,  or  pitching  moment  derivatives. 
Equations  All 5 through  All  7 and  A12U  through  A126 . Therefore,  the  circula- 
tory unsteady  components  of  the  generalized  aerodynamic  forces  are  obtained 
by  first  setting  Equations  A115  through  A117  and  A12U  through  A126  to  zero, 
setting  Cl,  Cp,  Cpf0,  and  Cp^  to  zero,  and  multiplying  Equations  A118 
through  A123  by  C'.  The  quantities  given  by  Equations  A115  through  A126 
are  then  substituted  for  their  equivalents  in  the  generalized  aerodynamic 
force  expressions.  Equations  A130  through  A13U. 

Noncirculatorv  Components 

In  addition  to  the  circulatory  effects  Just  discussed,  noncirculatory 
forces  and  moments  exist  that  do  not  depend  on  wake  effects . These  are 
produced  by  air  mass  accelerations  and  can  be  shown  to  be  given  by  expres- 
sions of  the  form 

Sc  • *°V[1;  ♦ < v2  ♦ “>«0J  <“37) 


Sc  " pUX’C1u  Vo  * ((,h  + i)k2  - (li  • “hlik)a0J  (“38) 
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• ' V*'. 


where  h0  corresponds  to  a blade  plunging  coordinate,  defined  positive  down, 
and  a0  Is  a blade  pitching  coordinate,  defined  positive  nose  up.  Lrc  end 
Mflc  are  defined  positive  up  and  nose  up  and  act  at  the  quarter  chord  and 
elastic  axis  respectively.  The  quantity  ahbh  ■ a*  c/?  ■ EA.  ho  is  measured 
at  the  semi-chord. 

These  equations  have  the  form 


“nC 


3lnc 

3a  360 


+ 


3%C  3a 
3a  360  0,0 


(A139) 

(A1U0) 


We  can  now  make  analogies  between  the  coordinates  hQ  and  a0  and  the 
equivalent  quantities  in  the  present  analysis.  ITq  is  the  equivalent  of  a 
change  in  inflow  velocity  and  a0  is  a blade  pitch  coordinate.  Therefore, 
we  can  write 


o 


(Alhl) 


(A1U2) 


(All*3) 


In  addition,  from  equations  A96  and  A97  we  cam  write 


6T  * 5Lcos$ 

(AlUU) 

6H  = $Lsin$ 

(A1U5) 

where  $ is  given  by  Equation  A98. 
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Substituting  Equations  A139  in  Equations  AIM  and  All* 5 . and  using  the 
equivalencies  AlUl  through  All*3,  the  follovlng  noncirculatory  unsteady  T 
and  H force  derivatives  can  be  defined. 


t* 
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t« 
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h' 
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h; 
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3 Up 


3a 


3HNC 
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3HWC 
= 3a 


-ipUT  | nk 

pUUp  | ir(a^k2  + ik) 

-ip Up  | wk 

pUUp  | ff(ahk2  + ik) 


(All*6) 

(A1U7) 

(A1U0) 

(AlU9) 


Using  Equations  All* 2 and  All*3,  the  noncirculatory  unsteady  moment 
derivatives  can  be  shown  to  be  given  by 


ml 


m; 


3MNC 

3UP 

3mNC 

3a 


-ipU  ^ ahTrk 


pU2£  [{ah  + ^)R2  " i(2  ' *h)k> 


(A150) 

(A151) 


The  noncirculatory  components  of  the  generalized  aerodynamic  forces 
are  obtained  by  setting  Equations  A115  through  A126  to  zero  and  substituting 
the  quantities  defined  by  Equations  All*6  through  A151  for  their  equivalents 
in  the  generalized  aerodynamic  force  expressions.  Equations  A130  through 

AI3I*. 


Total  Aerodynamic  Effect 


The  total  generalized  aerodynamic  forces  are  obtained  by  adding  the 
quasi-steady,  the  circulatory  unsteady,  and  the  noncirculatory  unsteady 
components . 
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All  of  the  above  derivations  have  been  dlacueaed  in  general  terms. 
Hovever,  It  should  be  remembered  that  the  aerodynamic  derivatives  and  lift 
deficiency  functions  depend  on  the  local  environment  and  geometry  of  the 
particular  blade  element  under  consideration.  Therefore,  these  quantities 
and  the  parameters  on  which  they  depend  are  functions  of  the  blade  radius. 

Flutter  Solutions 


Clearly,  the  values  of  the  unsteady  quantities  depend,  among  other 
things,  on  the  value  of  u,  the  flutter  frequency.  This  is  normally  an 
assumed  quantity.  If,  for  example,  we  are  employing  Iheodorsen  functions, 
we  would  assume  a flutter  frequency  w and  solve  the  set  of  equations.  If 
one  of  the  roots  in  our  solution  corresponds  identically  with  w,  we  have 
a valid  solution.  If  not,  we  would  try  more  values  of  u>  and  iterate  on  a 
solution.  When  Loevy  functions  are  employed,  only  one  additional  complica- 
tion is  added;  that  is,  in  addition  to  specifying  the  flutter  frequency,  we 
must  also  specify  the  mode  of  flutter,  namely,  reactionless,  collective, 
etc. 


Having  specified  the  mode  of  flutter,  we  proceed  as  described  for 
cases  employing  Theodorsen  functions . This  should  be  repeated  for  all  modes 
of  flutter. 

Coordinate  Transformation 


Hlngeless  or  Articulated  Rotors 

The  dynamic  system  equations  and  generalized  aerodynamic  force  expres- 
sions that  we  have  derived  are  expressed  in  terms  of  the  system  generalized 
coordinates.  The  coefficients  are  seen  to  contain  time-dependent  quantities 
These  could  be  solved  directly  and  the  transient  response  of  the  system 
determined.  However,  transient  responses,  although  useful  frora  the  view- 
point of  forced  response  phenomena,  are  less  satisfactory  for  determining 
whether  or  not  a system  is  stable.  This  ambiguity  derives  from  thr  fact 
that  transient  responses  contain  contributions  from  all  of  the  system  modes 
and  spectral  analysis  may  be  required  to  separate  individual  modal  contri- 
butions. This  is  a costly  and  time-consuming  procedure.  Experience  has 
shown  that  there  are,  in  the  majority  of  cases,  only  a few  modes  that  are 
cause  for  concern,  that  is,  in  regard  to  stability.  Generally,  these  are 
the  fundamental  modes . The  higher  order  modes  sure  normally  important  more 
in  relation  to  forced  response.  Consequently,  in  order  to  obtain  eigen- 
solutions  of  the  system  of  equations  we  have  developed,  other  than  using 
Floquet  type  theory,  it  is  necessary  that  we  remove  the  time  dependence  of 
the  coefficients . This  is  achieved  by  assuming  a form  of  response  of  the 
rotating  system  coordinates  and  making  the  required  coordinate  transforma- 
tions . 

In  this  analysis,  it  was  assumed  that  each  rotating  system  generalized 
coordinate  q could  be  expressed  in  the  form 
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1 " Vo  * Vd(-1)D  * %Vln*n  * Vc'0’*,, 


(*152) 


Thus , each  individual  blade  mode  is  replaced  by  four  rotor  nodes . 

^ defines  collective  notion,  q^  defines  differential  collective  or  reac- 
tionless motion,  and  q6  and  qc  combine  to  define  cyclic  motion,  The  reac- 
tionless motion  as  defined  here  can  only  occur  in  rotors  vith  an  even 
number  of  blades.  The  rotor  warping  motions  involving  terms  in  2ty,  3^» 
etc . , were  neglected . 

To  express  the  differential  equations  of  the  rotating  system  degrees 
of  freedom  in  terms  of  the  rotor  modes,  the  transformation.  Equation  A152, 
and  its  first  and  second  time  derivatives,  Equations  A150  and  A151,  ore 
substituted  directly-  in  Equations  A76  through  A79  and  in  the  generalized 
aerodynamic  force  expressions,  Equations  A131  through  AI3I*.  The  required 
time  derivatives  are 


4 - 1/n40  ♦ VjftjM)11  ♦ 2/N(4a  * nqc)8in0»n  + %(<!<,  + fiq8)cos<j)n 


(A153) 


$ * "'A  + + - nV8in*„  + Vv  + 2flSi 


- n2q  )cosi|;  (A15U) 

C XI 


In  making  these  substitutions , it  is  found  that  the  equations  still 
contain  time  dependent  coefficients.  This  time  dependence  involves  terms 
in  simp,  costp , sin2ip,  cos2\p,  sin2ip,  and  cos2\p,  which  con  be  removed  for 
rotors  that  have  polar  symmetry,  i.e., rotors  with  any  number  of  equally 
spaced  blades  greater  than  two,  by  employing  the  following  identities: 
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- • .t  zr  • * ' 

5,  - - n2^ 

% * 204.  - o2% 

N 

Z cosi|>u  ■ 0 
N 

Z cob  2^i  - 0 

Z cos2^n  - W/2  , N>2 

The  method  of  obtaining  these  vill  be  illustrated  by  example.  If 
equation  A152  is  summed  from  n*l  to  N,  and  it  is  recognized  that  because 

of  the  polar  symmetry  T (-1)°  * Z sini|)  ■ V cosip  » 0,  it  can  be 

. *■*  n n 

seen  that  Y q ■ 1/„  V q ■ q„ . Similarly,  if  Equation  A152  is 

M C.  o 0 

multiplied  by  sin<>n  followed  by  summation  from  n*l  to  N,  then  again  recog- 
nizing that  because  of  the  polar  symmetry  £ sinip  » £ (-l)nsin^  = 

sinT>  cos<»  ■ 0,  it  can  be  seen  that  Y qsinip  = ^/„q  Y sin2ip  . 
n n N l-  ' n IPs  L n 

It  can  be  shown  that  £ sin2^n  ■ /2  for  N>2,  which  leads  to  the  result 

£ qsin<»n  * qfl . Continuing  in  this  manner  using  Equations  All*9,  A150, 
and  Al^l,  the  remaining  terms  are  readily  established. 


£ vi^n 
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It  is  clear  from  the  above  that,  provided  the  polar  symmetry  require- 
ment is  satisfied,  the  time  dependence  of  the  coefficients  can  be  removed. 
This  is  accomplished  in  the  following  manner.  Suppose  we  are  dealing  with 
the  blade  rigid-body  flapping  equations  given  by  Equations  A77  and  A132. 

By  making  appropriate  substitutions  from  Equations  A152,  A153,  and  AI5I+  in 
Equations  A77  and  A132,  the  flapping  equations  will  be  expressed  in  terms 
of  the  rotor  modes  described  earlier.  However,  the  coefficients  will 
contain  time -dependent  terms . This  is  removed  in  exactly  the  same  manner 
as  described  above.  By  summing  the  equations  that  are  expressed  in  terms 
of  rotor  modes  fromn=lto  N,  the  equation  of  motion  of  the  collective  flapping 
(6coi^  *8  defined.  By  multiplying  the  equations  by  (-l)n,  then  summing 
from  n=l  to  N,  the  reactionless  flapping  (Bp)  equation  is  defined.  Multi- 
plying the  equations  by  sinT/n,  then  summing  from  n*l  to  N defines  the  cyclic 
flapping  (63)  equation-, while  multiplying  by  cos^n,  then  summing  from  n=l 
to  N defines  the  cyclic  flapping  (Be)  equation.  This  procedure  is  repeated 
for  all  rotating  system  equations  of  motion. 

The  fixed  system  and  control  system  equations  are  expressed  in  terms 
of  the  rotor  modes  by  substituting  from  Equations  A152,  A153,  end  A15b  in 
Equations  A75»  A80,  A87»  A88,  and  A89,  and  in  the  generalized  aerodynamic 
force  expression  A130.  In  this  case,  the  coefficient  time  dependence  is 
removed  by  simply  summing  each  equation  from  n*l  to  N and  applying  the  above 
identities  as  appropriate. 
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Because  of  the  length  of  the  equations  involved,  no  attempt  was  made 
to  revrite  then  in  their  transformed  form.  Instead,  the  transformations 
were  performed  in  the  computer  program  of  the  analysis. 

Gimbaled  Rotors 


All  of  the  foregoing  can  be  applied  only  to  rotors  that  are  not  gimbaled. 
However,  the  gimbaled  rotor  is  readily  accommodated  since  the  equations  are 
expressed  in  terms  of  rotor  modes.  The  method  employed  is  best  illustrated 
by  example. 

Consider  a four-bladed  gimbaled  rotor  system  whose  blades  can  move 
only  in  an  out-of-plane  direction.  For  zero  spring  restraint  on  the  gimbal 
the  first  mode  of  the  blades  will  be  rigid-body  flapping,  and  since  blade 
motions  are  transmitted  across  the  hub,  the  mode  will  be  purely  cyclic. 
Therefore,  for  this  mode,  the  coordinate  transformation  will  be  as  given 
by  Equation  A152  vith  qQ  and  qn  set  to  zero.  The  next  blade  mode  will  be 
a bending  mode  that  is  the  equivalent  of  the  first  bending  mode  of  a rotating, 
cantilevered  beam.  Tfcis  mode  will  be  purely  collective,  or  reactionlees, 
and  its  coordinate  transformation  will,  therefore,  be  as  given  by  Equation 
A152  vith  <18  and  <lc  set  to  zero.  The  next  mode  will  be  the  equivalent  of 
the  first  bending  mode  of  a hinged,  rotating  beam.  The  mode  will  be  purely 
cyclic,  and  its  coordinate  transformation  will,  therefore,  be  as  given  by 
Equation  A152  with  qQ  and  qD  set  to  zero.  Continuing  in  this  manner,  all 
of  the  rotor  modes  and  their  respective  coordinate  transformations  can  be 
defined. 

An  approach  similar  to  the  above  is  employed  in  this  analysis  co 
describe  the  gimbaled  rotor  dynamics.  Since  the  blade  modes  used  are  fully 
coupled  flatwise/edgewise,  mode  selection  is  slightly  different  but  !-he  logic 
is  the  same.  Mode  selection  and  coordinate  transformations  axe  carried  out 
automatically  in  the  computer  program  of  the  analysis  whenever  gimbaled 
rotors  are  being  studied.  This  is  done  in  the  following  manner. 

The  blade  root  boundary  condition  in  the  plane  of  the  rotor  is  always 
assumed  to  be  cantilevered.  The  gimbal  is  restrained  with  the  spring  kg, 
which  may  of  course  be  zero.  Excluding  torsion,  five  blade  modes  are  auto- 
matically used  which  are,  by  virtue  of  the  coordinate  transformations, 
described  in  terms  of  rotor  modes.  These,  together  with  the  coordinate 
transformations  vised,  are  shown  below. 


Out-of-Plane 

Mode  Root  Biry . Cond,  Type 


Coord.  Trans. 


1 

2 

3 

k 

5 


ke 

Cantilevered 

kg 

kg 

Cantilevered 


Rigid  Body  Flap 
1st  Flatwise 
1st  Flatwise 
2nd  Flatwise 
2nd  Flatwise 


^s » 

*o.  <lD 

^ » Qj) » » % 

‘Is.  <lc 

ck>* 
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Thus , we  have  two  modes  that  are  purely  cyclic,  two  modes  that  are 
purely  collective,  or  reactlonless , and  one  mode  that  has  collective,  reac- 
tionless , and  cyclic  components . It  should  be  remembered  that  the  reaction- 
less components  occur  only  in  rotors  with  an  even  number  of  blades . 


Solution  of  Equations 

To  compute  the  eigenvalues  of  the  second-degree  matrix  system,  it  is 
first  necessary  to  reduce  the  problem  to  the  standard  eigenvalue  form. 

This  is  accomplished  by  multiplying  through  by  the  inverse  of  the  stiffness 
matrix,  and  forming  the  companion  matrix.  This  companion  matrix  is  then 
balanced  to  in?) rove  the  accuracy  in  the  case  vhere  corresponding  rows 
and  columns  have  very  different  norms.  The  balanced  matrix  is  then  reduced 
to  upper  Hessenberg  form  by  stabilized  elementary  similarity  trams  formations , 
and  finally  the  eigenvalues  are  extracted  using  the  modified  LR  method. 

For  a given  eigenvalue,  the  corresponding  eigenvector  is  found  by  using 
the  eigenvalue  to  combine  the  matrices  into  one  NxS  matrix,  thus  yielding 
a redundant  system  of  N simultaneous  equations  in  N unknowns . This  system 
is  solved  by  setting  the  ith  component  of  the  eigenvector  to  1.0  and 
solving  the  reduced  system  obtained  by  placing  the  ith  column  in  the  solution 
vector  position  and  deleting  the  ith  row  by  the  Gauss-Jordon  method  with 
interchanging . 
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APPENDIX  B 


FORWARD-FLIGHT  ANALYSIS 


This  appendix  gives  details  of  a forward- flight  analysis  that  approxi- 
mates forward-flight  effects  by  neglecting  certain  periodic  terms  in  the 
equations  of  motion.  It  will  be  seen  in  what  follows  that  the  analysis  is 
essentially  the  hover  analysis  developed  in  Appendix  A,  appropriately  modi- 
fied to  account  for  both  dynamic  and  aerodynamic  forward-flight  effects . 

Important  assumptions  are  noted  below. 

ASSUMPTIONS 

1.  Aerodynamic  trim  loads  and  derivatives  with  respect  to  perturbations 

in  blade  section  total  velocities  and  angle  of  attack  are  obtained  from 
any  desired  aerodynamic  analysis  (linear  or  nonlinear).  Steady,  first 
and  second  harmonica  of  these  quantities  are  retained. 

2.  The  trim  inflow  velocity  distribution  is  determined  from  any  appropriate 
analysis.  Steady  and  first  harmonics  are  retained. 

3.  The  trim  values  of  the  blade  rigid-body  and  elastic  deflections  and 
slopes  are  determined  from  any  appropriate  analysis . Steady  and  first 
harmonics  are  retained. 

U.  Trim  values  of  steady  and  first  harmonic  pitch  at  each  station  are 
retained. 

5.  Trim  values  of  blade  velocity  as  they  affect  the  trim  dynamic 
forces  acting  on  the  blades  are  neglected  as  higher  order. 

6.  Products  of  trim  quantities  are  generally  neglected  in  the  aerodynamic 
force  terms,  but  they  are  retained  in  the  dynamic  system  equations 
where  deemed  appropriate.  Products  of  trim  quantities  involving  aero- 
dynamic offset  (which  is  assumed  to  be  of  the  same  order  as  the  trim 
deflections)  are  generally  neglected. 

7-  All  products  of  perturbation  quantities  have  been  dropped. 

8.  Terms  involving  products  of  a perturbational  quantity  and  a first-order 
trim  quantity  are  retained. 

9.  The  approach  used  implies  implicit  retention  of  some  'higher  order" 
terms  because  terms  involving  products  of  aerodynamic  derivatives, 
first-order  trim,  and  first-order  perturbation  are  retained.  To  the 
extent  that  the  aerodynamic  derivative  is  a function  of  trim  quantities, 
higher  order  terms  will  be  present. 


177 


DBVELOPMEJrr  or  FOHWARD-FLICBT  AHALY8I8 


In  Appendix  Af  ve  developed  a stability  analysis  for  rotors  in  hover 
or  conditions  of  pure  axial  flow.  The  equations  vere  written  with  all  peri- 
odicity retained,  and  it  was  shown  that  by  making  appropriate  rotating 
system  generalised  coordinate  transformations,  the  periodicity  could  be 
removed,  leading  to  a set  of  constant  coefficient  equations  in  a fixed-axis 
system.  Making  this  coordinate  transformation  necessitated  that  each  rotat- 
ing system  degree  of  freedom  be  assumed  to  respond  in  a prescribed  number 
and  type  of  modes.  This  was  done  to  allow  computation  of  system  eigenvalues 
and  eigenvectors  for  stability  analysis. 

Such  an  approach  is  not  essential;  we  could  have  solved  the  equations 
directly,  obtained  the  system  transient  response  characteristics,  and 
assessed  stability  from  these.  However,  this  approach  is  not  entirely 
satisfactory  since  the  transients  contain  components  of  response  from  all 
modes  and  it  is  not  always  apparent,  without  making  runs  that  involve 
excessive  computer  time,  whether  a system  is  stable  or  unstable.  We  could 
also  have  employed  analytical  techniques  that  determine  eigenvalues  of 
systems  with  time-dependent  coefficients . Floquet  theory  is  one  such 
technique.  Here,  again,  the  time  to  obtain  solutions  has  to  be  considered. 

In  multi- degree-of -freedom  systems  Floquet  solutions  can  be  costly,  and 
their  physical  implications  require  a degree  of  expertise  not  required  in 
solutions  for  constant  coefficient  systems. 

When  forward-flight  effects  are  added  to  the  equations  of  motion  of 
the  system,  it  is  not  possible,  by  a simple  coordinate  transformation,  to 
remove  the  resultant  periodicity  in  the  coefficients.  For  exact  solutions, 
we  have  to  revert  to  either  of  the  alternatives  mentioned  above.  However, 
approximate  solutions  can  be  obtained  without  recourse  to  these  alternatives 
if  certain  higher  harmonic  periodic  terms  are  neglected. 

If  the  same  coordinate  transformation  that  was  used  in  the  hover  analysis 
is  also  used  in  the  forward- flight  equations,  it  is  found  that  in  the  fixed 
system  equations,  a substantial  amount  of  the  periodicity  disappears  by 
virtue  of  the  summation  over  the  number  of  blades,  rotor  system  polar  sym- 
metry being  implied.  However,  in  the  rotating  system  equations,  terms  in 
all  rotor  harmonics  higher  than  the  first  can  appear.  We  do  not,  therefore, 
have  a system  of  constant  coefficient  equations . If  it  is  assumed  that  these 
higher  harmonic  terms  can  be  neglected,  constant  coefficiert 
equations  result  that  are  identical  in  form  to  the  hover  equations  except 
that  they  now  contain  terms  that  sure  functions  of  advance  ratio.  These  can 
then  be  solved  using  standard  eigenvalue  approaches. 
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Hohenemser(^)  suggests  that  the  approach  gives  valid  results  up  to 
advance  ratios  as  high  as  0.5.  Biggers(5)  reached  this  same  conclusion  but 
suggests  that  the  approximation  be  used  with  caution  when  higher  frequency 
modes  are  important . 

Considering  the  conclusions  reached  by  these  researchers  , and  remembering 
Biggero'  caution,  it  was  decided  to  apply  this  approximate  technique  to  the 
stability  equations  derived  in  Appendix  A.  The  approach  used  is  described 
in  vhat  follows . 


By  way  of  example,  let  us  examine  a single  aerodynamic  force  component 
on  a fixed  system  coordinate  q.  arising  from  an  inplane  bending  perturbation 
qg  of  the  blades  of  a rotor  in  forward  flight.  The  condition  is  illustrated 
in  Figure  Bl.  Considering  only  that  part  of  the  force  caused  by  changes 
in  the  blade  tangential  velocity,  it  can  be  seen  that  for  small  perturbations, 
the  moment  in  the  direction  0y  is  given  by 


6M  * 


(^R-e((3-^^-)6Ur(r  + e)sinip)dr) 


(Bl) 


From  Figure  Bl, 


5Ut  = [n(r  + e)  - VpSin^cos  (q'g  + YQ  )3  “ Cn(r  + e)  - VpSiwlO  = V^'g^siml* 

(B2) 


q'E  * Y.  the  moment  becomes 

is  NE 

6M*  51(<t#(Vr  + #ho^'in2*  X UE,iqT,i))dr) 


With 
N 


(B3) 


(U)  Hohenemser,  K.  H.,"Some  Applications  of  the  Method  of  Multiblade 
Coordinates,"  Journal  of  the  American  Helicopter  Society,  Vol.  IT, 
No.  3,  July,  1972. 

(5)  Biggers , J.  B.,  "Some  Approximations  to  the  Flapping  Stability  of 

Helicopter  Rotors,"  AHS/NASA-Ames  Specialists  Meeting  on  Rotoreraft 
Dynamics,  Feb.  13-15,  1971*. 
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The  force  on  the  coordinate 


qj  ie  then  given  by 


«F, 


WM 


i-1 


(Vp(r  + ^0^“fc^Bin2'l'qT,i^^dr^ 


(BU) 


jj  An  examination  of  Equation  BU  reveals  that  only  the  terms  under  the 

]T  are  time  dependent.  If  ve  know  the  time  histories  of  these  terms  for 
each  blade,  wo  could  compute  the  force  component  6Fj . However,  this  is  not 
our  immediate  objective.  We  are,  rather,  more  concerned  with  expressing 
Equation  BU  in  an  analytical  form  suitable  for  inclusion  in  an  eigenanalysis . 

If  the  transients  of  y0  and  vere  available,  we  could  represent  these 

as  Fourier  series  in  harmonics  of  <|».  In  addition,  if  the  response  of  the 
generalized  coordinate  q?  ^ has  known  analytical  form,  it  is  apparent  that 
6Fj  can  be  expressed  analytically . 

In  Appendix  A,  ve  made  a coordinate  transformation  that  expressed 
^ in  the  form 


^,1  " ^N<*TO,i  + /V-1^  ^TD.i  + //NqT8,i8in'1'  + ^N^Tc  .i008’*' 


(B5) 


Suppose  we  now  assume  that  the  Fourier  series  expressions  for  y and 

(b6) 


(~3U^")  can  k*  truncated  such  that 


Y = y + y ain<l»  + y cosib 
'o  '01  'si  ci  r 


and 


1 3Ut  ' 


*01  + t1s8in*  + Sc008*  + ti828in2*  + ^ic2c08^ 


(B7) 


At  this  Juncture,  let  us  also  assume  that  , q^,  q^,  q'pQ,  ^'eo* 

8q  , and  Up  are  all  expressible  in  the  form  of  Equation  B6  and  that  t2 , 
t3,  hj,  h^,  h3,  m , m2,  m3,  dT,  dH,  and  dM  are  all  expressible  in  the  form 
of  Equation  B7. 
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Substituting  Equations  B5,  B6,  and  B7  in  Equation  BU  gives 


j " H * £ *Vr  + ®^,E,i*ey,j 

H i-l 

x Z [(yoj  + Yslsin<»  + Yclcos<|0 


x (t0J  + t simp  + tlccos*  + t-ls2sin2*  + t cos 2ip) 

x ^ ^N^To  ,i  + ^TD.i  + //NqTs,iBin'*'  + V^c.i008^ 

x sin2^]))dr) 


(B8) 


After  expansion  of  the  harmonic  terms  in  thiB  equation  and  appropriate 
summation  over  the  number  of  blades  N,  it  can  be  shown  that 


6FJ  * ^ ( Jl  + e^,E,i^0Y,J^CoqTo,i  + Csq«Pe,i  + CcqTo,i 


NE 


+ terms  involving  products  of  q^,  1#  <LrD>i,  it  and  q^,  t and 

harmonics  of  i p equal  to  or  greater  them  2])) dr)  (B9) 

where  Cq  , Cfl , and  Cc  are  constants  that  erne  functions  of  YQ1,  Y j,  Ycl,  tQl, 
tis’  tic’  tlsi!»  and  tic2‘ 

If  we  assume  that  the  harmonic  terms  in  Equation  B9  can  be  neglected, 
then  performing  the  required  integration  we  can  write 


NE 


SFJ  1 (Wl  * Vt«.1  * Vtc.i’ 


(BIO) 


where  A , Ag  , and  A are  quantities  that  are  independent  of  time.  This 
equation  expresses  ?he  element  of  force  6F,  in  an  identical  form  to  the 
force  elements  in  the  hover  analysis . 

It  will  be  noticed  that  the  removal  of  the  harmonic  terms  in  Equation 
B9  results  in  the  elimination  of  the  reactionless  coordinate  9rpp  ^ from 
the  , or  fixed  system,  equations.  This  was  mathematically  correct  in 
the  hover  analysis  since,  as  the  name  implies,  the  mode  is  indeed  reaction- 
less. It  is  therefore  apparent  that  coupling  of  the  reactionless  mode  with 
the  fixed  system  modes  can  occur  only  if  harmonics  of  4/  greater  than  the 
first  are  retained  in  the  analysis. 
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The  total  force  on  the  coordinate  q<  resulting  from  the  blade  inplane 
bending  perturbation  qE  is  the  sum  of  all  of  the  6Fj  resulting  from  this 
perturbation.  We  derived  above  a single  component  caused  by  changes  in 
tangential  velocity.  Clearly,  other  components  arise  from  changes  in  inflow 
velocity,  changes  in  thrust,  changes  in  drag,  etc.  Also,  additional 
conponents  of  force  result  from  blade  out-of-plane  bending  perturbations, 
and  yet  others  arise  from  blade  inplane  and  out-of-plane  bending  velocity 
perturbations.  In  any  event,  by  following  the  procedure  outlined  above, 
all  of  the  elements  of  force  resulting  from  the  blade  bending  displacement 
perturbations  can  be  expressed  in  the  form  of  Equation  BIO,  while  those 
resulting  from  blade  bending  velocity  perturbations  will  have  exactly 
the  same  form  except  that  they  will  be  in  terms  of  qipo  j,  q^g  and 
qTC  i.Sumning  all  of  the  force  elements  resulting  from’the  dis’placement 
perturbations  gives  the  complete  aerodynamic  coefficients  of  q<j*>  i,  qq*  i, 
and  qTc  ^ in  the  q.  equation.  Likewise,  summing  all  of  the  force  elements 
resulting  from  the  velocity  perturbation  gives  the  complete  aerodynamic 
coefficients  of  q^  it  &pa  and  ,i  in  the  q^  equation. 

Any  coefficient  of  any  displacement,  velocity,  or  acceleration  of  any 
of  the  system  generalized  coordinates  in  the  q<  equation  can  be  obtained 
in  the  same  way.  This  is  so  whether  the  coefficients  arise  from  purely 
dynamic  considerations  or  aerodynamic  considerations,  or  both. 

All  of  the  dynamic  and  aerodynamic  force  elements  in  each  coefficient 
in  the  q.  equation  for  a purely  axial  flow  condition  were  derived  in 
Appendix*5 A.  To  include  the  forward-flight  aerodynamic  effects,  it  is  only 
necessary  to  modify  the  expressions  for  6UT  and  6UP  given  in  Appendix  A 
by  adding  the  quantities  given  below. 


6UT(fwd.  fit.)  = 
+ 


+ 


+ 


NE 

[ (VF(-  VE>i[Y0  - q'E0]sin*  + ♦,E,ic08*)qT,i>] 

[VF(*  eo[Yo  - q'Eo^08^)^ 

[vf(-[yq  - q'E0]sini|>  + cosily] 

NA 

t ^ex,iVF6o^Yo  ” q'EO^qi^ 


(Bll) 
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6Up(fVd.  fit.)  - 

Cl  (V<-  ♦'e>1,'pO  - - ■‘W181-*  * ♦,F,1C08*),1T,1):I 

1*1  ’ * * 

+ [vFcos*e] 


♦ [VF(-  4'ro»in<»  - q’ro[y0  - <TE0]cos*)y] 

NA 

+ [ Z (♦ex.iV^iJ 


(B12) 


Although  ve  have  discussed  only  the  construction  of  the  fixed  system 
equations,  the  procedure  outlined  is  equally  applicable  to  the  construction 
of  the  rotating  system  equations.  The  only  difference  is  that,  whereas  a 
summation  over  the  number  of  blades  appears  in  the  fixed  system  equations , 
as  seen  in  Equation  B8,  no  such  summation  appears  in  the  rotating  system 
equations.  Therefore,  the  force  components  in  the  rotating  system  equa- 
tions will  have  a form  similar  to  that  of  Equation  B8  with  the  ^ omitted. 

Again  by  way  of  example,  a typical  aerodynamic  force  element  in  the 
rigid-body  flapping,  6 , equation  resulting  from  a blade  out-of-plane 
bending  perturbation  is 

NE 

6FB  ^ ^ r*'F,iUP^o  + q,F0^  aUp^^.i^ 

(B13) 


Using  the  assumed  forms  for  the  time-dependent  quantities  as  defined 
earlier,  this  can  be  rewritten 


NE 

% " <£"'*<  ^ (-  ’•♦'f.iHo  * “ps81""'  * 

* uPCcos»)(Boi  + q'pjj  ♦ [Bos  ♦ q’ro.lsln* 

+ C’oc  + q,F0c^COS'l')(t02  + %2sain*  + t2cC08^ 

* Ssa8182*  * t2c2CM2*,(1/»,1Ito,l  * 

* 2/*«T,,l8ln*  + 2/»9Tc.1M8*)J))lr) 
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It  ia  clear  that  if  the  harmonic  terms  in  this  are  expanded,  the  result 
can  be  expressed  In  the  general  form 


a + bsin<i  + ccosip  + dsin2^  + ecos2^  + ...  (B15) 


where  the  a,  b,  c,  etc.,  are  functions  of  the  Up,  8,  q'p,  t,  and  q^. 

If  in  B15  we  ignore  harmonics  higher  than  the  first,  then  after  substi- 
tuting back  into  Equation  BlU  and  performing  the  required  integration,  the 
result  can  be  written 

NE 

6F0  * ^AoiW  + Aoi(-l)  ^TD,i  + As  1 ^Ts ,i  + AciV,i> 

+ ^si^o.i  + ^1^  %),!  + As  2^8,1  + Ac2qTc,i)8in’1' 

+ (Aci«*To.i  + Aci(-1,n«TD,i  + AC2^Ts  ,i  + kcfre  .i)c08^  (Bl6) 

where  the  A's  are  quantities  that  are  independent  of  time.  This  equation 
expresses  the  element  of  force  6F0  in  an  identical  form  to  the  force 
elements  in  the  hover  analysis . ftie  total  force  on  the  rigid-body  flapping 
resulting  from  the  blade  out-of-plane  bending  perturbation  is  the  sum  of 
all  of  the  6Fg  resulting  from  this  perturbation.  Similarly,  the  forces 
from  all  other  perturbations  of  any  generalized  coordinates,  whether  they 
be  displacement,  velocity,  or  acceleration  or  arise  from  purely  dynamic  or 
aerodynamic  considerations,  or  both,  are  obtained  in  the  same  manner  and 
can  be  expressed  in  the  same  form  as  Equation  Bl6. 

If  each  of  the  rotating  system  equations,  expressed  in  the  form  of 
Equation  Bl6,  is  multiplied  by  (-l)n  and  summed  over  the  number  of  blades 
N,  the  reactionless  mode  equations  are  obtained.  Summing  each  equation 
over  the  number  of  blades  gives  the  collective  mode  equations . Multiplying 
each  equation  by  sinij*  and  cos^,  respectively , and  summing  over  the  number 
of  blades  gives  the  cyclic  mode  equations . On  completing  these  manipu- 
lations, it  Is  again  found  that  the  reactionless  mode  is  completely  decoupled 
due  to  the  harmonics  of  ♦ above  the  first  being  discarded. 

In  some  of  the  system  equations , products  involving  sin0o  and  cosBQ 
are  encountered.  In  this  analysis,  the  blade  angle  90  is  not  considered 
to  be  small . Therefore,  when  substitutions  of  the  form  given  in  Equation 
B 6 are  made  for  0Q,  we  get 
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(BIT) 


Bin0. 

0 

CO*0 

0 


sin(0Oi  ♦ 0B1Bin*  + 0cjcob*) 
coa(0Qi  ♦ eSl8l-n*  ♦ 6c1cob*) 


These  were  expressed  in  terma  of  Bessel  functions  such  that 

sln0  ■ B 8in0  + B cob 6 
o l ol  2 01 

• CJ0(^l)J0(eci)  + (2(J2(^l,J0(eci) 

’ J2{eci)Jo(esi)))co82* 

- 2J  ( 0 )J  (0„  )sin2|h 

1 81  1 ci 

- UJ  (0_  )J  (0  )co822^]sin0 

2 8 1 2 c i Oi 

+ [2Ji(0Si)Jo(0ci)sin<, 

- UJ  0fl  ^)J^(  ec  i)sim|cos2^ 

+ 2J1(0C1)JO(9S1)CO8* 

+ la1(9cl)J2(eaj)coa^os2'(!]cose01 

cose  * B.ccsO  , - B„sin0  , (Bl8) 

0101201 

The  Bessel  aeries  was  trucated  at  J.  due  to  the  rapid  decrease  in 
the  magnitudes  of  the  ltmctions  beyond  this  point  at  typical  values  of 
the  arguments . 

It  should  be  clear  from  what  has  preceded  that  literal  expression 
of  the  forward- flight  equations  would  indeed  be  a monumental  task.  There- 
fore, this  task  was  left  to  the  computer.  Special  routines  wore  developed 
that  perform  all  of  the  manipulations  described  here,  the  end  product 
always  being  the  stability  matrices  from  which  the  system  eigenvalues  and 
eigenvectors  are  extracted. 


186 


LIST  OF  SYMBOLS 


- qgQ  ♦ (%  - EA)co.0o 

- qEQ  - ACcos0o 

- qEQ  - CGcob0o 

- EA/b^ 

Servo  piston  area 

Distance  from  blade  elastic  axis  to  aerodynamic  center — 
positive  toward  leading  edge 


Non-time  dependent  coefficients  arising  from  expansions  of 
products  of  harmonic  series  in  forward  flight  analysis 


Hub  roll  transformation  matrix 

Hub  pitch  transformation  matrix 

Blade  azimuth  tranr formation  matrix 

Blade  flapping  transformation  matrix 

Blade  lagging  transformation  matrix 

Blade  pitch  transformation  matrix 

Blade  inplane  bending  slope  transformation  matrix 

Blade  out-of-plane  bending  slope  transformation  matrix 

" “TO  * <%  - “)sirtlo 
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1 

*>2 

c 

C' 

C1 


2 * 3 *4 * S 


■ «F0  * AC“ln8o 
’ 1TO  - C0,ln8o 
Chord  /2 

Chord 

Lift  deficiency  function 

Damping  constant  at  pushrod 

Tail  rotor  control  system  damping  constants 


C Tail  rotor  control  system  damping  constants 

31*41,51 

C^  Lift  coefficient 

CD  Drag  coefficient 

CM  Pitching  moment  coefficient 

Cp  Damping  constant  at  forvard  main  rotor  servo 

CA  Damping  constant  at  aft  main  rotor  servo 

CTT  Dancing  constant  at  lateral  main  rotor  servo 

LL 

C Servo  valve  flow  gain 

Cp  Servo  valve  pressure  gain 


C. 

C 


'IS 


*L,C 


'*o 


"M 


•a 


'L,M 


"D,M 


Non-time  dependent  coefficients  arising  from  expansions  of 
produc+s  of  harmonic  series  in  forvard  flight  analysis 


Servo  leakage  coefficient 

Slope  of  lift  curve 

Slope  of  drag  curve 

Slope  of  pitching  moment  curve 

Partial  derivative  of  lift  coefficient  vith  respect  to  Mach 
number 

Partial  derivative  of  drag  coefficient  vith  respect  to  Mach 
number 
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V 


- ,-L. 


1 

i 

cu  M Partial  derivative  of  pitching  moment  coefficient  with  respect 
’ to  Mach  number 

CO  Distance  from  blade  elastic  axis  to  center  of  gravity — positive 

tovard  leading  edge 

dl  Elemental  lift 

dJD  Elemental  drag 

I 

dM  Elemental  pitching  moment 

dT  Elemental  thrust 

dU  Elemental  inplane  force 

{dK}  Blade  elemental  lift  and  drag  force  vector  in  absolute 
coordinate  system 

{dM}  Blade  elemental  pitching  moment  vector  in  absolute  coordinate 
system 

D Dissipation  potential 

e Blade  offset 

{E}  Blade  offBet  vector 

EA  Distance  from  blade  semi-chord  to  elastic  axis — positive  tovard 

trailing  edge 

F'  Real  part  of  lift  deficiency  function 

G Blade  torsional  modulus  of  elasticity 

G'  Imaginary  part  of  lift  deficiency  function 

Hankel  function  of  second  type  of  order 
r 

h Wake  spacing  parameter 

h Blade  plunging  coordinate  used  in  definition  of  unsteady 

° aerodynamics 

h Partial  derivative  of  drag  with  respect  to  local  blade  taugen- 

1 tial  velocity 

h Partial  derivative  of  drag  vith  respect  to  local  blade  vertical 

2 velocity 
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h 

3 


[I] 


Partial  derivative  of  drag  vith  respect  to  local  blade  angle 
of  attack 

Blade  flatvise,  edgewise,  torsional  mass  moment  of  inertia 
matrix 

Blade  torsional  mass  moment  of  inertia 

Blade  edgewise  second  moment  of  area 

Elemental  blade  flatvise  mass  moment  of  inertia 

Elemental  blade  torsional  mass  moment  of  inertia 

Elemental  blade  chordvisi!  mass  moment  of  inertia 

Total  blade  torsional  mass  moment  of  inertia 

Blade  mass  moment  of  inertia  about  lag  hinge 

Main  rotor  svash  plate  fore/aft  mass  moment  of  inertia 

Main  rotor  svash  plate  lateral  masB  moment  of  inertia 

Local  blade  polar  second  moment  of  area 

Bessel  function  of  first  kind  of  order 


r 


K 


Ke 

K 

1 


K 

2 » 3 ,4 , 5 
jr 

31,41,51 

*ma 

kf 

ka 


Blade  torsional  stiffness 
Airframe  mode  generalized  stiffness 
Blade  lag  hinge  spring  rate 
Blude  flapping  hinge  spring  rate 

Pitch  beam  stiffness  or  stiffness  at  main  or  tail  rotor  blade 
pushrod 

Tail  rotor  control  system  spring  rates 
Tail  rotor  control  system  spring  rates 

Stiffness  of  tail  rotor  pitch  actuator  for  pure  moment  applied 
at  pitch  beam  end 

Spring  rate  at  forward  main  rotor  servo 
Spring  rate  at  aft  main  rotor  servo 
Spring  rate  at  lateral  main  rotor  servo 
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Spring  rate  at  servo  support 
Reduced  frequency  parameter 

Distance  from  blade  elastic  axis  to  pushrod-- positive  toward 
leading  edge 

Radial  location  of  blade  pitch  horn 

Length  of  one  arm  of  tail  rotor  blade  pitch  spider  beam 

Noncirculatory  lift 

Length  of  tail  rotor  pitch  beam  arm 

Blade  elemental  mass  and  flutter  frequency  ratio 

Airframe  mode  generalized  mass 

Noncirculatory  pitching  moment 

Generalized  mass  of  blade  bending  modes 

Generalized  mass  of  fixed  system  modes 

Mass  at  pushrod 

Tail  rotor  control  system  masses 

Tail  rotor  control  system  masses 
51 

Main  rotor  swash  plate  mass 
Mach  number 

Partial  derivative  of  pitching  moment  with  respect  to  local 
blade  tangential  velocity 

Partial  derivative  of  pitching  moment  with  respect  to  local 
blade  vertical  velocity 

Partial  derivative  of  pitching  moment  with  respect  to  local 
blade  angle  of  attack 

Blade  number 

Number  of  blades 

Number  of  blade  bending  modes 

Number  of  fixed  system  modes 


Ns 

’t.J 

qF0 

qE0 

qeY 

qx 

qo 

qD 

«. 

qc 

J 

f 

^8  ( 

) 

qF0'  J 

qPX)s( 
qFOe  ' 

Q 

^6T 

“6 

S 


Bulk  modulus  of  servo  fluid 

Fixed  system  mode  generalized  coordinate 

Blade  bending  mode  generalized  coordinate— bending  up  and 
leading  positive 

Steady  blade  flatwise  deflection— up  positive 

Steady  blade  inplane  deflection — lag  positive 

Hub  pitch  coordinate 

Hub  roll  coordinate 

Hub  lateral  coordinate 

Hub  longitudinal  coordinate 

Hub  vertical  coordinate 

Collective  mode  coordinate 

Reactionless  mode  coordinate 

Sine  cyclic  coordinate 

Cosine  cyclic  coordinate 

Blade  bending  collective,  reactionless,  Bine  cyclic,  and 
cosine  cyclic  coordinates 


Blade  flatwise  bending  slope  steady,  and  sine  and  cosine 
coefficient  components 

Generalized  forces 

Generalized  force  on  J 'th  fixed  system  mode 
Generalized  force  on  blade  pitch 
Generalized  force  on  blade  flapping 
Generalized  force  on  blade  lagging 
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J3  oS° 


Generalized  force  on  J »th  blade  bending  mode 
Flow  through  servo  valve 


Q Generalized  aerodynamic  force 

J 

r Radius  of  local  blade  element  from  offset. 


rj  Radial  location  of  inner  snubber  of  crossbeam  rotor 

r Radial  location  of  outer  snubber  of  crossbeam  rotor 

i 

r Rotor  radius 


Radius  to  servo  connections  on  main  rotor  swash  plate 
Radius  to  pushrod  connections  on  main  rotor  swash  plate 


r Servo  linkage  ratios 

1,2,3 

SF  Servo  valve  feedback  factors 

v,x,Y,z,ex,ey 

SF  Servo  support  feedback  factors 

t ,X,Y,Z,0x,0y 

T Kinetic  energy 


t Partial  derivative  of  thrust  with  respect  to  local  blade 

1 tangential  velocity 


Partial  derivative  of  thrust  with  respect  to  local  blade 
vertical  velocity 


t 

3 


t 


is 


t 

t 


lc 
18  2 


t 


lc  2 


t 


02 


t 


28 


2C2 


U 


Partial  derivative  of  thrust  with  respect  to  local  blade  single 
of  attack 


Steady,  first,  and  second  harmonic  coefficients  of  the  peurtieil 
of  thrust  with  respect  to  tangentisil  velocity 


Steady,  first,  and  second  hsumonic  coefficients  of  the  partial 
of  thrust  with  respect  to  inflow  velocity 


Total  local  blade  inflow  velocity 
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UP0  ) 

Ups 

| 


V 


X 

2 3 4 5 

» 9 » 

X31  41  51 
» » 


Local  blade  vertical  velocity 

Steady  and  firet  harmonic  coefficients  of  inflow  velocity 

Local  blade  tangential  velocity 
Speed  of  sound 
Potential  energy 

Total  volume  of  servo  fluid  under  compression 
Rotor  axial  velocity 
Forward- flight  speed 
Airspeed  vector 

Loevy  unsteady  aerodynamics  factor — Equation  AI36 

Absolute  displacement  vector 

Displacement  at  blade  pushrod 

Tail  rotor  control  system  motions 

Tail  rotor  control  system  motions 


*F 


XA 

*L 


Motion  of  main  rotor  swash  plate  *‘t  forward  servo 
Motion  of  main  rotor  swash  plate  at  aft  servo 
Motion  of  main  rotor  swash,  plate  at  lateral  servo 
Motion  of  output  of  forward  main  rotor  servo 
Motion  of  output  of  aft  main  rotor  servo 
Motion  of  output  of  lateral  main  rotor  servo 
Feedback  motion  at  servo  valve 
Feedback  motion  at  servo  support 
Motion  at  output  of  tail  rotor  servo 
Motion  at  servo  support 
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i 


a 

a 

o 


a 

a 


1 

8 


®col  \ 
®n  [ 
0s  ( 
Bc  ) 

B°1 

B08 

B0C  ) 


Y 

Yo 

Yoi 

YS1 

YC1 

5 

3 

fiP 

6fb 


PR 


fir 

fiF 

5M 

6F6 


Local  blade  angle  of  attack 

Blade  pitch  coordinate  used  in  definition  of  unsteady  aerc dy- 
namics 

Pitch- lag  coupling — lead,  pitch-up  positive 
Servo  linkage  ratio 

Blade  rigid-body  flapping  generalized  coordinate — up  positive 
Steady  blade  coning— up  positive 

Blade  flapping  collective,  reactionless , sine  cyclic,  and 
cosine  cyclic  coordinates 


Steady  and  first  harmonic  coefficients  of  blade  flapping 

Blade  rigid-body  lag  generalized  coordinate — lead  positive 
Steady  blade  lag — lead  positive 

Steady  and  first  harmonic  coefficients  of  blade  lagging 

Pitch-flap  coupling — flap  up,  pitch-down  positive 
Servo  differential  pressure 

Angle  between  airframe  longitudinal  axis  and  main  rotor  fore/ 
aft  servo  axis — yaw  right  positive 

Angle  measured  at  center  of  rotation  from  blade  feathering 
axis  to  pushrod — positive  in  direction  of  rotation 

Element  of  blade  radius 

Perturbation  force 

Perturbation  moment 

Perturbation  force  on  flapping  coordinate 


I 
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Change  in  inflow  velocity  resulting  from  displacement  and 
•velocity  perturbations 


( Change  in  tangential  velocity  resulting  from  displacement 

T and  velocity  perturbations 


Motion  of  servo  valve  relative  to  housing 

Fraction  of  critical  structural  damping  of  blade  bending 
inodes— -based  on  modal  frequency 

Fraction  of  critical  structural  damping  of  blade  pitch  mode — 
based  on  rotor  speed 

Fraction  of  critical  rigid-body  lag  damping— based  on  uncoupled 
lag  frequency 

Fraction  of  critical  structural  damping  of  fixed  system  modes — 
based  on  modal  frequency 

Blade  pitch  generalized  coordinate — leading  edge  down  positive 
Steady  blade  pitch  angle — leading  edge  down  positive 
Blade  pitch  normal  coordinate 

Geometric  blade  pitch  angle — leading  edge  up  positive 


Steady  and  first  harmonic  coefficients  of  blade  pitching 


Blade  flatwise,  torsional,  and  edgewise  angular  velocities 


Absolute  blade  angular  velocity  vector 
Servo  parameter — Equation  AbT 
Servo  parameter — Equation  A68 

= (f  *,E,iq'EO  + <>,F,iq,FO^dr^  “ “ri 

° ^ ^E.i^E.j  + *,F#i*,F,J^dr^  = "ri 
Air  mass  density 
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* 

*X,Y,Z 

*ex,eY 

♦p 

^E 


p0PR 


pFPR 


'ET 


n 


“ti 

UT2 

“p 

“y 

“en 

(1) 

q 

“y 

“a 

01 


n 

Subscripts 

A 


F 

L 

i 


Local  blade  inflow  angle 

Fixed  system  translational  mode  shapes  at  hub 
Fixed  system  rotational  mode  shapes  at  hub 
Blade  flatwise  bending  mode  shape 
Blade  inplane  bending  mode  shape 
Blade  torsional  mode  shape 

Blade  torsional  mode  shape  at  pushrod  radial  station 

Blade  flatwise  bending  mode  shape  at  pushrod  radial  station 

Blade  inplane  bending  mode  shape  at  tip  radius 

Wake  phase  angle  parameter 

Blade  azimuthal  angle 

Blade  asymmetric  torsional  frequency 

Blade  symmetric  torsional  frequency 

Hub  pitch  frequency 

Hub  yaw  frequency 

Blade  edgewise  natural  frequency 

Frequency  of  blade  bending  modes 

Uncoupled  rigid-body  lag  frequency 

Frequency  of  fixed  system  modes 

Flutter  frequency 

Rotor  speed 

Refers  to  aft  servo 

Refers  to  forward  servo 

Refers  to  lateral  servo 

Refers  to  blade  element  or  mode  number 
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J 


Refer*  to  node  or  force  number 


Refers  to  mode  number 

Refers  to  blade  number 
Refers  to  pushrod 

Refers  to  hub  lateral,  longitudinal,  and  vertical  directions 
Refers  to  hub  pitch  and  roll  directions 

Superscripts 

c Means  coupled 

u Means  uncoupled 

Differential  Notation 

• Differentiation  with  respect  to  radius 

Differentiation  with  respect  to  time 

•• 

Second  differential  with  respect  to  time 
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